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THE MOMENT ON A SUBMERGED SOLID OF 
REVOLUTION MOVING HORIZONTALLY 


By T. H. HAVELOCK (King’s College, Newcastle upon Tyne) 
[Received 20 February 1951] 


SUMMARY 


The moment, due to surface waves, on a submerged solid of revolution moving 
axially at constant depth below the surface of the water is examined in detail. 


1, A SUBMERGED solid of revolution moves axially with uniform velocity 
and with its axis at a constant depth below the surface of the water. If the 
solid is such that the motion in an infinite liquid can be represented by a 
known source-sink distribution along the axis, the horizontal and vertical 
forces on the solid due to the wave motion can readily be obtained to the 
usual approximation; however, for the moment about a transverse hori- 
zontal axis it is necessary to obtain the velocity potential to a higher degree 
of approximation, a point which was noticed in an early paper on the 
circular cylinder (1) but which has sometimes been overlooked. In the 
present note we consider a prolate spheroid, for which this extension can 
be carried out; the form of the additional term in the moment in 
this case suggests an approximation applicable to other elongated solids 
of revolution, such as a Rankine ovoid, generated by an axial source 
distribution. 


2. We suppose the spheroid to be held at rest in a uniform stream of 
velocity c in the negative direction of Ox, the axis being at a depth f below 
the free surface of the water. We take O at the centre of the spheroid, 
Oz along the axis, Oy transversely, and Oz vertically upwards. Using the 
known axial distribution for motion in an infinite liquid, the velocity 
potential is given by 

( k dk 
1 (a k)e's 


‘-+-K,sec"6 
KT Ko§ e-K(2f—2)+ikw dye 


k dk [ dé [ 


kK —Kk, Sec? +-ipu sec 6 


T 0 


e*) log{(1 e)/(] e)}, (a—k)cos 6 t y sin 6, 


ky = g/c?, and the limit is taken as p> 0. 
[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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The first two terms in (1) satisfy the condition at the surface of the solid, 
The third term, which we shall denote by ¢y, is the first approximation to 
the wave motion; its form is determined so as to ensure that the three 
terms together satisfy the condition at the free surface of the water (2). 

The second term in (1), which is the velocity potential for the axial 
motion of a prolate spheroid, is usually given (3) as 2AcaeP,(u)Q,(f) in 
terms of coordinates specified by x = aeul, y = ae(1—p?)!(f2—1)! cosa, 
z = ae(1—p?)*(f?—1)! sinw; it can readily be verified that the two forms 
are equivalent. This equivalence is a particular case of a general relation 
which does not seem to have been stated explicitly, and the opportunity 
is taken of recording it here in view of its use in problems dealing with the 
motion of a spheroid. The relation expresses prolate spheroidal harmonics 
in terms of axial distributions of poles or multi-poles. Using the appropriate 
form of the known general expansion of reciprocal distance (4), it follows 
at once that a 
F P, (k/ae) dk 


2 J f{y?+224+(a—k)2 
—ae 


P,(u)@,,(0) 


For the general case, forming the corresponding expansion for the potential 
of a multi-pole, it can be shown that 


ae 


P8 (u)Q5 (Lei 4( 1»( 0s ) | (a2e2— k2)§s P8 (k/ae dk 
c 


Y cz fy? +224 (x—k)?\4 


ae 
We use the theorem that the forces on the solid can be obtained as the 
resultant of forces on the internal sources, the force on a typical source m 
being —4zpmq, where q is the fluid velocity at the point other than that 
due to the source itself; in fact, we may omit the part of the velocity due 
to all the other internal sources and sinks. Thus for the horizontal force, 
or wave resistance R, we have 
es 
R 4p Acx(éd,/éx) dx, (2) 
—ae 
taken along the axis y = z = 0. 
Taking ¢, from (1) and omitting terms which, on account of the integra- 
tions in x and k, give no contribution to the final result, this reduces to 





ae ae tar x 
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R 16px2c2A2 | a da k dk | sec30 dé [ as dk, 
} : K—Kk,sec?0+-ip sec 6 
ae ae 0 0 


where the imaginary part is to be taken. 
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The integration in x may be transformed in the usual manner by treating 
«as a complex variable and integrating round a suitable quadrant accord- 
ing as x—k or 0. Finally we obtain 

R 32779 pare 342 sec J? (eg ae sec O)¢ 2ko f sec? dé. (4) 
which is the known expression for the wave resistance. 

The vertical force, Y, apart from that due to buoyancy, can be obtained 
similarly from 


ae 


r 477 p ( Acx(éd3/éz) dx. (5) 


This involves the real part of the contour integrals in « referred to above, 
ind leads to double integrals; the expression for Y can easily be written 


down, but it is not very suitable for numerical computation. 


3. The moment of the forces about Oy requires more consideration, and 
we shall take it in two parts. 

We calculate first the moment on the initial source distribution arising 
from the vertical component of the velocity derived from the term ¢3 
Thus 


in (1); we denote this part by G;. 


G, tp | Acx?(éd,/ez) dx, (6) 
taken along the axis. 

3ut we have to proceed to a further approximation to the velocity 
potential, because the uniform stream produces on this second approxi- 
nation a contribution to the moment of the same order as G,; we denote 
this second part by G,. Let 4, be the term to be added to (1) for the next 
pproximation. This term represents some distribution of sources and sinks 


within the spheroid; if M is the total moment of this distribution resolved 


parallel to Oz, then we have 
Gs trpcM, (7) 


ind the total moment on the solid to this stage is G,+-G,. 


4. From (6) and (1), we have 


K(k-+K,Sec"@) 





K—K,sec?0+-in sec 8” 


x 2«f+ix(e—k)cosé dic. (8) 
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Treating the integration in « as before, and carrying out the integrations | 
in x and k, this leads, after some reduction, to 


hor 


P sin p 
sec (=? — cos r) x 


Pp 


G, 647 pare?A?(c?/K) 
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0 
9 Sin Pp 


et f sec?@ dé, (9) 
p 


(psin p+2 cos p— 


with p = x,aesec 0. 
We now determine the next approximation to ¢ so as to satisfy the 


condition at the surface of the spheroid, namely that the normal component * 


of velocity from ¢,+-¢, must be zero over the spheroid. 
We use coordinates py, f,w given by 
x aepl, y ae(1—p?)'(C?—1)sinw, z ae(1—.?)#(C2— 1)! cosw, 


(10) 
the spheroid being given by ¢ = {, = IL/e. 
If, in the neighbourhood of the spheroid, 4, is expressed in the form 


$3 


Ms 


¥ (. 4* cos sw+ Bs sin sw) P8() P(C), (11) 


r=] s=0 
then the required next approximation is given by 


ri 


= <7 pe: Q"( (A; cos sw + By sin sw) Pr(u)Q;(6). (12) 


By considering the behaviour of the terms in ¢, as {> ©0, we see that 
the only one which contributes to the moment M referred to in (7) is the 
term in P}(1)Q}(f)cos w; this latter quantity approximates to — 2a*e?z/3r° 
as €-> 0. Alternatively, we may get the same result from the expression 
of this term as a line distribution of doublets, parallel to Oz along the axis 
of the sphecola between the two foci. Hence, putting in the value of the 
factor P}'(f,)/Q}'(f,), we have 


M = ja’e?BA}, (13) 
with B-1 4 log{(1+e)/(1—e)}+-e(2e?— 1)/(1—e?). 


To determine A} we take from the expression for 4, in (1) the term in the 
integrand involving the coordinates, namely 


exp «(z+ 7x cos 6-+- ty sin 8), 


and expand the value of this on the spheroid in the form 


— 
SY > (Cicossw+ D§ sin sw) P§(p). 
r=0 s=1 
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The coefficient C} is then given by 
Arr : 4 
= C} | (l—p?)* du | exp{«b(1—,?)*(cos w+-7 sin w sin 6) + 
1 0 


+ixap cos 6}cos w dw. (15) 
The integrations can be reduced to known forms, and we obtain 
Cc} 3(7/2«a%e*)'b sec!@ J; (Kae cos 8). (16) 
Hence, from (11), the corresponding term in the integral for A} is 
3(7/2xae)* sec!@ J; («ae cos 6). (17) 


Using (7) and (17) in (1) we obtain the expression for G,, which may be 
written in the form 


ae 977 x 
am & . ‘ * | (sin 4 oer ys 
(i, = 16paec?AB | kdk | sec6 dé | D4 — cos q\e-2e/—ixk €088 dic, 
: , ‘ ‘ q 
ae 0 0 
(18) 
where the real part is to be taken, and 
gq = Kae cos 8, D («+ x, sec?6)/K(K— Kk, sec76-+- iu sec A). 
\fter carrying out the integrations in « and k this leads to 
as . * (sin p . on sad 
G, 647 pa*e2A B(c?/K,) P_ cos p} e~2xe/ see“ sec A dO, (19) 
: . pP 
0 
with p K, ae sec 0. 
For computation it is convenient to express these results in terms of the 
-called spherical Bessel functions, of which tables are available. If we 
write ’ 
S,43(p) (77/2p)' J, 43(p), 
I, } S (p US (pye 2kof sec? sec’é dé, 
0 
I, | S?/ pye 2xo f sec*@ sec? dé. 
0 
e have : 2 
R 647g pK,a*e ‘Ai, 
G, 647g pare 4 A?( Ke, ae i, 2I,), 
G, 64rgpatetA BI,. (20) 


5. These results may be checked, to some extent, by taking the limiting 
ise of a sphere. In the first place we may calculate directly the case for 

the sphere by the same method. For a sphere of radius a, we obtain 
Gy Acrpc2a®k? | sec®@ e—2kof sect? dQ, (21) 
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For G, we expand the corresponding term ¢, in spherical harmonics, and 
we find that G, reduces to the same expression (21) with a negative sign; 
thus the total moment is zero, as should be the case. Turning to the 
expressions in (9) and (19) for the prolate spheroid, we find their limiting 
values for e > 0 reduce to the correct values for the sphere. 

6. Returning to the spheroid, we notice that G, may be positive or 
negative according to the speed; on the other hand, G, is essentially 
negative. Further, from (20) we have 

G, BRix,A. (22) 
If /, and k, are the inertia coefficients for axial motion and transverse 
motion respectively, we have 2e3A (1—e*)(1-+-k,) and a similar relation 
between B and k,. Hence (22) may be put in the convenient form 
G, (1+-k,)R/«o(1+4h)). (23) 
The ratio (1+-4,)/(1-+4,) is unity for a sphere, and approximates to two 
for a spheroid of large length-beam ratio. When c— «, or «,— 0, the 
integrals in the expressions for R, G,, and G, all reduce to the integral 
given in (21), which can be expressed in terms of Bessel functions: hence 
we may find the limiting values of these quantities as the speed increases 
indefinitely. It appears that as c > oo, R becomes zero of order c~*: on the 
other hand G, and G, approach finite limiting values, with 
G, > Sagpae’ A2/f?, G,—> —$mgpa*e®A B/f?. (24) 
Thus the moment G approaches the limiting value 
G G, + Gs > emg pa*b*(1 + ky)(ke ei if", 
and this is negative for a prolate spheroid. 

Some numerical values have been calculated from (20) for a spheroid of 
a length-beam ratio of 10. The moment at low speeds may be positive or 
negative and is small numerically; after a Froude number, c , (294), of 
more than about 0-4 the moment remains negative and increases rapidly 
towards its limiting value. 

[It may seem unexpected, as compared with surface ships, to find the 
moment remaining negative at high speeds. The model of a surface ship 
is usually allowed to trim and at high speeds it takes up a position with 
bow up and stern down, corresponding to a positive moment; the attitude 
of the model is then roughly parallel to the mean line of the water surface 
in its vicinity. But the submerged spheroid we have been considering has 
its axis maintained horizontal; so we may describe it roughly as being in 
a stream whose effective direction in the vicinity of the spheroid is inclined 
to the axis and this provides a moment tending to increase this angle, that 
is, anegative moment. For a numerical case take a spheroid with a 10b 
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und immersed to a depth f = 2-5b; we calculate the part G, of the 


, small angle 5 to the axis the moment tending to increase this angle is 
te pa 2(k, h 


we find that G, would be accounted for in this way by an angle 6 of about 


,)5. Comparing these two moments in this particular case, 


(03, which seems not unreasonable. However, this comparison cannot be 
pressed far; it is only intended to indicate a possible physical explanation 


{the negative moment at high speeds. 


7. Consider now any solid of revolution which, so far as axial motion in 
n infinite liquid is concerned, can be specified by a known axial source 
listribution. The part G, of the moment can be obtained at once by the 
method used in the prey ious sections; but it is not possible, in general, to 
ilculate the part G,. Turning to the connexion between G, and R for the 
spheroid given in (23), it is proposed now to use this as a suitable approxi- 
mation for any solid of revolution, and in particular for one of large ratio 
flength to beam. The inertia coefficient /, can be calculated; if M is the 
total moment of the given axial source distribution and V is the volume 
of the solid, we have 47M (1+k,)V. It is not possible, in general, to 
caleulate k,. However, for a long slender solid, k, is small; on the other 
hand, /, approximates to the value unity which it has for the transverse 
motion of a circular cylinder. Thus, in such a case, it is sufficient for a 
fairly close approximation to take 

Gs 2R/ ky, (26) 
where R is the wave resistance of the submerged solid. The simplest case 
sthat of the solid specified by a single source and sink. If m is the strength 
ifthe source or sink, 2A the distance apart, 2/ the axial length of the solid, 


nd 2b the maximum beam, we have 
tonlh c(l2@— h2)2: 4inh ch?(h?+-b?)!. (27) 


Taking the axis at depth f, the velocity potential can be written down 
to the same approximation as for the spheroid in (1). The process of 
letermining R and the part G, of the moment is the same as before, and 
the details need not be given. Using (27) to express m in terms of the 
limensions, we obtain 

R = 2ngpx,b*(1+-b2/h?) | {1—cos(2«yhsec 0)\e—2xof see*? sec39 dO, (28) 
j 


G, 2argpK, hb4(1 |. h2/ fh?) sin(2K,h sec @)e—**o/ sec*d sec dé. (29) 
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j 
For G, we should work out the next approximation for the velocity | 
potential, as in the case of the spheroid; but this does not seem feasible 
for the given solid. Meantime, as already indicated, we shall take (26) as 
giving a sufficient approximation and thus we assume 
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G, 4ngpb*(1+-b?/h?) | {1—cos(2x,yh sec 0) te-2x0f se*8 sec39 dA. (30) 
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(29) and (30) either by direct quadrature or by asymptotic expansions ou 
suitable for large values of the parameter 2«,h. To show the nature of the ha 
results calculations have been made for an ovoid with h 10b, giving a a | 


length-beam ratio of about 10-5. Two depths of immersion were taken, and } ed 


the results are shown in the figure with values of G/zgpb*(1-+-b?/h?) graphed in 
on a base of Froude number c/,/(2gl). Curve A is for f = 2-5b, and curve B fr 
for f = 5b. Curve A shows the typical oscillations at low speeds due to Wi 


interference between bow and stern waves: these would no doubt be 


damped by viscous effects in an actual liquid. For curve B at greater E 
depth these oscillations are too small to be shown on the scale of the or 
diagram. tc 
ti 
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ON THE LINEARIZED POTENTIAL THEORY OF 
UNSTEADY SUPERSONIC MOTION. II 


By K. STEWARTSON (The University, Bristol 8) 


[Received 2 May 1951] 


SUMMARY 

The velocity potential on a nearly plane wing executing simple harmonic motion 
in a uniform supersonic stream is obtained in a form suitable for systematic com- 
putation. Near subsonic trailing edges it is more convenient to use the acceleration 
potential, and if the wing has no supersonic leading edges the expressions become 
very complicated. 

Convenient forms for the solution when the flow is steady are also obtained, and 
if the wing has straight edges lying entirely within the Mach cone from the tip, a 
transformation is given correlating unsteady and steady flows. The assumptions of 
the linearized potential theory are used throughout. 


1. Introduction 

[IN a previous paper with the same title (1) the lift distributions on some 
straight-edged laminar wings, oscillating in a supersonic stream, were con- 
sidered on the assumption that the disturbances in the stream were suffi- 
ciently small to linearize the equations of motion. The problems were 
treated by means of the Laplace transformation and explicit solutions were 
given for the symmetrically swept back plane wing whose edges are outside 
the Mach cone from the tip, and for a semi-infinite plane wing, whose edges 
are either parallel or perpendicular to the direction of flow. It was pointed 
out that the method could also be used if the wing were semi-infinite but 
had edges inclined to the direction of the stream, provided that neither was 
a subsonic trailing edge, and that one of them was a supersonic leading 
edge. Of the straight-edged wings excluded by this condition, the most 
important is the ‘delta’ wing which lies entirely within the Mach cone 
from the tip, and a method for obtaining a series solution for the particular 
wing was given in the paper. 

Since that paper was written a new method has been found by 
Evvard (2) and developed by Ward (3) to determine the lift distribution 
ona nearly plane wing of arbitrary plan form, fixed at small inclination 
toa steady supersonic stream. By means of this method the lift distribu- 
tion on all types of nearly plane wings may eventually be found, although 
if the wing does not possess a supersonic leading edge (so that at no point 
are the flows on the upper and lower surfaces of the wing independent), 
the method leads to an infinite series. 

In addition, several methods have been proposed for tackling the 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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corresponding unsteady problem. Evvard (4) has generalized the method 
which proved so successful in the steady problem to cover unsteady motion, 
but his solution appears to be incorrect. Galin (5) has transformed one of 
the space coordinates simultaneously with the velocity potential, and in his 
new coordinates his new velocity potential satisfies the linearized potential 
equation for steady supersonic motion, so that the unsteady problem is 
reduced to a steady one. Unfortunately, the method is only applicable to 
a certain class of straight-edged wings, as may be seen from a more detailed 
discussion in section 5 at the end of this paper. The most important wing 
which may be solved by this transformation is the rectangular wing, but 
in this case the lift distribution may be determined more concisely by 
another method (1). 

Gardner (6) introduced a new coordinate into the differential equation 
satisfied by the velocity potential and divided the equation into two, which 
he proposed to solve successively by means of Evvard’s method. No explicit 
solutions were, however, given. Miles (7) used the Laplace transformation 
and the Wiener-Hopf technique to determine the lift distribution on some 
straight-edged plane wings. His method is analogous to the one used in (1) 
and is subject to the same limitations. Finally Stewart and Li (10) have 
developed Evvard’s theory (4), but their solution in the particular case of the 
rectangular wing does not agree with those obtained using other methods. 
These methods agree within themselves (7) and we shall see later that 
Stewart and Li have omitted a term from their solution. 

The advantage of Evvard’s method when the flow is steady is the great 
simplicity of the solution to the fundamental problem, which in turn 
depends on the solution of a certain integral equation, connecting the 
normal velocity on the wing, where it is known, with the normal velocity 
at other points in the same plane, where, in general, it is not known. Un- 
fortunately the corresponding integral equation in unsteady flow does not 
appear to possess such a simple solution, although Evvard (4) and 
Stewart (10) do claim to have found one. 

Instead of obtaining the solution of this equation it was decided to 
express the velocity potential at any point P, not only in terms of its given 
normal derivative, but also in terms of the values it, itself, takes in a 
certain region upstream of P. This region is to a certain extent arbitrary 
and so it was chosen to be as far as possible from P, subject only to the 
condition that the corresponding region where the normal derivative was 
required, was entirely on the wing. The method is suitable if a solution is 
desired as a power series, or if the velocity potential is required at all points 
of the wing. In the latter case, the solution may be started either from the 
steady solution, or from those given in the first paper (1). 
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The determination of the lift distribution on plane wings may therefore 
be carried out if the motion is unsteady, just as it could be if it were steady. 


The ‘delta’ wing, of course, stands apart from other plane wings and, in 


the steady case, leads to an infinite series. An alternative method of 
solution is suggested in section 5 of this paper, especially suitable when 
the lift distribution over the entire wing is required. Corresponding to any 
point P two domains of the wing are specified, and the velocity potential 
at P is expressed in terms of the distribution of the normal derivative of 
the velocity potential over that nearer P, together with the distribution 
of the velocity potential over the other, more distant from P. If the flow 
is steady then the expressions are simple, but generalizing to the unsteady 
problem makes them very complicated. Another method is therefore pro- 
posed when the wing has straight edges. This method is analogous to 
Galin’s (5), but instead of transforming 2, a space coordinate, we trans- 
form the hyperbolic distance r from the tip. This is much more convenient 
because in the plane of the wing we either know the velocity potential or 


its normal derivative all along the line of integration in the transformations. 


2. The equation of motion and the boundary conditions 

The notation used in (1) followed that adopted by Gunn (8) and is 
different from that commonly used at present in papers on the linearized 
potential theory of supersonic flow over thin wings. In this paper we shall 
use the following notation which will bring it more into line with present 
usage 

Let O(a, y,2) be a right-handed system of Cartesian coordinates, such 
that Ox is in the direction of the main stream, in which the undisturbed 
velocity is L’, the velocity of sound is c, and the Mach number is M. The 
presence of a thin oscillating wing will cause the velocity of the fluid to 
change slightly and we shall suppose that the velocity potential of the fluid 
is then Ui (a dh) 


The differential equation satisfied by ¢ is derived in (1) and is 


a2d, _ o2 __ o2 
02024 oy &¢ , yi e (2.1) 
ot cxot ox” 
and the excess pressure at any point is given by 
y Ch 7 
Ap pl ( 1 og : (2.2) 
ot Ox 


where p is the density of the undisturbed stream. 

We may assume, firstly, that the wing is executing simple harmonic 
motion about a stationary mean position in the plane z = 0, and any 
results we obtain under these circumstances may then be generalized by 
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means of Fourier’s integral theorem, to give the solutions of general 
problems involving gusts and quasi-periodicity. It will be assumed that 
the motion of the wing takes place in a direction normal to the plane z = 0 
and that the direction of the wing normal is always nearly the same as the 
direction of the z-axis. If the motion is simple harmonic with period 27/w, 
¢ will depend on ¢ only to the extent of an exponential factor, and we may 
write 


c(M2—1) 


subsequently y will be referred to as the velocity potential. The differential 


. «M 
f(x, y,2,t) = x(x, y, cJexp ia ; mi (2.3) 


equation satisfied by x is 


a2 2, (22 : 
X+—X = Bs X42), (2.4) 
oy? 62? ox? 
where M?—1 B2 and w= cAB?. (2.5) 
The excess pressure at any point is given by 
Ape iAz/M — | -pU? c (xe iAx/M) (2.6) 
Cx 


The changes of notation from the first paper are therefore that 
Y>2, z>4, x>y, VoU, and a>B. (2.7) 
It is hoped that the change will not prove confusing. 

The boundary conditions for such a wing as we shall be envisaging have 
been exhaustively discussed by Ward (3) in the corresponding steady 
problem. The treatment in the unsteady problem is the same except that 
there is also a time factor entering into the boundary values of ¢ and @¢/éz. 
Since the mean position of the wing in the plane does not change with time, 
we can write down the boundary condition from Ward’s paper, using Fig. 1 
and defining the characteristic coordinates « and B by 

x = x— By, B= «+ By. (2.8) 

On the wing Q, M, M, T, Q, Q, the normal velocity, and therefore ¢y/éz, 
is prescribed. 

Off the wing the pressure must be continuous, and hence, since y is an 
odd function of z, Ap = 0 in the rest of the plane z = 0. Therefore, 


x(v,y,0+) = A(—2By)exp a (a+ By) = A(a—p)e*P'™ (2.9) 


M 
in the rest of the plane z = 0. This new function Y is non-zero only in 
the region S, 7, Q, 7, S,, by continuity, and its value may be determined 
there by a generalization of the method used when A = 0 (5), if we assume 
the Kutta—Joukowski condition that at a trailing edge (for example 
T, Q.T, in Fig. 1) the velocity must remain finite. It will, however, be 
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more convenient to use the acceleration potential @ = —Ap/pU?, for this 
function is zero everywhere on the plane z = 0, except on the wing. 

Let us choose the origin O of coordinates in the plane z = 0 so that the 
lines a = 0, 8 = 0 just touch the periphery of the wing at M, and M,. In 
the region bounded by the lines a = 0, 8B = 0, and by the wing, there is 
no disturbance and merely for mathematical convenience we shall suppose 











By S; 
Fic. 1. 

that the wing occupies this region as well, (éx/éz),.) being equal to zero 
there. Since a disturbance at a point P cannot penetrate outside the 
downstream Mach cone at P, this will imply that » = 0 there too, and 
so there is no inconsistency. The integration of any relevant function over 
this region will be zero. We shall also assume to begin with that the wing 
proper is not pointed at the leading edge, but this may usually be allowed 
lor by using the steady solution near the tip, and then extending the 
solution downstream using the method we shall develop here. 


3. The effect of a subsonic leading edge 
In terms of the characteristic coordinates, « and f, the partial differential 
equation satisfied by x is 
c* : ory ‘ 
1B2—X 1 B22, — —X, (3.1) 


Q ~ 
CXxCp Cx 
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and it is easy to verify that a particular solution of this equation is 
cos{A(aB — B2z?)!} 3.9) 
- — (3.2 
(aB— B2z?)! 
Since @y/éz is prescribed and is continuous across the wing and every- 
where else on the plane z = 0, it follows that y is an antisymmetrical 
function of z. If, for z > 0, we write 


* [ cost : B’ 252]! 
x(a, B, 2) {| - {A| (a—a’)(B—P’)— B 


ue oe ee 
‘ la’dp’, 3.3 
((a—a")(B—B’)— B27} f(a’, B’) da'd oy 


the integration being taken over the area bounded by 


7 oe 





roa. £ B’, and (a—a’)(B—f’) > Bz, 
then if z < 0 we merely have to reverse the sign of y. 
We shall write 
X1(a, B) x(a, 8, 0 |) and x(a, 8) ‘3 x(a. B,2)]] ; (3.4) 
Cz 2-0 


The complete solution is known if we can determine f(a, 8) which is simply 
related to y., for 


x: 2 Bx f (x, B). (3.5) 
Once we know yx, in the region a’ x, B° < B the value of y,(«, 8) may 


be computed and also the pressure fall from the underside to the upper- 
side of the wing, which is 
ro) C X ir ! 
») Al (* > 
2pU?{- oh > (3.6) 
, (x M x} 


The solution is therefore known in the region OM, CM, O of the wing 
because xy, is known on the relevant portion of the plane, but the 
potential at points on the rest of the wing depends on a region in which 
x. is not known. To calculate the potential at P we should first have to 
determine y, in the curvilinear triangle P, M, P;. In the steady problem 
the contribution to y,(P) from this curvilinear triangle may be accounted 
for in a very simple way, because it is equal and opposite to the contribu- 
tion to x,(P) from the curvilinear triangle P, P| OM, P,. Hence the con 
tribution to y,(P) from the parallelogram P, P| OP; is zero and we only 
need the contribution from the region EF, of Fig. 1 to determine y,(/P). 
When the motion is unsteady Stewart and Li (10) assert that the same is 
true, but we shall show that this is by no means the case. In this paper 
we shall express the contribution of y, in the parallelogram P, P| OP, in 
terms of the distribution of y, in the region F, of Fig. 1. 

Let 8 = A(a) or «a = B(f) be the equation of the wing boundary from 
M, to Q, and let us evaluate y, at the point P(a«,f) in Fig. 1, where 
vy > a > 0 and B, > B > By. This region excludes all points affected 
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either by the subsonic leading edge M,T, Q,, or by the subsonic trailing 
edge 7; Q, Tj, or those points at which x, is already known. 
Denote by a bar the Laplace transform of any quantity with respect to 8 


so that, for example, 


X(a, Pp, z) | ¢ PB (x, B,z) ap. (3.7) 
Then x satisfies 1B2p“X4 Bye — SX (3.8) 
CX Cz* 
since if B 0, x Q. 
We shall express the value of y,(a~) where a > a* 0 in terms of the 


ilues of x, in the interval (a*,«) and of x, in the interval (0,a*). The 

ost convenient value of «* will be taken. 

The differential equation for ¥ is an equation of the type commonly 
wcurring in problems of heat conduction, and the solution given here of 
the following problem may be found in Carslaw and Jaeger (9). 

[f ov/ot x 0*v/ax", where « is a constant, (3.9) 


j 


0 when ¢ 0, and v A(t) when x 0, then, if t > 0, 2 > 0, 


x § by exp} x?/4n(t—t’)' dt’ 


- (3.10) 
2(7K) (¢t—t ) 


In particular, if x Q, 
cv | F | 


Cx (7K)* | (t t’)? 


0 


“ {h(t’)\ dt’. (3.11) 
ot 


We can deduce the solution of the problem when the boundary con- 


litions are changed to 


v 0 when f 0, and év/ex %(t) when x 0. 
In this case if a Q, 
K : it’ 
Y ue). (3.12) 
(7K)? | (t—t )? 


These results may be applied to our problem if we replace v by xe™*?, 
by z, t by a, and « by 1/4B*p. Then, if z = 0 and 0 < a’ < a*, 


\ r \i da” c ” 2 ‘ ‘ 
rere fa} 2B{-)° | (x’— a" )b 5 {Xi (a” eo" 47}, (3.13) 
77 | x 
U0 


x )? Ox 


nd if a * 


e" VAD Y (cx) - X,(a’ lero Ap , (3.14) 
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which may be written as 


] - da’ ‘ 
¥,(a) = — | —~ ¥, "\o—?/Ap(a—a’) 1 
X1\% 2 B( pz)! . (x ” yi Xe ) 
] . da’ “ dx” C a re ee ae 
+: | (a—a’)! | Zomae eens eae 
Wy JA 
0 0 


The second integral may be simplified if we interchange the orders of 
integration; it becomes 


r da’ [a—a*\t__ , nya 
( * =) Ral Je-(a-a' MAD (3.16) 
xX x — 


x 





7 
« 


0 


Hence in terms of £, if J, is a Bessel function of the first kind 


B : 
rs f} + pr Cos[A{(x—a’ )(B—B’)}4] 
(a, 8B) = —— ip da’ y.(a’, 4 
xl B) = — ap, | AB | a! xP Pp 
0 
l i da’ X y* 1 A : - , 
| i % J? xl) 5 | dp’ | dx’ x;(a', B’) > 
7 xX x ~~ x “7 . 


0 0 Bip’) 
ms (= ‘ “| J,[A{(a—a')(B—B’)}* | (3.17) 
eo {(x—a’)(B—B’)}# 
since y,(a’,B’) = Oif a’ < B(p’). 

The best value of «* is the least possible consistent with the necessary 
requirement that y,(«’,8’) should be known everywhere in the region 
v* < a’ <a, 0 < Bp’ < B, so that on one hand the influence of the last 
two terms of the right-hand side of equation (3.17) should be as small as 
possible, while on the other, it should be possible to calculate the first term 
directly. If a* < B(8) the region «* < a’ < a includes points at which yx, 
is not given and so the best value of «* is B(8). This has an additional 
advantage because then the second integral is zero in virtue of the 
boundary condition that if « B(B) then y,(a, 8) = 0. With this value of «* 


cos{A| ( r—«’)(B —B’) |} 


1 ore 
(a, B — - dx'dB’ x.(a’, B’ 
_— 2Bz | | OP Pa BBY 
ai ew » an{ x—B(B)\$ A [A{(a—«’)(B—B’)} | ‘ 
. da'dB’ y,(x’,B aie : 3.18) 
ar | | oe RNS Fa 2 {(a—a’)(B—B’)}! 


As A> 0 the last integral tends to zero and we recover Evvard’s solution 
(2) for the corresponding steady solution. Stewart and Li’s solution (10) 
of the unsteady problem omits the second term on the right-hand side of 


equation (3.18) but this term may be verified to be non-zero in general. 
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The solution may be regarded as an integral equation for x, of a kind 
suitable either for systematic computation or successive approximation. 
The method to be used when y, is required everywhere is as follows. 
Suppose x, is known in the region F,. Then by means of equation (3.18) 
it may be calculated along the line PP; as far as the line M, Q, and, of 
course, at points upstream of these. Since every point of the line PP, is 
outside the region £, except for P, where x, = 0, no complications arise 
due to singularities in the integrand. Having determined x, along P, P, 
we may now move P, a little way towards 7; and since y, vanishes on the 
edge of the wing, it is known everywhere in the region £, defined by the 
new position of P,. The process may now be repeated until we reach 7, Q,, 
after which we must use the method given in the next section. For suc- 


cessive approximation let the nth approximation to y, be x{, where 


| # » pr, COSA{(a—ax’)(B—B’)}# , , 3, 
0 0, x.(«’, B : la’d 
XI m 2Br | | XAG SP ((a—a’)(B—B’)}3 ali 
! 
ind 
A (ft 9) 
x1 Xi" , fe | (x? = 
"Es 


/ > 1 f ‘ B’)\4 
\ ed eco v' \(B—B')} dx'dp’. (3.19) 


B(B)—« {(a—a’)(B—B’)}! 


Since the lowest power of A occurring in the expression for y{”— x{"~ is 
2n—2, x will tend to y, rapidly near M, but more slowly further down- 
stream. Equation (3.18) is not convenient when the potential at one 


particular point is required. 


4. The acceleration potential 
The solution given in the previous section is correct provided that 


v, and B < f,. If a > a», then we must take account of the second 
subsonic leading edge M, 7, and its effect may be computed from equation 
3.18) after making slight changes, the chief of which is to interchange « 
ind 8 on the right-hand side. In fact, we proceed in an analogous fashion 
to that when A = 0 (5). 

If, however, 8 > B, the solution will be different in character owing to 
the possibility of there being a vortex sheet (a discontinuity of x across 
the plane z = 0) downstream of the trailing edge 7, 7, T,. The excess 


pressure must be zero of necessity there, and so in this region 


x1 = CMF (xB), 


where .¥, is at present an arbitrary function, from section 3. When the 
flow is steady the Kutta—Joukowski hypothesis that the fluid velocity 


5092.18 


L 
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must be finite at a trailing edge is used to determine A, and hence also 
the excess pressure on the wing. 





This method may also be used when the flow is unsteady if we assume, 
as we shall, that the Kutta—Joukowski hypothesis still holds, but the 
complication is very much greater. It is much more convenient to use the 
acceleration potential 6, defined by 

Ap ex WU 


a 4. 
pU? ox M* - 


which has the advantage that it gives immediately the excess pressure on 
the wing, whereas if we determine y, then @ has to be obtained by differen- 
tiation before we can discuss the aerodynamic forces. The acceleration 
potential satisfies the same differential equation as x, and the boundary 
conditions on @ are that on the plane z = 0, 6 = 0,sgnz, where 6, is zero 
except on the wing, and where @6/éz = 6, is prescribed. 

These conditions are formally simpler than those on y because 0, is 
identically zero off the wing, but we have to be rather careful about possible 
discontinuities in y, which would imply that 6, would be undefined, and 
possible singularities in @, at the edge of the wing. At supersonic leading and 
trailing edges, and also at subsonic trailing edges if we make use of the 
Kutta—Joukowski hypothesis, 6, is bounded, but at subsonic leading edges 
such as MV, 7;, it has an algebraic singularity with exponent —}. Both of 
these difficulties may be overcome provided only that we replace the 
discontinuity in x, by a region of rapid change. 

The appropriate solution for 6, is, from equation (3.17), 





pena 
0,(,8) = — wel da‘ag’ OOM(o—o' WBBM g 
. {( x X )(B -B )}? ” 
0 
if A heme 
+ nad os Jéx(a’.)— 
7 P "= ero 
B(p) 
rfp x=" )t Hoof) 
: | dp’ | it ———F° 6, (4.2) 
2r J ! x* — oy (a—«’)(B—f’)} 
0 B(B’) 
if « < a,and P lies on the wing. The region of most interest here is 8 > £,, 


so that P, is on the subsonic trailing edge, because 8 < f, has already been 
considered. In this case 0, is bounded at P, and so as a*-> B(f)-+, the 
second integral will tend to zero and the third will be integrable in the 
limit. The only complication will therefore arise from the discontinuity 


in xy, along the line P, P,. If we replace this discontinuity by a rapid change 
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in y, and carry out the integration with respect to x, across the line P, P, 





then 
B X a 
- ] : | dp’ | da’ 9 cos At (a | x \B—B } 
2Br J “3 F {(a—a’)(B—B )\ 
0 B(p) 


2 OS NM , — , 14 
a l da'dg’ COM -2'\(B - ; lp 
2Br J J {((x—a’)(B—p’)}#  * 
Ey 
l ; is 1, COS| A{(a— x’ )(B—A(a’))}4 
. xia’, A(a’)} ; Avi 
2Br. | (a—a’)(B—A(a’))}! 
B(B) 
since the equation of the line P, P, is «’ = B(p’) or B’ = A(a’). 
When A = 0 equation (4.3) is equal to —Apsgnz/pU? and hence we have 


haa’) I} da’, (4.3) 


the solution in a simple form, comparing favourably with Ward’s (3) for 
the same problem. In general, if «a < a», the method of systematic com- 
putation given in the previous section may now be extended right down 
to the trailing edge. 


As an example of the use of the method which we have developed let 





us consider the semi-infinite wing « > 0, y < 0 on which the normal 
velocity is X, = e“/M, so that 0, a. (4.4) 
The solution of this problem has already been given in substance in (1), 
and it may be deduced that if 0 < z - By, 
l l 72? 
6 Jo(Ax) ] — + O(A‘)), 4.5 
on a ;.o (4.5) 
while, if x By > 0, 
2 7 
0, a Jy{A(x?— By? secty)*|} dy, where sec*y, x/ By, (4.6) 
0 
2 A2x*\ . By\} 
] }sin 2) 
7B 1 x 
A“ > 5 4 4 | 
Ths 2 By)(a+ By)!(— By)*+-O(A4)]. (4.7) 
Since 22 (= a+) is zero on the line P, P,, we have for the first approxi- 
tion to 6, in the region 0 By < x, from equation (4.2), 
a1) 2 [ da’ cos|A{(a—«’)(B+-«’)}*] (4.8) 
, 27B \(a—a’)(B+-«’)}4 , 


= | Yaersin ( 2} Se |-2 By)(a+- By)*(— By)! 00). 


(4.9) 
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This first approximation which is correct only as far as terms of O(A) 





may now be substituted back into equation (4.2) to obtain the second 


approximation. We have 


ap 2 1 { , __ A’ rs 
a) — ara I ( <n BEM Ta da’ap’ (4.10 


) 
«TI 


y ((a—a’)(B—B’)}* 
since « = Bf on M,T,. Hence, 
x 


¥ Ffo-p,1.,2f x .Je’—P 
(2 QD =d } sin! dp’ | +- O(A4 
1 : a! (; °) si = (=F) ) 
0 


(4.11) 


i — ~e +- By)?(— By)!+ O(a‘). (4.12) 


Substituting for #{? from equation (4.9), we obtain 6 which is correct 
as far as terms of O(A*) agreeing with the known solution, equation (4.7), 
as far as this point. The procedure may be continued as far as required, | 
and each successive approximation will be correct to a higher power of )’. 

Although the solution may be computed from equation (3.18) when 
x vg, it is not in a convenient form when 8 > f, as well, for then we 
need y, to compute the effect of the subsonic leading edge M, T, and 6, 
to compute the effect of the trailing edge 7, T}. We shall overcome this 
by obtaining a formula for 6, when P, is on the subsonic leading edge M, T,. 
This is given by equation (4.2) except that the first integral has to be 
modified in the way already discussed, and the second integral is no longer 


2 


zero in the limit as «*—> B(8)+. Near P, 6, is singular and may be 





written as C(B) 
{1 — B(B)}’ 


where C is to be determined, and the contribution from the second integral, 


(4.13) 


~ 


1 (£ da’ r— a* \3 ; C(B) 
1" 6. (x' Bb) > as a* > B(B)+. (4.14) | 

FE | oe ae 8 

BiB) 

We have therefore to add this term to the solution appropriate to the 
subsonic trailing edge. Returning to equation (3.18) we may easily show 
that near the subsonic leading edge y, is given by 
X11 B(B)}! | 

1 [( dB’ x,(a,B’) A x1(a’, B’) da’dp’ J,[Af( B(B)— a’ )(B—B’)}"] 
Br) (B—B’) 27 )) {BiB)—a’}! — ((B(B)—a’ )(B—B'} 
Es 
0 7 | : (4.15) } 
Since 6, OX, OX _ Xp 
Ox op M | 
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O(A) | then the right-hand side of this equation is equal to 
cond 20(B)/(1— B’(B)). (4.16) 


The method of systematic computation is now clear. Starting off at the 

4.10 most forward point of the wing we systematically calculate either x, or 6, 
it all points as far as the line 7; 7, using the method given in the previous 

section. Either of these functions has its disadvantages, for if we calculate 

y,, then 6, must be found afterwards to determine the excess pressure on 

t) the wing, while if we compute 6, then we must also compute x, simul- 


taneously because we need it to determine the singularity in 6, at the 


(4.11 subsonic leading edge from equation (4.15). However, since it is easier to 
(4.12 integrate than to differentiate numerically, it is probably better to calcu- 


ite 6, as the basic function. Downstream of 7', 7, we calculate 0, only, 
1 ite 1 3 


orrect making use of the modifications to equation (4.2) already discussed. Since 
(4.7), | @. is finite at a subsonic trailing edge the function C() will be zero except 
uired, ipstream of 7; 7, where it is already known. 
* of A. When the flow is steady A = 0 and the results we have obtained simplify 
when | considerably. We have, if « ve, B < B,, 
en we , 
ind 6 |p —_ _ l 5 dat (A () l)y.(a’, A(a’)) | [ [ da'dp’ 6(«’, B’) 
e this | 2Ba ((a—a')(B—A(a’))j# 2B J J {(a—a’)(B—B')}* 
M1, | on 

a LBB) xl BEB’) agro. 17 
longer 27(a—B(B))' J (B—Bf’)* 
lay be ’ : ae 

ind if « vo, B > f, the last term on the right-hand side is to be omitted. 
(4.13 These expressions may be compared favourably with those obtained by 


Ward (3). They are more convenient than his for computation because 
tegral, | they do not introduce any intermediary functions and because they give 6, 
_ directly in terms of @, and y, which are given, whereas in (5) 6, has to be 

(4.14) | calculated by numerical differentiation. 


to the | » Wings without a supersonic leading edge 

v show These wings are of quite general shape except that at their most forward 
| } point O, the wing edges form a triangle which lies entirely within the Mach 
one from O. The remaining part of the edges may be either subsonic or 


supersonic, leading or trailing edges, but we shall only consider subsonic 


>/\) | . . . . 
p ); eading edges; the solution near other types of edge may be deduced from 
3") 5 this one, taking account of the previous work in this paper. 


(4.15) } In accordance with the practice in this paper we shall build up our solu- 
tion step by step, so that to calculate the potential at P(«,8), it will be 
issumed known at all points in the upstream Mach cone through P. The 
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calculation is started near the tip by using the known results given, for 
example, in (1) where a series expansion in powers of A was developed. 
The wing is shown in Fig. 2 and the equations of its edges are 
x B,(B) or B = A,(a«) (5.1) 
and x= BB) or B= A,(a). 


Let us begin by considering the steady problem in which A = 0 and let 
the velocity potential be ys, on the upper side of the wing. The contribution 





to w%,(a,8) from the values of (é/,/éz),.) = 4, may be divided into three 
parts. Firstly there is a contribution from the parallelogram PP, QP, 


which is . B 
] i gilake ; us, ( x’, B’) il 
: - 1 | ek a Ae 5.2) 
27B | - | ” {(xa—a’)(B—B’)}4 " 
B,(B) Ao(x’) 


Secondly there are the contributions from the parallelograms P, M, 0Q, 
and P, M, OQ, and these, by Evvard’s theory, are zero. The third contri- 
bution arises because the parallelogram OQ, QQ, has been surveyed twice 
and so we have to subtract the contribution from it. The third contribu- 


tion is therefore 


A(x) B,(B) , B 
l a : + : h(x’, , il 
4 1B | 1 YAS >P) 5.3) 
mB) 8 | ye PP ” 
0 0 


Since ys, is not known everywhere on the parallelogram QQ, OQ,, the 
solution is not complete, but it may be completed by expressing #, in this 
parallelogram in terms of the values taken by ys, on the same parallelogram. 
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This function is known everywhere on the parallelogram since it vanishes 
off the wing and we have 





*8 x 
B a . d , eo , 
us.(a, B) =— | dp’ | — tila +B) rt (5.4) 
“TT Y J Wa—-a )(B B )5 
0 0 


where the star (*) denotes that the finite part of the improper integral is 
taken. Substituting this expression, in equation (5.3), inverting the order 
of integration, and adding to equation (5.2) we get eventually 
; } da'dp’ yb x’, B’) 

f Q 


27B . J i (ax x )(B p’)} 





; ; - : 0) 
(a— or \(B B ) (B,(B) — OL )(Ag(a) —B ) 





' [ b,(«’, B’) da’dB’ | (a— B,(B))(B—A,(«)) | (5.5 
S: 

where the regions S, and S, are shown shaded in Fig. 2. Thus the potential 
it any point P of the wing has been expressed in terms of the normal 
derivative of the potential over the region S, together with the potential 
wer the region S,, and is in a form very suitable for systematic compu- 
tation. 

The solution may also be obtained in the corresponding unsteady 
problem, but we shall not give it here, since in the first place the argument 
is similar, and secondly the result is complicated, involving ten integration 
signs. The reason for this is that the contributions from the parallelograms 
UM, P,Q.O0 and OM, P,Q, do not vanish and have to be calculated in the 
same way as in section 3, and the contribution from the parallelogram 
0Q.QQ, does not simplify. 

Another method, however, may be used in the most important wing of 
this type, in which the edges are straight and symmetrical with respect 
to the direction of flow. 

It was shown by Galin (5) that %(x, y,z), defined by the equation 

2\14 

ub(a, y, 2) x(x, y,Zz)+A ve. ¥, 2) a) dé, (5.6) 

satisfies the differential equation of steady supersonic flow, and further, 
that the transformation from y to % could be inverted merely by inter- 
hanging y and % and replacing A by 7A in equation (5.6). If we are given 
yéz when z = 0 and y > 0, and y when z = 0 and y < 0, we can deter- 

ine the corresponding values of és/éz and %. Hence, from the theory of 

steady supersonic flow we may determine ¢% when z = 0+ andy > 0, and 
finally deduce y when z = 0+ and y > 0 by inverting equation (5.6). The 
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method depends for its success on knowing x for all x when z = 0, and T 
y < 0, and éx/éz for all x when z = 0 and y > 0, so that the method may _ | 

be used to determine the lift distribution on the semi-infinite wing lying 

in that part of the plane z = 0, defined by y > 0, 2 > Béy where @ < 1. 

If we rotate the axes about the axis of z the lift distribution on wings lying | 





in that part of the plane z 0 defined by B6,y < x < BOé,y may also be Le 
2! ib si 
determined in a similar way, provided that #6, < 0,and @, > 1. In general, 
however, it is not possible to divide the plane z = 0 into two regions, on , 
one of which % could be specified from the boundary condition on y by Th 
means of equation (5.6) and on the other éy/ez could be specified from the 
vr P - ° e ° mn ' 
boundary condition on éy/éz in a similar way. The steady problem cannot 
be specified, therefore, in a way that is suitable for solution, and so Galin’s 
method is not applicable to the nearly plane wing of arbitrary plan form, 
and in particular it is not applicable to the ‘delta’ wing. We may, how- 
ever, proceed as follows, using Galin’s idea of transforming from the 
unsteady to the steady equation. 
Let us introduce new coordinates r, u, v defined by 
x r cosh u, By rsinh wu cos v, Bz r sinh w sin v wl 
and let the wing lie in the plane » = 0 and be bounded by the lines 
u Lu,, where wu, is a constant. The boundary conditions are that 
x= 0 when |u| > |u,|,v = 0 (5.7) - 
‘ | sic olive "} . ‘ ‘ rf) ki 
and @y/év is given when |u Ur|, 7 
Further, y satisfies 
5OX . ee. eee l ae ie C l o*y ” 
~ x 1 a °X ry (sinh wo +-.—.. x (5.8) 
er or sinh w cu cu} © sinh? ev? 
Since y vanishes when r = 0), it may be verified by direct substitution, o1 T 
by performing the operations successively, that 
X(p, u, v) | x(7, uv, v)r?(p?+-A*)le-Pr dr (5.9) 
0 a 
— ee cx . . : 
satisfies (A*+- p?)_“+ p 1x + L(x), (5.10) | 
ep* op } 
f 


where L(y) has been written for the right-hand side of equation (5.8). 


Moreover, if p s—A* 4s, ¥ also satisfies 


(5.11) 
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and This equation is the same as the previous one with A 0 and hence, if 
nay a+io 
; l ‘ ds 
ing r sR ~ ; a nD 
5 b( R, u,v) = est 5 X(P, U2 ) (a 0), (5.12) 
] aml i: & 1. 
ing 
2, m 
F . 9 COTW Cw = ‘ 
» be | satisfies k2 = 2R ; Li). (5.13) 
ne ck? oR 
ral, 
Deas Thus & is the solution of a correlated problem in steady supersonic flow. 
by The relation between #% and y is 
the 
not 1 ¢£ dsekt , wd 
ca di( R : dae s\rhe—PT( n2_L 2)5 dr 5.14 
ub( Uv (7, U, v)ree ) ar. ). 
lin Ss ' Dai e R} X ) (} ) ( ) 
rm, ° ° 
OW The order of integration may be inverted and the inner integral is then 
the a i 
l ‘ 2 | dr er 
ds(1 ,jexp (R—r)s-4 : (5.15) 
2mt . 4s- 4s 
a 
which vanishes if 7 R, and if r < R, is equal to 
lines \R ] \S 2 > ~\V4 
Ah / r(R—r ° ~ > 
6(R—r)- lA M | (5.16) 


2 {r(R—r)} 


\?-() where 6(R—r) is the delta function, and J, a Bessel function of the first 


kind with imaginary argument. 


Hence, 
R 
AR* f x(r, u,v) 
1( PR , (PR , fy( RB “\h . e497 
(5.8 ib( R, u,v) x(k, u, = (Rr)! L|A{r( R—r)} | dr. (5.17) 
a t r)? 
0 
n, O1 Thus in our correlated steady problem 
ous ° 
is given when |u|! - Uy, 2 0 
Oz 
9.9 
and y% vanishes when u| > u,), v 0. 
= 10) This problem may be solved by means of the method given in the first 
>, ‘ . 
part of this section. To complete the solution we must invert the trans- 
formation from y to y%. This is done by reversing the steps of the argument 
above and we find 
a t= F 
ey R\} . ee . 
5.11) x(7, u,v) — d( R. u. o| “dR | = 2 8 dp: (5.18) 
. mt. / A (pr T x, 
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and since 2s = p+(p?+A?)* this reduces to 


. ) Ri 
x(7, u,v) = wr, u, "3 | a a 


0 


J[Mr(r—R)}*]. (5.19) 


The steady solution determines % when |u| < |u,| and v = 0+ for all R, 
and then y may be found from equation (5.19). 
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SUMMARY 

The laminar motion of a plane symmetrical jet of compressible fluid has been 
investigated previously by several authors, but no closed analytical solution has been 
obtained. This paper shows that such a solution is possible if the Prandtl number 
f the fluid is assumed to be unity, and the dynamic viscosity pw in the motion 
lepends on the absolute temperature ¢ according to the relation pz « t. In terms of 
the stream function % and the axial distance 2, a very simple function is obtained 
for the velocity parallel to the axis; ¢ then being obtained from Crocco’s relation. 
Graphs of velocity and temperature distributions across jet sections are given for 


1 wide range of jet Mach numbers and temperatures. 


1. Introduction 

THE motion in a laminar jet of fluid issuing from a long narrow orifice has 
been investigated by several writers. The solution of the problem when 
the fluid is incompressible has been fully given by Schlichting (1) and 
Bickley (2), and in the case of compressible fluids Pai (3) has dealt with 
a supersonic jet exhausting into a parallel supersonic stream, and 
Krzywoblocki (4) with the more general case. No complete closed 
analytical solution (under the usual boundary layer conditions) for the 
compressible jet has, however, been given. This paper shows that such 
a solution is possible under several simplifying assumptions, the use of 
which does not, however, limit the application of the solution excessively. 
The main difficulty, in fact, would be the attainment of a flow with 
sufficiently low Reynolds number and high Mach number to realize the 
conditions of laminar compressibility in the jet. It is possible that some 
application of the results could be made to low density motions and as a 


basis for stability calculations. 


2. Notation 
Ox, Oy axes parallel and perpendicular to the axis of the jet; O being 
some fixed point on the jet axis, 
u velocity of fluid parallel to Oz, 
v velocity of fluid parallel to Oy, 
p density of fluid, 


~ 


absolute temperature, 
uw coefficient of viscosity, 


(Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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v kinematic viscosity, 
C,, specific heat of the fluid at constant pressure, 
coefficient of heat conduction, 
L a characteristic length, 
Pr Prandtl number (= C, p/h), 
Sutherland’s constant (for air, approximately 117° C.), 
d length of nozzle in direction Oy. 
M, Mach number of jet at nozzle, 
y ratio of specific heats of gas in jet, 
J mechanical equivalent of heat. 
Suffixes: co denote values in the ambient undisturbed stream, 
j denote values at jet orifice, 


a denote values on jet axis, 


3. Basic equations and boundary conditions 

On the assumptions that the boundary-layer approximations can be 
applied, and that no pressure gradient exists throughout the jet, the 
relevant equations are 

equation of motion in the x-direction 


ou ou c ou 
pu—-t pv- (1 ) (1) 
a cy oy\ ey 
equation of energy 
pu MC!) ov Ar!) oe Sh BS (2) 
a oy oy\ ey) J\ey 
equation of continuity 
- (pu) -+ te (pv) 0, (3) 
0x Cy 


equation of state for a perfect gas 


pt Din be p;t;. (4) 


These four equations involve seven dependent variables, u, v, p, t, uw, Cy, 


and k, so that three more conditions or equations are needed to render a 
unique solution possible. These can be supplied by the assumptions 


C,, = constant, (5) 
Cp 
Pr pt — 1 6) 
k 





where C, m are constants. 
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Assumptions (5) and (6) are often made in boundary-layer problems, the 
difference obtained in the velocity distributions by taking Pr = 1 and 
Pr = 0-73 being small, as shown by numerical results obtained by Emmons 
and Brainerd (5). 

In equation (7), C is usually assumed to be unity, but the introduction 
of the constant C can be extremely useful, as has been shown by Chapman 
(6), for it is then possible to ensure that the approximate relation (7) gives 
in exact value for u at any required position in the field of flow. For 
example, if Sutherland’s equation is taken to be the exact form, i.e. 

: (( 5 (8) 
Bj tj} tH Ss 


then u can be made exact at any temperature ¢ by means of the relation 


(i) = (ies © 


It is most important in the general case to have » exact where the 
shearing stress gradient is largest, i.e. in the present problem, on the jet 
axis. This, unfortunately, is not possible at all points as ¢ varies along the 
axis, and for the calculations made for the graphs given below, » has been 
made exact at the jet orifice, i.e. 

C = 1. (10) 
The boundary conditions to be satisfied by the above series of equations 
are 


on jet axis (y = 0), 


v=0 and - 0; (11) 


at jet edges (y — +0), ne 


The first two of these condition follow from symmetry considerations and 
the third from the assumption of a static ambient fluid. 


4. Reduction of the basic equations 
As the Prandtl number has been assumed to be unity, equations (1) and 


) 


2) are satisfied by Crocco’s relation 


| og 
t— A+ Bu ay where A and B are constants. 
With the values 
t f when u U;, 
t ‘s when u 0, 


this relation gives 
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The introduction of the stream function ys, defined by the equations Wit 

pj Op , ; p5 Op “ 

Y= =; v= —-?_, (13) 
p dy p dx 

ensures that equation (3) shall be satisfied, leaving only equation (1) to For 
be solved. hav 
A change of the independent variables (x, y) in equation (1) to (x, #) is , 
now effected by the transformation formulae . = 
Fu 


(ea). = Sr oph is 


This operation reduces (1) to the form 


= = (uw 2 =) (15) + 


(=) eam ms +(5 : (14) mo 
ox), pj \Ob), \ex}y 


Cx Ob p; © As 
The introduction of the new variables 
, u t - x ; I 
l ae z : X 7 : 1 = carat 
U; 4 U,V; L 
\ ' 
J J ia (16) wh 
7 
A=f=0T™, a= | 
Bj Pj 
then gives the dimensionless form 
oU c , ,oU z so 
mee os a PM as OT — 1. (17) 
oX a ov | 
Equation (17) may now be fully solved if two further assumptions are 
made. These are that the linear dependence of yu on ¢ holds (i.e. m= 1), In 
and that the velocity profiles in the jet exhibit similarity at all points 
downstream of the orifice. 
This latter assumption is reasonable provided that the boundary-layer 
and mixing-region in the jet both vanish at the orifice, i.e. that the jet 5 
derives from a point-source. U will then, at any cross-section of the jet, 
depend entirely on a single dimensionless variable { (depending on X, ), in 
such that { vanishes at the orifice. It will be observed that the last 
condition is fulfilled in this case. 
It is then possible to set 
U = XP F(t), (18) 8 
where { = FUR". . (19) 
p, q, and r being numerical constants. 1. 
a ay a 
a 0 rear o¢ — ‘ 
Then oe qX"¥e-1 and — ryaXr-1, (20) 


ot ox 7 





15) 


16) 


7) 


ire 


3) 


9) 
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With these substitutions, equation (17) takes the form 


. - -aF ewe 9. . o dF m 
XP nF 1 rf | _— q?.X *P+eryt ‘a i Pp q-1 q (21) 
\ dt} of\ dt 
For equation (12) to be a differential equation involving ¢ only, we must 
have -] 
(22) 
and p+2r —1. 
Further, as there is no pressure gradient in the jet, the rate of flow of 
momentum across transverse sections of the jet is constant, i.e. 
2 | pu* dy = constant, M, say. (23) 
0 ; 
As pu = p,(0p/ey), equation (23) gives 
y . ; . . 
M 2 | p;udyp 2p;u; | Ud 2pul(v; LC)! | X?-F(C) dl, 


y=0 4=0 y=0 


7] x uv x 


whence 9 =f. (24) 
Solution of equations (22) and (24) gives 


p=r 4; q= 1, (25) 


so that (21) is finally reduced to an ordinary differential equation in @, 
d/{,dF\ lf dF 
F +-{F+{——\—0. 26 
z\ ai si a (26) 
Integration with respect to ¢ then yields 


,aF 


F — 
dt 


1¢F = a constant, D. (27) 
5. Application of boundary conditions 
In anticipation of section 7 below, it is seen that the transverse velocity v 
in the jet satisfies the relation 
ae OF v,C } 
na j-i~|? — —Qu, (28) 
3 L 
so that the first boundary condition of (11) is satisfied if and only if 


¥ 


C¢=0 wheny = 0, (29) 
i.e. if the jet-axis is the streamline 4 = 0. This condition is quite clearly 
acceptable, as the axis is, by symmetry, a streamline of the flow, and the 
numbering of the streamlines is, at the moment, purely arbitrary. 
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Further, 
ou c oU oF 
U; a, Qu - Quu, X-* OX -iuu, F’ u,v; LC), 
oy ” ey 7 Gp ’ ous i P'(O)/(ujv; LC) 
(30) 
Hence the second boundary condition of (11) is satisfied if 
F’(0) = 0. (31) 


6. Solution of the basic equations 

Application of the conditions (29) and (31) establishes the constant D 
of equation (27) to be zero, so that as F(¢) oc U does not vanish anywhere 
inside the jet, the equation (27) yields, after further integration, 


F(¢) = N—}2, (32) 
where NV is a constant of integration. 
If the position of the orifice is taken at 2 = L, then at the orifice, 
X 1, ¢=0, and « U;, 
so that NV 1, and F (Cf) — =. (33) 


The characteristic length LZ has then to be determined from equation 
(23), which gives 


46 


M P; u? d 2p; uj (v; LC): | (1— 3G?) dZ, 
0 
3M? 
i.e. L = ——._., (34) 
32p; UF 1; ¢ 

7. Transformation of solution to physical coordinates 

The transformation of the solution into the physical coordinates can be 
achieved by simple quadrature. By definition 


dub ( 4 dx 4 (“*) dt ( ") dy- (7) dx, 
ox C \OS] OY) x dx y 


, Tay CV ; ; 
1.e. : ( "i )x 'dxa+X,/(u;v; LC) dy e dy — PY dx. (35) 
o N L Pj Pj 
If x is held constant, then 
q 
lip tev. £ F . 
y / "i f xi | — a (36) 
N u; J F(¢) 
0 
as y = 0 when € = 0, and TQ l. 


Now if ¢) is the stagnation temperature of the undisturbed jet fluid, it 


can be shown that 


t y—l aro 
°= T = 14+-— M}, 


j ~ 





so tha 


Equa 
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so that. from (12), 


T = T,—Y—* M3U%4(7,—T.)U 
7, —Y—* M} X-4F%(0)4(T,—T.)X-4F (0). (37) 


Equation (36) then gives 


if ou , (v6+-C y—1 c3 3 
1 / J 1 Xx: 3 oo _<, wae M2 a ae A 3 4 T_T ' 
“ 5 ah XW Host 7) (2 ) si Niles 
(38) 


oc 


The streamlines of the flow (“’ = constant, R) are therefore given by the 


equations 


y “) ] T.,, X §(3) log(* = , 
AJ \v Lit " V6 RX-3 


3 
- uy( ex , OR 


ISX 


It will be observed that the flow is symmetrical about the z-axis, as 


) R(T,—T,). (39) 


expected, as 


y(R) y(— R). 
Further, when ¢ v6, y © (hence the limits used in the evaluation 
of L), so that as y—> +oo, uO, and the final boundary condition of 


(11) is satisfied. 


Equation (35) also gives the formula for the derivation of v, viz. 


Or C (“ v; ¢ } X 
5) A L 
ie. ; =x (“"} T. (40) 


\s Cu/ey <0 when y > 0, and F’(¢) has opposite sign to ¢, then from 
30), y and ¢ are positive or negative together. Hence as 7’ is positive, 


the transverse velocity in the jet is directed towards the jet axis. 


8. The growth of the jet 

The rate of flow of mass crossing any transverse section of the jet is 
given by 

/ 16 
2] pu dy 2p; | db 2p /(u;v; LC) | X' dl 2 pjx (Gu; v; LC)X}. 
0 ~9 0 (41) 

Hence the jet continually entrains ambient fluid as it proceeds, but the 
constant rate of momentum flux expressed by (23) is maintained as the 


jet will widen progressively. 


5092.18 
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An estimate of the jet width can be obtained if the distance from the 








; , : . In cas 
axis at any cross-section at which vu = }4(u),,;. is calculated. 
; form 
y Y 7 r 1 Cc 
Now uU u;X-*F(l) = u;X- (1-3 " 
) 
and therefore at half-velocity points, 
C = +v3, 
and from (38) the breadth of the jet is proportional to 
7 7#(3\4 « »,/9 5(y—1) - 9 rr? le 7 4b 
b= T,, X4(3)* log(34+-2v2) — 7 M3 + X*v3(%)—T,.). (42) 
9. Results 
The equations (33), (37), and (38) have been used to obtain velocity and 
temperature distributions at a representative cross-section of the jet; these 
are represented graphically in Figs. 1, 2, and 3. The jet gas was taken to 
Uw ] a 
Uai:0| 
' 
be air, and equation (34) used to determine the appropriate values of L. | - 
For other gases, of course, further calculations would have to be made, 
but the work is not prohibitive. The effects of temperature difference and 
varying jet Mach numbers as shown in the figures are self-explanatory. 
Tm 10. T 
rhe constants used were 
C l, In: 
d = 2-5 cm. } small 
Y 1-405. 
and distributions were calculated at a distance 200d downstream of the 4 
nozzle for the cases 
a xe | whe 
(1) ¢, =, 288° A.; M; = 4,1, 2,4 (Fig. 1); ? 
(2) t 288° A.; M; 3; t,t. = 0, 100° A., 200° A., 400° A. (Figs. | 


2 and 3). | s0 th 












LAMINAR MOTION OF A PLANE SYMMETRICAL JET 163 








he In case (1), as TZ, 1, the transformation formula (38) reduces to the 
form 
U; a V6+0\ . y—lapo ce 
xX!,/31 o( — M4 CX*+—C4+-—}. 43 
v|(-¥) V2 19 6 t le a 2 aC aT ( ) 
2) ‘ 
400°A 
=200°A 
0 
id 
se 
tO il i fg 
——— a a aT hE ht CUP 
Fic. 2 
- 400° A 
300° A 
) 200° A 
100°A 
“ 1 2 5 by mn.) ee 
a“ -3 -?2 -1 0 y.(cm 
le, | y¥( 
id Fic. 3 
10. The degenerate case 
In the degenerate case when f t,., and the orifice velocity U; is very 


mall, equation (43) reduces to the approximate relation 


Pai; 
y |("£,) = X#/$log(—>}, 
he 4 “A \v; Le V6—C 
, 2 (yLi /[2u;p; 
vhence : tanh?| / 5 Pj \ 


r 


re 3 7 Oe 
,' that : l—— sech2!Y a Pi \ 
it. ) \2 ai 3Ly; 
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Now, using (34), 


Li 2U; p; 9) | 
y MiPi\) __% 4 where « : af 
xt 2 3 Ly; 3 x 16p; v3 


so that in the degenerate case the solution derived here becomes 


at 
sech?é, where é = - y 
u 3x3 


a 


(44) 


which is identical with that obtained by Bickley (2) for the incompressible 
jet. 


11. Conclusion 


The problem of the plane compressible laminar jet has a closed analytic 


solution, and it is possible to calculate, for any given flow parameters, the , 


velocity distributions, temperature, and density distributions at any 
required point of the flow by means of the equations 


wi 2 
Te Gy iy Semen 45 
u j ( | ( )) 


where ¢ can be obtained in terms of x and y from (38); 


ee 


v 4 px | "i “) a. (46) 


and ee a ry() (1,7). (47) 
J 


Q 


The attention of the author has been drawn recently to a paper by | 


[llingworth (7) in which a similar analytical approach has been made. 
Illingworth assumes a different form of Crocéo’s relation (12), however, 
and although his derived values of u/u, have the same functional form, 
the transformation formula for y (38) is functionally different. No caleu- 
lated results have been given so that no direct comparison with the 


conclusions of the present paper is possible. 
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sible SUMMARY 
[he linearized theory of supersonic flow is used to determine the wave drag on 
slender pointed bodies, both of revolution and with an arbitrary cross-section. A 
mple deduction of von Karman’s expression for this drag is given and two alter- 


vtic itive expressions for the drag are obtained by distinct methods. For a given body 
. the ich of the new forms for the drag is simpler to manipulate than is the case with 
the classical expression due to von Karman. 
any 
1. Introduction 
(45) | THE linearized potential theory for supersonic flow is used to calculate 
the drag on a slender pointed body moving with uniform velocity U 
, parallel to its longitudinal axis through a perfect, inviscid gas. The body 
s taken to be slender in the sense that its transverse dimensions are small 
(46) | compared with its length / and that the normal to the surface of the body 
s everywhere nearly perpendicular to the longitudinal axis. The nose is 
(47) | taken to be pointed so that the shock wave is attached to it and there is no 
stagnation point. The tail may be pointed or consist of a straight cylindrical 
tr by | portion. In the first six sections the cross-section is assumed to be circular, 
nade. | but it is shown thereafter that the results obtained are also valid when the 
ever, » cross-section is irregular in shape. 
form There are three distinct ways of calculating the drag on a body placed 
aleu- | ina fluid flowing with a velocity [ parallel to the axis, namely, 
1 the i) by integrating the resolute of the surface pressure in the axial 


direction : 


li) by calculating the rate of transfer of momentum across some control 
surface entirely surrounding the body ; 
iii) by determining the rate of transfer of momentum across a control 
surface entirely within the body. 
Method (i) was used for bodies of revolution by von Karman and by 
Lighthill (1) who obtained for the drag the expression 
ee l ae 
D = 4p, U? — log ——a (x)S"(£) dadé, (1) 


. . > 


0 0 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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where S(x) is the cross-sectional area at distance x from the vertex, and] 

is the length of the body. Method (ii) was used by Ward (2), who obtained | 
the result (1), and showed that it was also valid for any body, whether of} 
revolution or not. 








- 
pom anne enn a - -- - - ~~ - ----, 1A 
Ax.r) 
‘B 
1 
T(h,o) ' 
s, 
. 0 z 
i ' 
i ' } 
! 3 
' c 
i H 
is, 
: 
Fic. 1. Section through the axis of the } 
body. AOD is the cone x Br; BTC the 
cone «—l = Br. denotes the large 
cylinder. 


In sections 3 and 4 of this paper an alternative form for D is obtained 
by method (ii), and it is shown in section 6 how method (iii) can be used to 
derive still another form for D, namely, 

1 “ 
D boy U2_ - | nQ(n)Q(—n) dn, (2 


a 
« 


0 
where Q(n) is the Fourier transform of S’(x). Finally in section 8 a correc- 
tion is given to the drag for a body whose cross-section is S(x) Lat*sin*7z, 
where ¢ is the maximum diameter of the body. 


2. Basic equations 

Throughout the paper (x, y,z) denote rectangular Cartesian coordinates, 
with the x-axis in the axial direction; the corresponding velocity com- 
ponents are (U+-u,v,w), U being the free stream velocity parallel to the 


a-axis. Introducing a velocity potential ¢ such that w= U(é¢/éa), | 


v = U(éd/ey) and w U(éd/éz), and using the linearized approximation, 
¢ satisfies the usual wave equation, 
~ 
a) 
V26 = M? 
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9? 
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where M is the free-stream Mach number. This equation may also be 


Cc 2h C 24 ‘ 0h 


written 


> 
+ PR, (4) 
cy~ 02° Cx" 
where B? = M*—1 > 0, or expressed in cylindrical polar coordinates 
(r.@, a), 24 J a9 a9 
\ C-@ l ¢ @O l ¢ 24 a CO" ~ 
91 7. —— B? am” (5) 
rer Tr For ox" 


In the case of axial symmetry, (5) simplifies to 


a a (6) 
or“ r or ox" 
3. Expression for the drag using a large surface surrounding the 
body 
The justification and range of validity of the linearized supersonic 
theory lie outside the scope of this paper. In this theory the ‘surface 
pressure’ integral for the drag 


D { (p- Pil, dA, (7) 


where p and p, are the local and free-stream pressures respectively, 1, is 
the cosine of the angle between the normal to the body and its axis, and dA 
is an element of area of the body, is transformed into an integral over a 
control surface consisting of the curved part of a right circular cylinder 
coaxial with the body, together with the flat top and base. 
If the cylinder has a large radius, second-order terms must be retained 
in Bernoulli's equation 
P—Py p, Uu+}p, Bou? —$p,(v?+-w*) + O(9°), (8) 
where p, is the density of the free stream and q? = u?+-v?+-w*. A careful 


investigation is needed to establish the final result that 


D im bo, | (B?u?+v?+w?)(dS,+d8,), (9) 
where the integral is taken over the base of the cylinder at x X.so that 


the effective area of integration is the circle Br < X (Lighthill (3)). For 
convenience this circle is subdivided into the annular region S, lying 
between the cones 2 Br and x—l 3r, and the remainder, which we 
denote by So. 


4. Solution of the basic equation 
A solution of (6) which is valid between the cones x = Brand a—l = Br 
we call d,, and the coordinates of a field point P, namely (x, r), are changed 


to (7, s), where B 


sr «t—S8, (10) 


so that, in the region under consideration, 0 : 
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A solution of (6) is 
dy(x.r) = — 


from which it follows that 


Ob, | f'(@(a—£) dé 


((@—e2— Br (12) 


cor r 
0 
In order to satisfy the boundary condition that the body is a streamline, 
the flow must be tangential to the body. If the radius of the section at 


station x is R(x), this condition implies that 


dR éd, (1 Ody 
dx or ox |? 


or approximately, dR/dx = é¢,/cr on the body, assuming the linearized 
approximation. However. 


1S ; 
; S’(2) 27 RR 27R oP 7 (13) 
dx or 

where @¢,/ér is evaluated at r R. Since R(2) is O(t), to a first approxi- 


1) 


od é 
mation, we can replace R—=! by lim IR 
F 0 cor 


, that is, Roe is approxi- 
oO R—-> cr 


mately equal to f(z). Thus 
J (x) S’(a)/2z7. 


Hence a solution of (6) which satisfies the boundary conditions is 


] *f S’(E) dé 
$,(a, 1) | c = (14) 


) 
=7i 


Lighthill (1) has given a detailed investigation of the order of the terms 
neglected in this expression; the analysis is fairly lengthy and the reader 
is referred to Lighthill’s paper for further details. 

Changing the variable from r to s, we find 


l S'(€) dé 


2n(2x)t J (s—é)! 


[1+ O(«-)] 
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An integration by parts gives 


. *)a | S'(E)(s &)? \, ent | S"(€)(s— £)3 d&é+ O(a-?). 
0 


dh, (x, 8) 


For a body which is pointed at the nose S’(0) 0, whence 
d(x, 8) S”(E)(s—E&)! d€+ O(a-*). (15) 


Since s is related to 2 and r by (10), the velocity components are now 


found to be 


d l 
u <P wh + O(x-1) (16) 
C4 27r( 22°)? 
ad BUI, . , 
und l = In(2x)! + O(a-*), (17) 
where I, : (5) de (18) 
(s—€)? 


We denote the solution of (6), valid inside the cone x—l = Br, by ¢y. 
However, in this case 1 < s < x, so that 


x  — eo 


from which it follows that (22—s) is always large compared with € (which 


is such that 0 < é /). By a reasoning similar to that above, using also 
the fact that S’(/) 0, we can show that 
| "tA - 1 pe 
d(x. 8) S”(&)(s—é)? dé Of (2a s) |. (19) 
7 7(2x—S) 5 ” : 


0 
[t should be noted that this assumption does not necessarily imply that 
the body has a pointed tail and may mean that the body ends with a 
cylindrical portion. Thus, shell-shapes are not excluded. In the present 


case the velocity components are 


Cd, UL, 
72 ; I, + Of (2a—s)-*] (20) 
Cx 2m(2xa—s) 
, Of SUL, 
and US BUI, of O| (2x- s) |. (21) 
a) 2n(2x~—s)? © 


where 
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Bearing in mind that we are assuming axial symmetry we now calculate 


: t 
the drag from (9); we have 
j 
l y We 
—" , (BeU2zR?  BtU272\ 2n(x—s) 
. D = }p,lim | ———* ——.—: ds + and 
ro] y 4? Qa 47? 2x B 
0 
x " 
[ BU BUR aale—6) 
+ - in — ds 
ar m?(2a—s) 477? (2a —s) B? 
i Int 
Pi ~ te J I? ds+ i [2 ds}. 
“x 
Since S(€) is only non-zero for 0 < € < lif 
{ S"(é) dé (23) 
(o—£} . Bu 
0 
it follows th: D =4p,u2 [ Pas 2 
it follows that 3P4 _ as, (24) 
as } Fc 
which is the first alternative form for the drag. 
- eo ‘ . al 
5. Proof of von Karman’s formula using a small cylinder 
Equation (14) gives an expression for d which makes the body a stream- 
line and which is valid at all points in the fluid not on the axis. If we now 
consider the body removed, we may take the control surface to be a small A 
cylinder, radius 7, length /, whose axis is parallel to the axis of the body. | hi 
The rate of transfer of momentum across the curved surface of such a 
cylinder is , 
. {| Ob od 
.U2 | p $ ° Qar dx. 
: Cx or 
0 
. . . . . V 
The density p may, to our approximation, be replaced by p, in this expres- } 
sion; hence ' 
mae * Od Od 7 
D 27p, U? lim | r dx, (25) | 
r0 . Cx or t 
0 


since in the limit the integral over the ends of the cylinder will contribute 
nothing. } 
From (14) we find 


ob f (x—£)S"(t) dé 
fas é)?— Bp 


or 27r . [(a 
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(24) 
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L oy 
so that, as before, im (r=) — — §'(z). (26) 
or 


Ou» 
rv 


aT 


We now change the variable € in (14) toa, by the equation = x— Brcosh yf, 
and we write « = cosh—(x/ Br); we then find 
x 


d(x, 1) . S'(a— Br cosh ys) dip. (27) 


© 


0 


Integrating this equation partially with respect to % we have 


d(x, 7) 5-YS'@ Br cosh ) |, | 
l : , ; 
SS" (x— Br cosh &)(— Br sinh) di. 
Sut S’(0) 0: therefore 
x— Br é 
l r ~r— 
A(a.r xosh-! ——> §"(&) dé. 28 
i= —-_ | oo ree (28) 
0 


For large values of y, 
cosh-! y log 2y+ O(y-*), 

and thus Br 
A( r.7r) a log 


0 


2(x—£ 


) i 
=“ 1 O(r?2)| S°(E ; 
+O ] (€) dé 


An integration from 0 to 2 only introduces another term which is O(r); 


hence 


l i ae . 2 oo 
(x. 7) log(a—£€)S"(€) d+ log —- S"(x)}| +O(r). (29) 
27}. br 
0 
Substituting for r(é¢/ér) from (26) in (25) and integrating once partially 
with respect to x we obtain 
l 


D Py U? lim S'(x)b), T Py U2 lim | S"(x)d dx, 
r—>0 r—0° 
0 


the first term of which vanishes, so that substituting for ¢ from (29) we have 


l 
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) ° 7 Wie 

D f - lim S"(x) da log(a—&)S”(E) dé 
“tt 7 ~{) 

0 0 


9» re 


_P- 


U2... 
1 ~_lim log S"(x)S'(x) dx. 


27 r—0 B 


~ 

















172 G. N. LANCE 
The last term vanishes and the first term is equal to exactly one-half the 
symmetrical double integral 
“as 
= al | S"(x)S"(E)log |a—E€! drdé. 
a8 


This is the result (1) quoted in section 1, namely 


D = 4p, U2 2 | | log ! _8"(x)S"(€) dedé, (1) 
x 


) -- 
=7 


uN 


6. Drag in terms of the Fourier transform of the cross-sectional 
area 
Since ¢(x,r) vanishes when 2 is negative, we can employ the right- 
handed generalized Fourier transform of 4(x,7) with respect to ., namely, 
. 
O(n, 7) | e'ntd (xr) dx, (30) 
0 
in which the imaginary part of n must be greater than some positive 
number c and all subsequent integrals with respect to n are to be taken in 
the half-plane imn > c. For conciseness we write O(n, 7) = d(x.r). (The 
lower dot denotes the transform and the upper dot the original.) It now 
follows from (6) that the equation satisfied by ® is 
72 nd Lr ” +(Bnryod = 0. (31) 
or- cr 
The relevant solution of this equation, which represents outgoing 
waves at infinity, is the Bessel function which has the integral representation 


[ ¢ iBnr cosh 0 dé. 
0 
When ¢ is small this expression is approximately equal to 
log($ Bnr)+y-+ $i, 
where y is Euler’s constant. Hence, if 7 is small, 

O(n,r) = A(logn+logr+C), (32) 
where A(n) is independent of r and C = log(}B)+-y+4ni. Defining Q(n) by 
Q(n) = S’(x), (26) shows that 

, o® 
Q(n) 2a lim {r — 27r(A/r), 
r—0 cr 


whence A 277. de 
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If we define S(a) to be zero outside the range 0 < 2 < I, (25) may be 


written : 
i * Od é : 
D p, U* lim — 2ar us dx. (34) 
r->0 Cx cr 
Parseval’s theorem states that 
ia l = : o 
f(x)g(x) dx F(n)G(—n) dn, (35) 


where f(a) <= F(n) and g(x) = G(n). Using the relation —in®(n,r) = é¢/éx, 


(34) may be written 


ed U2 ¢ . 
D p, U2 S'(x) — da ed | | O(n)inD(—n)]| dn. 


Substituting for ®(n) from (32) and for A from (33), we find that 


ro? t r " 
D (4p, U?) = nQ(n)Q(—n)flog(—n)+ CC"! dn, 
where (’ is a term independent of n. By considering the ranges ren < 0 
and re x 0 separately, it is seen that 
> 1 d j ve S 
D(tp, U? | nQ(n)Q(—n)(log n+ C’ +17) dn+ 


| nQ(n)Q(—n)logn+C" dn, 


. 


9 
= 


which gives the required result, namely, 


. 


D = 4p, l 2 | nQ(n)Q(—n) dn. (36) 


« 


0 
This is the second alternative form for the drag. 
Since (24) and (36) were obtained by two entirely distinct methods, 


it is of interest to prove independently that 


- en e]2 
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By definition Q(n) = S’(x) and 1/(—in)” = a-1/[(v), where v is not 
necessarily a positive integer. The convolution theorem states that 


Q(n) _. [ (w—é/* g, 
Set | se aes 


0 


an integration by parts of the right-hand side shows that 


Qin) ( (x—£)” - 
; ——. § lé. 38 
(—in)’ ° T(v+1) (6) dé ie 
0 
When v = —4, (38) can be written 


nw r gn “\d 
(—in)?Q(n) = a-? | a 


ry 


0 
On writing F(n) = G(n) = (—in)'Q(n) and f(a) = g(x) = m-*T in (35), and 
remembering that J = 0 when 2 is negative, we find that 


= wt], (39) 


2 | (—in)'Q(n)(in)*Q(—n) dn = a7} | I? dx 
2a 
a“ 0 
or } nQ(n)Q(—n) dn = | [2dx - [ | aeat| & (37) 
0 0 0 0 


This completes the identification of the two new expressions for the drag. 


7. Extension to bodies of arbitrary cross-section 
In Fig. 2 Cis a point at a distance x from the vertex O on the line through 
the pointed nose, O, parallel to the free-stream velocity U. The condition 
that the body is slender implies that CP is O(t) for all 0. 
The radius CP, denoted by R, will be a function of x and @ which we may 
write 
R(a, 6) Xo(x)-4 


{x,,(a)cos nO+-B, (x)sin n6}. (40) 


? 


ePA8 


The velocity potential d(x, 7, @ 


— 


, which is also expanded in a Fourier series 
of the form 


d(x, 7,8) = d,(a,7r)-4 (a,, cos mO+-b,, sin mé)¢d,,,(x. 7), (41) 


HhA8 


satisfies (5), from which it follows that ¢,(x,7) satisfies (6) and 4,,(x,7) 
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The drag force is, by (9), 
— ad b < 2 ] 0 2 
D = 4p, U? lim | [B*(= | () +- (3) |r arao 
? ro J J OX or r?\ 06 
0 0 
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. m4 r a 2 
ee »o| CD “ ‘ C 
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ed, ]2 
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al — | 
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Fic. 2. Plane perpendicular to 
Ox at station x. 


When we perform the integration with respect to 0, all the cross-product 


» 


terms of trigonometrical functions vanish and we are left with the result that 
D = mp,U?lim | 


B2(o) (“eo | dr 3 
OX \ or 


ao. U2 » oO 2 é 2 242 
“Pr tim | > (a2, +02 | Br =) +(e) +" er dr. (43) 


2 Ja Ox or r 
0 1 





is 


The first term is just the result (24), obtained in section 4 for a body of 
revolution, which is correct to O(t®log?t) (Ward (2)). 

We now show that a,, and b,, are O(t+) so that the last term in (43) is 
U() and so can be neglected to our approximation. The direction cosines 


{ the outward drawn normal to the body are nearly, from (40), 
lek oR 
=<» —}, 
R 0b ox 
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whence, expressing the fact that the component of velocity of the fluid | 





normal to the body, on the body, must be zero, we have the condition that whie 
l F F 
on the body eb 1 ed eR. * éb\ OR ‘ - 
er R206 08 | ' Ox) Ox 
also 
Hence ” 
+5 = +> @, _cos mO+-b,, sin mé) oo x The 
1 givel 
< [ao(z) LS {ai,(x)cos n+ 8), (a)sin n6}} - 
1 
a — ne 3 OP», 
+ > t,, cos mO+-b,, sin m@) it whos 
. : that 
L > m(—a,, sin md+5,, cos m6), | x TI 
I - 
aSsIS 


: y n{ —«,,(a)sin n6+-B,, (a)cos n0}] | R2(cx, 0) 0. (45) 
1 
For small values of 7 the dominant term of the solution of (42) is 7-”, thus, 
on the body, 4,,(x,7) is O(t-”), so that é¢,,,/ex is O(t-™) and é¢,,,/ér is O(t-™*»), | 
However. by hypothesis, «,(a), 8,,(x), «/,(x), and B)(a) are O(t) for all n, | 1. M 
including n 0. So it can be seen that, since the term a(x) of (45) is O(t), | 
all the terms must be O(t) or higher. In order that this should be so a, 


m m 


te 


and b,, must be O(t”+?). Hence, as explained above, it follows that (24) 3. M 
and (36) give the drag on a slender body, pointed at the nose, whose cross- | 
sectional area at station x is S(z). 


8. Example: Evaluation of the drag for the case S(x) = }xf*sin*xx 
using (36) 
For a body with this cross-section the function Q(n) is given by 
Q(n) bart? | e'"? sin 2ra dx, (46) 
0 
which can be easily evaluated, and we find 
Q(n) Aar3{2(e* — 1)/(n?- ag 
sin? 


hence Q(n)Q(—n) = xt! (ne? Se 


and equation (36) now gives 


~nsin® dn dn 


D = hp, U? 4xXt4 , 
sal . (n?— 47”)? 


























THE DRAG ON SLENDER POINTED BODIES 


> fluid a4 F sin’ 
n that which reduces to D Sp, Ly2 = | a dt: 
(44) , - ; sia 
' mt* ) at? ( sin 2t mt? ( sint 
also Up 7 dt | dt (47) 
+ 3A 0° + J ft 4jJjit 
U0 0 


The numerical value of this expression is 10-99f?. Lighthill, however, has 
siven the result _ 


| ° si t 
, Dart? ie 


whose numerical value is 11-64¢?. The revised value is only slightly above 
that for the symmetrical parabolic profile. 
The author wishes to express his thanks to Professor G. Temple for his 


wsistance in this research. 


(45 
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SUMMARY 
The forward escalator method is a technique for the computation of the latent 
roots and vectors of a matrix in terms of those of a submatrix of any lower order, 
The backward escalator enables the latent roots and vectors of a submatrix to be 
calculated in terms of those of the complete matrix. The properties of partitioned 
matrices are used to establish the basic theorems and equations of these m« thods, 


for both symmetric and unsymmetric matrices. 


1. Introduction 

WITHIN the last ten years J. Morris has introduced, in a series of papers 
and a book (1), a technique for the computation of the latent roots of 
matrices to which he gave the name ‘escalator method’. 


Considering a general matrix of the form 
(A— BA)x = 0, (1.1) 


in which A and B are square matrices, x a column matrix or vector, the 


process is used to compute the roots of the determinantal equation 
A— Br 0), (1.2) 


and the non-trivial solutions x of (1.1) which exist fer the values of A 
derived from (1.2). 

The X’s are the latent roots. the x’s the latent vectors of the matrix (4 Bd). 
In the particular case when B is the unit matrix J it is customary to refer 
to the latent roots and vectors of the matrix A. 

The term ‘escalator method’ derives from the fact that the required 
solution is obtained in steps, the roots and vectors of the complete matrix 
of order n being obtained from a knowledge of those of a leading minor of 
order m (< n), which were previously deduced from those of a minor of 
order less than m, and so on. 

In the method of Morris, which might be called the ‘one-step escalator’, 
the order of the minor is increased by one at each stage. Morris also uses 
the ‘escalator in reverse’, whereby the roots and vectors of a matrix of 
order » are used to determine those of a leading minor of order m < n. 
The terms ‘forward escalator’ and ‘backward escalator’, used in the sequel. 
are self-explanatory. 

N. Aronszajn, in an unpublished article, produced an independent treat- 
ment for the forward escalator only, in which the minor can be increased 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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by any order p, and which might be called the ‘p-step forward escalator’. 
He showed in particular that the two-step escalator was a very practicable 


omputational proposition. 


lee = a . 
The account of the method given by Morris is not easy to follow, the 
equations and results being written in full algebraic detail, and the notions 
f Hilbert space used by Aronszajn are not widely known. The purpose 
ff the present paper is to establish the theory of the general escalator 
lat pay ) g 
order, processes, both forward and backward, by the use of matrices. In the 
to | ourse of this work some interesting and elegant theorems were discovered 
- which seemed to arise only fortuitously in the special cases given by Morris 
ng ; : ; 
nd Aronszajn 
2. General theorems on latent roots and vectors 
apers Consider the matrix equation 
ITS 
1— BA)x 0, (2.1) 
where A and B are real square matrices of order n. The determinantal 
equation . 
1] 1— Br 0 (2.2) 
~ the has n roots. real o1 complex und to each root corresponds a non-trivial 
solution 2, a vector with real or complex elements, of equation (2.1). If x 
13 is any solution any multiple of x is also a solution, and in practice the 
bi vector is usually ‘normalized’ in some predetermined fashion. We assume 
of A throughout that all the roots are distinct, since this is nearly always the 
case in practical problems. 
Bd). The modal matrix X is a square matrix whose n columns are the n 
refel solutions x of equation (2.1). The latter equation can then be written in 
| | 
the form , , ‘ 
ine 1 \ BXA, (2.3) 
wired 
atrix where A is a diagonal matrix whose elements are the distinct latent roots 
or ol } A of (2.2). 
or of lf A and B are unsymmetric matrices we consider together equation 
2.1) and the allied equation 
tol > é 
| B A)y 0), (2.4) 
uses 
ix of the dashes denoting matrix transposition. The latent roots are the same 
n is before, the vectors are different. Associated with (2.4) we have the 
»] matrix equatio1 oar - 
“a is at aaias A'Y = BYYA. (2.5) 
_ Transposing (2.3), postmultiplying by Y, and using (2.5), we find 
X’B'YA AA BY. 
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It follows that X’B’Y is a diagonal matrix, and hence X’A’Y is also | If be 


diagonal. If x,, y, are typical vectors satisfying (2.1) and (2.4), this result posit 
can be written in the form (ii 
2, B'y, = 0=2,A'y, forr #8. (2.7) of th 
—_— , . ss equa 
This is the orthogonal property of latent vectors. We also ‘normalize’ the | 
' ‘ ; e equé 
latent vectors, making them wnit latent vectors, by the equivalent conditions ~s 
take 

a, B’y, = 1, t, Ay, = da, (2.8) 

It follows that xX’ py I ) (2 , 

x'A'Y = AJ : 
From the first of (2.9) we also have (i 
X'B'Y = BYX' = YX'B' = Y'BX = BXY'= XY'B=I. (2.10) | bes 
The matrices XY’, YX’ can be written in the respective forms - 

ru us , ry ud , 
AY’ = > «,¥,, YX" = > 9.4, (2.11) 
r=1 r=1 
and we can then obtain from (2.10) the equations 
n F n : /n ; L , Nn “a } 

(> 2,y,)B = Bl > x,y) =| > y,2;,)B Bi > 9%) = I, Ita 
r=1 r=1 , r=1 r=1 (2.12) , 
results which will be needed later. ¥ 
Special cases gen 

(i) The application of Lagrangian equations in the theory of small 
oscillations leads to equations of the type (2.1), where A and B are If . 


symmetric and usually B is positive definite. In such a case only one } _ Ulli 
equation need be considered, and in equations (2.4)—(2.10) the dashes can 


Ak 
be omitted from A and B, and y and Y throughout are equal respectively - 
° i} 
to x and X. 5 ia 
\U.« 
The latent roots A are then all real. For on 
Ax Bdx, (2.13) 4 " 


x’A x’ BY, (2.14) 
and using a star to denote the conjugate, 
Ax* = Bd*x*. (2.15) 
Premultiplying (2.15) by x’ and using (2.14), we find 
x’ Bx*(A—A*) = 0, 
and since x’ Bx* is real and non-zero, \ = A* and is therefore real. The 
vector x is therefore real, and the root can be obtained in the form 


A = 2 Az/x' Ba. (2.16) 





also 


‘esult 


(2.7 


2.10) 
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If both A and B are positive definite, the roots are therefore real and 

positive. 
(ii) If B is the unit matrix the problem reduces to the determination 
of the latent roots and vectors of A. Then B and B’ are replaced by J in 
Bo). 


equations (2.1), (2.2), and (2.4), and can be omitted altogether from 


equations (2.3) and (2.5)-(2.10). In particular the orthogonality relations 


take the simple forms ’ yes 
| | HYe= 9 (r #8)\ (2.17) 
] (7 s)) 
MW n 
> a y, 2 y, x, 7. (2.18) 
r=] r=1 


iii) In many problems of vibration B will be the unit matrix and A will 
be symmetric and generally positive definite. The vectors x and y are then 
equivalent, the orthogonality relations having the forms 


L_2 0 > s) | 


, (2.19) 
l (¢r s)) 
> «2% =. (2.20) 


It also follows from (i) that the roots are real, and positive if A is positive. 


3. The general forward escalator 
We now consider the application of the forward escalator process to the 
general equation 


(A—BAJE = 0. (3.1) 


If A and B are not symmetric we have to consider simultaneously the 


ullied equation 


(A’— B’r)yn = 0. (3.2) 
\ knowledge is assumed of the m latent roots u, and latent vectors 2, and 
/, of the leading minors of order m of the respective matrices of (3.1) and 
3,2), and from these we deduce expressions for the roots and vectors of the 
complete matrices of order n. 

To this end we partition the matrices as follows: 


mn Pp m P 
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In equations (3.3) and (3.4) A,, and B,, are square matrices of order m, 


A,, and By are square matrices of order p, where n = m-+-p. The column é 
has been split into the columns &,), &) and similarly with 7. 

We already know the roots p, and vectors 2,, y, for r = 1,...,m, given by 

the equations Aut, = pf, Byt, or 2A’, = pt, Bi, (3.5) 

’ y e , , > ‘ » 

Aud, = My Buy, 8 YA = 1 Yr Bu. (3.6) 


The z,, y, are furthermore unit vectors, satisfying equations (2.7), (2.8), and 
(2.12). 


From (3.3) and (3.4) we obtain the four equations 


Ay, & Aj &) A( By y+ By2 & a). (3.7 
Ay, &) A yo &y) A( Bay & a) ! Boy Ea)) (3.8 
Ay, Ha) + Adi Me = (By nay+ Bar 1). (3.9) 
Aj2 Ha) + Ad2 Ne) A( By Nay + Boz 7). (3.10 


Premultiplying (3.7) by y/., (3.9) by a’, and using (3.5) and (3.6), we obtain 


> 


the two equations 


(A L,Y} Bi, €a) Y'(Ayo ABs) 2) 0, (5.11) 
(A—p,)a, By, nqy—2,(As,—ABo,) 2) = 9. (3.12) 


Using (2.12), write equations (3.8) and (3.10) in the forms 


> (Ag,—ABo,)(2, y, By )Eq)+(Acs—ABos Eo = 0, (3.13) 
r=1 


= (Ajp -A Bis)(y,%, Bi) na) ee 
r=1 


292) N(2) 0. (3.14) 


Now eliminating the scalars y’, B,, &) from (3.11) and (3.13), and 2%. Bi, ng 


from (3.12) and (3.14), we obtain the two equations 





* (Ay, —AB,,)2, (A yo—AByo : 
S ( 91 a r Yr A y—AB jo) ie ~X Boo) | Eo) 0. (3.15) 
oot / br 
r=1 


* (A,.—AB,.)y, x,(As,.—ABs ; ; 
\ 12 919)Y,Up(Agy 01) (Abe AB») No) (). (3.16) 
Setting, for convenience. 
(A,,—AB,,)x, = 4, ) 
(A,.—ABj.)y, v, J 


7 


both being vectors with » components, we can write (3.15) and (3.16) in the 
simpler forms 


pee, oe 





. U,U, 4 (Ags ABs) 9 0 
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The brackets contain square matrices of order p, each of which is the 


transpose of the other. The vanishing of the determinant of either of these 


matrices then gives the condition that the vectors &) and 7.) are not null. 
The m+ p n roots of this determinantal equation are the required latent 
roots of the complete matrix (A— BA), and the solutions of the linear 
equations corresponding to (3.18) give, for a particular A, the ratios of the 
» elements of the vectors &, which are parts of the vectors €, 7 finally 
required. The remaining parts ¢ 
3.11) and (3.12 


ition. in the forms 


): 7q) Of these vectors are obtained from 


after premultiplication respectively by x,, y, and sum 


; mn , 
‘2,9 1 AB,,) . ‘ £., 
\ Vi nS Oe 
pad . r A By 
(3.19) 

 y 2'(A’.,—AB,,) 2. a’ 
S\ Ung = > 2 mw 
— A fant A B, 
r=1 r=] 

th being vectors with m components. 

Condition }O unit vectors 


Equations (3.18) and (3.19) give the ratios of the components of the 


mplete vectors € and To make these unit vectors we have to satisfy 
+] (2 

le conditiol Leo ) , ¢< 
Je, é B 1) dE l. 


terms of our partitioned matrices and vectors, this becomes 


, , > , pe , - 9 <¢ 
na) Br a n2) Boy & 11) Bye a) + 2) Boe Fy) l. (3.20) 
Substituting now for 2 »: €4) from (5.19), we have 
m 6 
4 ? " 2, t.tia 
N\ to" By? > ee 
— \ i os A i a ——s A By 
1 ’ 1 
» UY . P m 
4 TS By at 12) Bor &e 1. (3.21) 
hamd A— WU, sei ™ artes 


"7 
‘ 


From ) and the first of (2.8) the first term in (3.21) is capable of 


considerable simplification, and the whole can be written in the form 


ue , 
| iT . . Dean ~ P 
1S i > S 2 Byet Beslfe = 1. (3.22) 
* — \ ‘as — ‘ ee ae b 5 a e 
1 r=) r=1 , 


Now the matrices in (3.18) are functions of A, and to emphasize this we 


can write those equations in the forms 
F (A) Eo) 0 | 
Ff (A)? 


(3.23) 


I(2 
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the dash denoting matrix transposition in the usual way. The determi- 
nantal equation for the latent roots is then 


F(A)| = | F’(A)| = 0. (3.24) 
From (3.22) and (3.17) it easily follows that the new vectors, to be unit 
vectors, must satisfy the condition 


; a eee 


the bracketed quantities being square matrices of order p. 

This rather elegant result can be further extended to determine the 
actual magnitude of the components of &) and yj), and therefore of , 
and 7), necessary to produce unit vectors. Denoting the respective p 
components of &(), 72) by €;. n; (¢ = 1, 2.,..., p) it then follows from (3.25), if 
f,; is the element of F in the ith row and jth column, that 

SY Sing (3.26 
O77 = ds 0.20) 
—_ dxr li 
v =f 
Now if F,; is the cofactor of f;; in the determinant |F|, the matrix [F;;|, 
the adjoint of F, satisfies the equation 
Fadj.F = 0, since |F 0. 
From (3.23) we find, therefore, by postmultiplying by 7’ and €’ respec- 
tively, that , ia 
? Ey = q adj. I | 
n&’ = qadj. F’ J 
where g is any constant. These equations imply that 
¢  , 
fens = Oi \. (3.28) 
g; Ni qF;; 


(3.27) 


It then follows that the value of q which will satisfy (3.26) is given by 


p p d : 


$2, qr sia eis ~, 

t=1j=1 

d 
or =I F| : 3.29) 

17 ( 
Finally, we have the result 
re" i a; 

En = FGF o | (3.30 
This result shows how to choose the components of £), 7.) so that the 


equations (3.23) are satisfied, and the complete vectors te n are unit 
vectors. One element of &,) 07 7) is at our choice, the remaining elements 
are determined without ambiguity. If any corresponding pair of elements 


of &), 7) are made equal, everything is uniquely determined. 
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Spe cial cases 

Special cases can be treated separately by analysis similar to the above, 
or can be deduced directly from these results. 

(i) (A—BA)E = 0, A and B symmetric. Here A,,, Ago, By, and Boys are 
also symmetric, and A,, = Ay,, By = By,. Only the one equation (3.1) 
need be considered, 2 and € being throughout identical with y and 7 
respectively, and wu, is the same as v,. The main results are then as follow: 





F(A)é | NOU Ur 4 (Ags ABys) | Ew 0. (3.31) 
L— A br ‘ | i 
1 
‘ — v, U, i (3.32) 
aH A—_ ps, 
r=} 
u, = (A,,—AB,,)z,, (3.33) 
—) 2 ’ 

E(9 E F(A)\é 9) l (3.34) 
» » F fd F . (i ° 1.2 ) (3 35) 
€,€ y M ¢, 7 ann ?), O.v0e 
wns Aan . ' oi 


where F is the matrix in (3.31), |F| its determinant, and F,; the cofactor 
of the element f,, of |F 

(ii) (A—JA)E = 0, (A’—JA)n = 0, A unsymmetric. Here B, B,,, and By, 
ire unit matrices of order n, m, and p respectively; B,. and B,, are zeros. 
The main results are as follow: 


F(A) Ee F’(A) nH = 9, (3.36) 
where ; _ 
F(A) = |S —t+(4,.—Al)], (3.37) 
——* by J 
r 1 
ie Rae (3.38) 
Ur AinY, 


It will be noted that w,, v, are here independent of A. 

Equations (3.19) for &), nq) in terms of & ), 7), and (3.25) and (3.30) for 
the normalization constants, are unchanged in form. 

(ili) (A—JA)a = 0, A symmetric. Here case (ii) is further simplified by 
the fact that only one equation need be considered; £ and 7, u and v are 
respectively equivalent. The main results are then as follow: 


F()~o = |S ot 4+-(Aso—IA)| Ee) = 0, (3.39) 
hom Br 
r=] 


a, = Ae 


on %> = Mand. (3.40) 


the other results being identical in form with (3.32), (3.34), and (3.35). 
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4. The general backward escalator addit 
Treating the same general equations (3.1) and (3.2), partitioned as in vant 


(3.3) and (3.4), it is now required to calculate the roots and vectors of 
(Ay,—p By, )a = 0, (Ay—phBy)y = 9, (4.1) 


from a knowledge of those of the complete matrix equations (3.1) and 

(3.2), that is from a knowledge of the vectors &, and », partitioned respec- 

tively into &),, €),. and ),. 7),, and of the corresponding roots A,, for Fron 

r i Se n. 
Corresponding to (3.11) and (3.12) we have the two equations 


(u—A,)y’ By, yp t+ y'(Ayo—A, Bix Ei, = 9 ) (4.9) 
(u—A,)x’ By, Mp +2' (Ao:—A, Bai) 1, = 0 


Now equation (2.12) for latent vectors, expressed in terms of the parti- 


tioned matrices gives, among others, the equations a 
In 
” n 
) < c » ¢ { 
> Bu Sar Ne» > Bir Seren 0] \P 
_ = ' (4.3) vect 
vt ne 
> Bitar kart 2 Bar ter ker = 0 calc 
r=] r=1 
: | 
whence, postmultiplying equations (4.2) by 7.),, 9), respectively and | ( 
summing, we obtain the two equations A 
n of u 
. J i 
Y'(Ayo—p By») S Ea H2)r 0 com 
; ee mal 
r=1 ak 
. (4.4) 
” . 
» | ; t0r ' 
v'(Ay,—pBy,) S Marr S(2)r v unit 
one ft “%y) 
r=1 
Following the notation of (3.17), we write 
(A, OF 2 1 u ie 
21 I 21) | (4.5) 
(Aj.—pByo)y v J The 
. . e » 
both being vectors with p components, and equations (4.4) then become t 
n l n l ! 
wn ae att % Mar £ia)r 0. (4.6) 
an pt— A, of hn pt— A, 
r=1 r=] 
The square brackets again contain square matrices of order p, and we can He 
write — ae ” n 
v'| F(p)| 0 u'| k ()]. (4.7) X 
The m latent roots » are then obtained from the condition | t= 


F| = |f’ 0. (4.8) 


For the determination of the latent vectors we need conditions, 
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\dditional to equations (4.3), which make ‘.. n unit vectors. The rele- 


i 
vant ones are given by 
n 
> Bu sar na > Byer Nar Y | 
eis _— (4.9) 
Z By, 1(1)r a > Boy Near Far l | 
l ’ l 


‘rom (4.9) and (4.2) we find expressions for the vectors in the forms 


[1 Vo 
y > ( a 
(4.10) 
_ = 


Sr (2)r 
Lo p—d, 





which the negative signs can clearly be omitted. 
In (4.10) the square brackets contain (mx p) matrices, and v and wu are 


1) matrice The resulting products give (mx 1) matrices, that is the 


vectors x and y have m components. These are the required latent vectors, 


iculable when uw and v have been obtained from (4.7). 


Conditions for unit latent vect 

\s in the case of the forward escalator, only the ratios of the components 
1 and v can be obtained from (4.7), and hence only the ratios of the 
mponents of 2 and y from (4.10). This indeterminacy is removed by 
king w and y unit latent vectors, and as before, we can find expressions 
r the components of w and v which satisfy (4.7), and which make 2 and y 


nit vectors. The last condition implies, from the first of (2.8), the identity 


4 un , " a ¢ wo m 
> aoe ge, > Se |. (4.11) 
fon [L—A in ph A, 


he known vectors €,, 7, are also unit vectors, and this implies, again from 


? 8). the conditions 


of , e , er , a au 
Cr Pir Na &(2)r Bis na €(1)r Bor Narr t+St2yr Boz Nar | (4.12) 
S(r By, l (2), By. E(1)r Boy (2)s 7 &(2)r Bar Ne)s 0) 
Hence, substituting in (4.11) we have 
j } Foal nt ¢’ 
d C(9)> ¢ * U Yo), ¢ » 7 S(2)s U 
N 2 | K2 x NO Mads $02)s UV _ 
a4 = (u—A,)* ian pLr—A inn 92 pt —A, 
1 1 
y ! F n Py n 4! ¢’ n £ ? 
+ U No), é , . re S(9)e U ~U Ne ‘ < Noa)e Ea) V 
x ! 1 By, *\ H(2)s $(2)s pa = )(2)r $(2)r x. S 1(2)s $(2)s a 
— p—Aa ans LA, a sik A. a pk —x, 
l s=1 r=1 s=] 7 
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From (4.6) it follows that the last three terms of this equation vanish, and 
we have, by analogy with (3.25), the theorem 


«| d F)| “| 5 F(n) | lL. (4.13) 
du du 


The elements of wu and v must be chosen to satisfy (4.13). If these 


elements are denoted by u;, v; (i = 1, 2,..., p) we must have, by analogy 
with (3.30), the condition 
f¢ ik 
U;v; F,, I . (4.14) 
du 


The vectors x, y are then unit vectors. 
Special Cases 

Special cases can be treated along the lines of the above analysis, or 
deduced from it. 

(i) (A—BA)zx 0, A and B symmetric. Only one equation need be 
considered, and the results are as follow: 


p2 Sr Sr! y — 0, (4.15) 
pb A 


[F(u)]e 





r=1 rs 
i . £1) | > 
7 ne | u, (4.16) 
> p—A, ; 
r=1 
+ a ore . “ 
U;Uu; = - Fr. I (v, 9 BB cony @)- (4.17) 


(ii) (A—JA)x = 0, (A’—TA)y = 0, A unsymmetric. The forms of equa- 
tions (4.6), (4.10), and (4.14) are unchanged. 

(iii) (A—JA)x = 0, A symmetric. Only one equation is considered. The 
results are identical in form with equations (4.15), (4.16), and (4.17). 


5. Alternative treatment of (A— BA)x = 0 

The general case (A—BA)x = 0 (5.1) 
involves much more computation than the case in which B is the unit 
matrix J. The former can always be reduced to the latter by premultiply- 
ing (5.1) by B-!, obtaining the equation 

(C—IA)x = 0, C = BA. (5.2) 

The roots and vectors of (5.2) are the same as those of (5.1), but the roots 
and vectors of the leading minors are, in general, different. 

Even if A and B are symmetric, C is not generally symmetric, and the 
roots of the minors of (5.2) are not necessarily real. The roots of the minors 
of (5.1) are real, however, and it is possible in this symmetric case to effect 
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a different transformation, resulting in an equation like (5.2) but with 
C symmetric. 
The symmetric matrix B can be expressed as the product of a lower 
triangular matrix J and its transpose L’, in the form 
B= LL’. (5.3) 
Full details of the computation involved have been given by Fox (2). We 


can then easily express (5.1) in the form 


[L-1A(L-1)’—IJA]L’x = 0 (5.4) 
or 
(C—IA)y = 0 
C = L-1A(L-1) }. (5.5) 
y= L'x 


Since A is symmetric, C is essentially symmetric, so that the roots of the 
minors of C are also real. It is indeed easy to see, from the way in which C 
is constructed, that the roots of its minors are equal to those of the minors 
of (5.1). 

The vectors are also equivalent, and unit vectors y of (5.5) produce unit 
vectors x of (5.1). For if y is a unit vector of (5.5), then 


VeYe= 0 (r #8) (5.6) 
] (r 8) J 
From (5.5) it follows that 
2, LL'x, x, Bx, 0 (r~s) | (5.7) 
[ oes J’ 


which is the condition that x is a unit vector of (5.1). 

This treatment of equation (5.1), when A and B are symmetric, is 
superior in every way to that given by equation (5.2), or to the direct 
application of the escalator process to (5.1). The work of forming LZ and 
L-)’ is less than the difference between the escalator process for (A — BA) 
and that for (C—JA), when A, B, and C are symmetric. 


6. Computational note 

In practice the forward escalator process is used to compute, in a series 
of steps of length p, all the latent roots and vectors of a matrix of order n. 
Each stage gives the latent roots and vectors of a leading minor or sub- 
matrix of the given matrix. If we are proceeding from a submatrix of 
order m to a larger one of order m+-p:p, the latent roots of the latter are 
obtained from the condition of the vanishing of a determinant of order p. 
This is more suitable for computation than the condition of the vanishing 
of the determinant of order m+ p. In the determinant of order p each 
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element is a simple function of the required root. Both the function and 
its derivative are easy to compute, and the roots can be calculated (very 
conveniently for symmetric matrices, when the roots are real) by the use 
of Newton’s method. The corresponding unit vectors are then obtained 
from simple formulae of which some terms are already available in the 
determination of the root. 

The backward escalator process has a more limited field of application. 
One of its uses is to determine the effect on the roots and vectors of 
changing some coefficients in the original matrix. Suppose, for example, 
that all the roots and vectors of the original matrix of order n have been 
obtained, but that subsequently some coefficients have been altered in 
p rows and columns. Using the backward escalator, we can calculate the 
roots and vectors of that part of the matrix which has suffered no change, 
and using then the forward escalator we can compute the roots and vectors 
of the complete new matrix. 

Full details and examples of the computational work involved in the 
one-step and two-step forward and backward escalator processes can be 
obtained from the author at the National Physical Laboratory. The theory 
given here, however, applies to any step number and it may be found 
possible and desirable to use a p of three or even more, if the order of the 
given matrix is very large. The evaluation of the type of determinant in- 
volved, for p greater than two, is rather a formidable undertaking on desk 


machines, but may be feasible using high-speed computers. The roots and 


vectors of successive minors are usually of little interest, and fewer of 


these unwanted quantities have to be obtained as p gets larger. 


The work described above has been carried out as part of the research 
programme of the National Physical Laboratory, and this paper is pub- 
lished with the permission of the Director of the laboratory. 
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ALTERNATIVE DERIVATIONS OF FOX’S ESCALATOR 
FORMULAE FOR LATENT ROOTS 


By GEORGE E. FORSYTHE? 
(National Bureau of Standards, Los Angeles) 


Received 6 February 1951] 


SUMMARY 
‘ox’s equations for the forward and backward escalator processes are obtained 


ore consistent use of partit ned matrices. 


[His note is an appendix to the preceding paper by Fox (1). Fox’s 
essential equations for the forward and backward escalators are here 
lerived rather concisely by a more consistent use of partitioned matrices. 
\ typical abridgement is the use of (A,,—AB,,)~' instead of its Schmidt 
expansion 

SS trth 


La p,—A 


With such use of matrices one can express some of Fox’s formulae as 
special cases of more general relations for partitioned matrices. This is 
true in particular of Fox’s elegant but rather mysterious formulae (3.24), 
25). and (4.8), (4.1 3). 1 
This note will not directly shorten the work of the practical computer, 
for in his work he will actually use Fox’s formulae rather than those below. 
It is hoped, however, that the matrix derivations will indirectly help the 


computer by increasing his understanding of the escalator methods. 


Notation 

Fox’s assumptions and notations are adopted. For any matrix /, the 
transpose of M is denoted by 
denoted by M 


when J is square, its determinant is 


The identity matrix of order k is denoted by J. Let C 


be a square matrix of order n m-+-p, which is partitioned as follows: 
( C1,)™ 
( i 11 12 : (1) 
[21 Coo] p 


It is well known that such partitioned matrices can be multiplied as though 
the C.; were numbers, as long as the proper order of factors is preserved. 
The preparation of this note 


earch 
Formulae numbered decimally appear in the preceding paper by Fox. 


sponsored (in part) by the U.S. Office of Naval 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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The forward escalator 
When C,, 4 0, define the square matrix F of order p by the equation 
7 x Y Y Y—l1/7 »\ 
k Cyo—Co Civ’ Cre: (2) 
The following interesting lemma and proof were given by Schur (2) for the 
case mM = p: 
Lemma 1. If |C,,| 40, |Cj = |C,|.|F]. 
Proof. The lemma is proved by taking determinants of all matrices in 
the following readily verifiable identity: 





Y— 7% £77 Y Y—-1. 
| c in | Cn a K Ch io (3) 
Y ¥—j Y Y a ' 
—Cy Ca" Lp 1G Coe 0 F 
LEMMA 2. Assume |C,,| 4 0. A necessary and sufficient condition for the 


existence of non-zero vectors 
m m 
é - ra ”: bee 
E() Pp (2) P 
such that Cé 0, C'n 0 (4) 
is that |F| = 0. When |F 0, let &(), mg) be p-vectors such that 
F&() F'n) 0. 
Then (4) will be satisfied when 
Fi) — CF Cie Ee), (5) 
7) —(C3;")' Cor Ne). (6) 
Proof. The first statement follows from Lemma 1. When |F| = 0, 
premultiplication of both sides of (3) by the row [0, 7/(.)| yields the relation 
[— 12) Co Cy’, 12)]C = [9, 0]. 
The transpose of the last equation shows that the 7 defined by (6) satisfies 
the second of (4). 
The proof of (5) and an alternative proof of (6) follow easily by expansion 
of equations (4) in terms of partitioned matrices and vectors. 
LemMMA 3. When C O(A) = A—AB and F = F(A) is defined by (2), we 
have d F(A) 
dx 


Proof. Premultiplying (3) by [0 J,] and postmultiplying by 


Yn -C7 Ce " 
| Cn Cath F,]B| - a (4) 


Dp 





we get the identity 
C. C C37 
Y Y 11 12 11 12 u 
| C5; C oa i| 7 | | r. (8) 
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Differentiating (8) with respect to A yields dF'/dA as the sum of three matrix 
products, of which only the term in (7) is different from zero. 

To derive the equations for the forward escalator, we let C = A—AB. 
With the assumption (implicit in Fox’s work) that |C(A)| and |C,,(A)| have 


no common zero,t Lemma | states that |C| = 0 for just those values of A 
for which : 

lor F 0, (9) 
where 


F = F(A) = Agg—AByy— (Ay, —ABg)(Ay,—AB,1)-“"(A1p—AB ya). 


Equation (9) is the desired condition (3.24). To see that F is indeed the 
matrix of (3.15) and (3.24), we have only to use the identity 
rae. 
Cy! (A,,—AB,,)4 2 Ae (10) 
r 1 Ur “A 
where yz, are the zeros of |C,,(A)|, z,, y, the corresponding latent vectors 
if C\,, Ch,. This identity is easily verified, for 


ao 2 . i a 
5 Sees Ss B,,2,y, (since A,,2, = p, By, 7,) 
ri A r=1 
| (by equation (2.12) of Fox’s paper). 
The vectors £, 7 defined in Lemma 2 are the latent vectors of C, C’.. When 
expression (10) for Cjj} is substituted into (5) and (6), we derive Fox’s 
expressions (3.19) for these latent vectors. 


’ 


’ —" ; — CFC, , 
By Lemma 2, 7 N2| —Co, Cy, I], while é | 11?) &. When 
A 2 p > ] >(2 


p 
7) is premultiplied by 74) and postmultiplied by & we find that 
: 7 7 dF (Xx) 
Bé N19) (9)- 11 
) al dX li ( ) 


Hence, for the normalization condition 
n BE l (3.20) 


to hold, it is necessary and sufficient that 


, [dF (A) — 
nia) a }f0 L. (3.25) 


The backward escalator 


Let C be a matrix partitioned as in (1). 


Lemma 4. Suppose |C| 4 0. Let 
(7-1 b Dy De = 
Py D2] p 
m Pp 


A necessary and sufficient condition that |C,,| = 0 is that |\T,,| = 0. When 


+ See Vinograde (3) for consideration of the case of common zeros when p = 1. 
5092.18 


Oo 
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le 0, let u, v be p-vectors such that Ty. u [yu = 0. Then x = Tyu, 


ld 


y = [, v are non-zero vectors with the property that C,,2 = Ci,y = 0. 


Proof. Sufficiency. Since CT = I 


nu? 


Ch Tot Ce Te = 0 (12) 

Coy Ty + Coo Te = L, (13) 
If |T5. 0, there is a vector u ~ 0 such that D,,u = 0. By (13), 

Cy Tu = Cyx F 0, 


so that x 4 0. By (12), C,,x = 0. 

That y ~ 0 and C\, y = 0 are proved from equations corresponding to 
(12) and (13) for the matrix product C’T’ = J,. 

Necessity. The necessity proof is the same, with the roles of C and I 
interchanged. 

LEMMA 5. Suppose that C = C(u) = A—pB. Let T = C-, and suppose 
that p is fixed so that |C(u)| ~ 0 but \Tyo(h) 0. Let the vectors u,v, x, y be 
defined as in Lemma 4. Then 


[dT , 
v “ju = y'B,, 2. (14) 
du 
Proof. By differentiating the identity Cl!’ = J, and remembering that 
dC/du B, we see that dP’ /du = TBI. Hence 
dT» 


‘Ra. h.+t &.8,+0, +t, 2. r (15) 


>. 


du 
Since v'T,, = [T,,u = 0, Dau = a, v'TZ, = y’, (14) follows from (15). 
Lemmas 4 and 5 give the basis of the backward escalator. Let 


C A—pB, 


and assume that |C(u)| and |C,,(u)| have no common zero. Then Lemma 4 
states that the zeros of |C,,| are the zeros of |[,,|. Now it follows imme- 
diately from (10), expressed in terms of partitioned matrices, that 


n 


Po(u) = YS Ser Mey (16) 
a A, B 


where A, are the zeros of |C(w)| and &),, 7), are the last p components of 
the latent vectors of C, C’. Since, by (16), T,(u) is the negative of the F(u) 
of (4.7), Lemma 4 gives the equations (4.8) and (4.10) for yu, x, y. Lemma 5 


then states that the orthogonality condition y’B,,2 = 1 is equivalent to 
the condition aC. , 
v == ie E. (17) 

du 


But (17) is Fox’s condition (4.13), since [T,,(u) = — F(y). 


} 
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THE REFLECTION OF A TRANSIENT PULSE BY A 
PARABOLIC CYLINDER AND A PARABOLOID OF 
REVOLUTION 


By W. CHESTER (Dept. of Mathematics, The University, Bristol) 


[Received 29 May 1951] 


SUMMARY 
The reflection of a sound pulse by a parabolic cylinder and a paraboloid of revolu- 
tion is discussed. It is found that the disturbance lies wholly within a cylindrical 


or spherical wave front respectively, across which there is a discontinuity of pressure | 


if the original pulse has a discontinuous front. 

At the common junction of the incident wave front, the reflected front and the 
barrier, the discontinuity reproduces the familiar doubling of pressure obtained in 
regular reflection. In other directions the discontinuity tends to zero inversely as 
the square root of the distance from the focus in the case of the parabolic cylinder, 
and inversely as the distance for the paraboloid of revolution. The excess pressur 
distribution along the boundary is calculated when the incident pulse is a simpk 
step function, and the results are exhibited in Table I and Fig. 3. 


Introduction 

Two-dimensional problems involving the reflection of sound waves by a 
semi-infinite barrier or wedge have received considerable attention since 
the publication of Sommerfeld’s classical discussion of such problems (1, 2). 
In particular Lamb (3) showed that Sommerfeld’s analysis could be 
simplified by the use of solutions of the wave equation expressed in terms 


of parabolic coordinates. A natural extension of such solutions was to | 


reflection problems involving parabolic cylinders and paraboloids of 
revolution, and this Lamb did for an incident harmonic wave train. 

The reflection of a transient sound pulse by a bluff two-dimensional 
body, or any body of revolution, has, however, not previously been con- 
sidered, and an adaptation of Lamb’s analysis is used in this paper to yield 
the corresponding transient solutions. An interesting feature is that a 
direct comparison can be made between the reflection patterns obtained 
with a body of revolution and the corresponding cylindrical barrier. 

While the results apply only to contours of a particular class, it is hoped 
that they will yield some general information as to the relation between 
the diffraction effects induced by two-dimensional and axially symmetric 
barriers. 

Sound pulses only are considered, though the results hold generally for 
transient pulses satisfying the wave equation and the particular boundary 
conditions here considered. 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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Reflection by a parabolic cylinder 
The scale of time is chosen so that the sonic velocity is unity. The 
ncident pulse, in which the pressure is denoted by F(t—zx), approaches 
the cylinder from the left with its front perpendicular to the axis of 
symmetry; along this axis the x-coordinate is measured from an origin 
it the focus (Fig. 1). The front of the pulse arrives at the origin at time 
0. 


) Yo 


F(t-x) 
= rc 








Fic. 1. 


We use the Heaviside operational notation (see 4), in which F(t) and 
ts operational representation f(p) satisfy 


1 [ ef(p) 
27 . Pp 


dp, 


f(p) = p | e”"F(t) dt, 


nd the pressure in the reflected wave satisfies the equation 
5 = pd (1) 


together with the condition that its normal derivative should cancel that 
f the original pulse along the boundary. 
If we seek a solution of (1) of the form 
d e-Pry (2) 
ind transform to parabolic coordinates (£€, 7) defined by 
x £2— 7”, y 2&n (3) 
in terms of polar coordinates (7, @), 
é r? COs 40, n r} sin 30, (4) 
then the equation to be satisfied by w is 


Cus Gu eu Cu , 
1 SM ale 0, (5) 


o&* on* | 
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Such a solution is 
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This is, in fact, the form of so 
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Fic. 2. 


the incident wave is suitabl 
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and the boundary of the cylinder may be written 


1) = No- (6) 
ug only one of the variables clearly exist, 


x 


C | e-**' dz (7) 
7 
Ce-px | e—2p2* dz. (8 
7 


lution required, for the boundary condition 


= —~—{e-Prf(p)} (9) 
aa) 


at 1 = 7, is satisfied by (8) if the constant factor has the value 


2pm f(p) 


(10) 
2Pn — 2 nq | e-2v2* dz 
No 
dp, € P(x+27* 276). (1] 


yf (per) | e-2P* dz 
x) (12) 
2p —2pm | e-2p2* dz 
No 
behaviour of ¢, for large |p| shows that its 
interpretation represents a function which 
is zero for t < 0, implying that ¢ is zero 
outside the region 
t—x—2n?+-2n2 0, (13) 
or r 2n2-+t, (14) 
since 7” t(r—2). 
The reflected wave front is thus part of 
a circular cylinder whose radius increases 
with unit (i.e. sonic) velocity and travels 
along the barrier at a rate equal to that 
of the incident wave front (Fig. 2). 
The interpretation of (12) is not possible 
in terms of elementary functions unless 
y chosen. We can, however, obtain some 
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nformation from the asymptotic behaviour of ¢,. We require the 


following expansions : 





} e-* dz oa ” 57°p : #n>p?...| (15) 
‘ — 
e~2) ] 
a | er (p> 0). (16) 
tp (2pn)° 
Substitution of (16) in (12) shows that 
d, His 10 f(p) (p > 00), (17) 
which implies d,~ 40 R(t) (t +04) (18) 
nd, consequently, 
bw 2 F(t—a 2y?+2n2) {(t—x—2n?+2n%) > 0+}. (19) 


Thus any discontinuity at the front of the incident wave is reproduced at 
the circular front of the reflected wave multiplied by a factor y/n. In any 
fixed direction not parallel to the positive a-axis the discontinuity will tend 
to zero inversely as the square root of the distance from the focus. The 
familiar doubling of pressure which occurs when a sound wave is completely 
reflected at a rigid barrier is, however, also obtained on the parabolic 
vlinder at the junction of the incident and reflected wave fronts, where 
No: 


Substitution of (15) in (12) gives 








d, ~ aon P f(p) (p >). (20) 

| Thus if, for example, F(t) l (t—+0o) (21) 
so that f(p) l (p+ 0), (22) 
then d, ~~ 10 ({—> ox ) (2% ) 

2t)? 
} ] j lo > 2 
f and bm {(t—x— 2n*+ 275) > 20} (24) 
4 t bE og dy - n= 


This last relation shows that, at any fixed point, the pressure due to 
reflection tends to zero inversely as the square root of the time. 
| The above results are illustrated by a particular example in which the 
} pressure in the incident pulse is chosen to give a simple expression for that 


in the reflected pulse. Explicitly 


F(t—2) 


H(t—x) no H(t—2) (25) 
2?(t—a- 


i 


270) 
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where H(t)=1 (t>0)) 
; (26) 
-0 (t<0)) 
This pulse has a discontinuity in pressure of magnitude } at the wave 
front and tends steadily to unity at points far distant from the front. Its 
operational representation is 


x 





f(p) 1—2p7yy ce?" | e~*n?" dz, (27) 
no 
so that a 
d —— 2p e2P% ve re 2p2* dz (28) 
7 
“lo H (t—a—2n?+ 27). (29) 


24(t—2-4 22)! 


The resultant pressure on the surface 7 = yn, is thus simply H(t—<) for 
this particular pulse. 


Refleccion of unit pulse 
When the pressure in the incident pulse is represented by 
F(t—2x) = H(t—z), (30) 
equation (19) shows that there is a discontinuity at the front of the 
reflected wave of magnitude 7/7. For a more detailed study of the pressure 
after reflection we may interpret equation (11) by means of the Mellin 
inversion formula. Since 


f(p) i, (31) 
it follows that 
d=” dp ”" (32) 
7 J _— nn, l iid 
ho 


When (t—a2—2n?+ 27?) > 0, the contour of integration may be trans- 
formed into the two sides of the negative real axis. The following expres- 


sion is then obtained: 


r nr 


stp | e dz 
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where 
(t—2x ‘ 
T J (34) 
2n5 
7 = 
A : (35) 
"o 
In particular, on the contour 7 No: 
, . e-(t+Ir? dr ,; 
= 5 | (36) 


4 9 
0 (4 m)(1 re-"" | ez" iz) + y2e—2r* 


which depends only on the parameter r. 

The pressure variation along the surface of the cylinder arising from 
reflection has been evaluated numerically from equation (36) for 0 <7 < 2. 
It decreases steadily from its value of unity at the wave front and is O(7-!) 
for large r. For the purpose of the computation a table of the function 

: 
e-™ | ev dy 
0 
by Miller and Gordon (5) was used. The results are exhibited in Table I 
ind Fig. 3. 


Paraboloid of revolution 
When there is symmetry about the x-axis, the differential equation to 
be satisfied by ¢ takes the form 
24 2d, ] d 
CO C-®d c a 97 
to y pd, (37) 
ox op” Pp cp 
where p denotes the perpendicular distance from the z-axis. 
Again we derive a solution of the form 
d e-Pry, (38) 
ind write x €°— »?, p 2&n 


so that the equation for u becomes 


( 4pé) " Pu ; f ; tpn) on 0. (39) 
ce NS og en \Yy oF 
A possible solution of this equation, involving only the variable », is 
u A ei(2p7"), (40) 
which gives d Ae-P* ei(2p7n?), (41) 
where el(y) ( edz (42) 


« 


j 
is the exponential integral. 
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This solution also satisfies the boundary condition (see equation (9)) on 
the paraboloid » = 7, if 


4 nb PS (P) - 
7 e~2P "3 — n& p ei(2p72) . 
so that hb = e-Mixt2?-209)4,,., (44) 
2 mp2? —75) e} 2 2\¢ 
where ‘a — f(r) (45) 


. e~ 2PM — v2 ei(2p7z) 
The interpretation of ¢ represents a wave with a spherical front defined 
by t—72x 27? } 2% 0, 
or r 2ne+t, (46) 
where ¢ is the radial distance from the focus. 
With the help of the following expansions (4) 


ei(y) logy—y+y... (y>0) (47) 
l ] 
mw evUi_-—__—...] (y>o), (48) 
yy 
we find that by ~ 10 f(p) (p—> oo). (49) 
Un 
Thus d~ “No P(i—zx 2y?4+-2y2) {(t—ax 277-4 272) > 0+}, (50) 


which shows that a discontinuity at the incident wave front is reproduced 
at the reflected front multiplied by a factor 73/7. 
When p is small we find, with the help of (47), that 


b,~ —nspf(p)logp (p>). (51) 
If, as in the case of the cylinder, 
f(p)~ 1 (p->9) (52) 
equation (51) implies that 
d~ 70 , {(t—a—2n?+ 272) — oo}. (53) 
t—x—27*+ 275 


An example which can be evaluated explicitly, analogous to (25) for the 
parabolic cylinder, is 


F(t—a) = H(t—2) — x) 


o..3° 
“=p 


I (p) l 2pnj e27% dz. (55) 
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THE REFLECTION OF A TRANSIENT PULSE 
so that 
a r a 2p2? 
d 2p e=! ae dz (56) 
a H(t—a 27? 279) (57) 
t—x+- 279 
Reflection of unit pulse 
For an incident unit pulse, 
F(t) = f(p) a (58) 


the transient representation of the reflected pulse is 


: ev r [ — dz sie 
ne f[ : l re +)) ei(H2A) dA 
i a 4a 1— dre A ei(A) 

? e—2Pn 2p7 - dz Y 





' t—zx 
where, as before 7 ‘ T =, 
Y) 270 

If the contour of integration is deformed into the two sides of the 


negative real axis, equation (59) becomes 


ei 1+-lre-{Ei(n?r)— Ei(r)! 
rap) e-"nr+1 _ a \ Y ) — Me dr. (60) 
2. 11 — 4re-" Ei(r)}?+- 427? r?e-*" 
0 
where Ei(x) P ei(—z). 
In particular, on the boundary of the paraboloid (7 1), this reduces to 
l g ¢ (T+1)r dr ' 
p | N° 9 9.9 oy" (61) 
2 J {1—4re— Ei(r)}*+ 42?rte-*” 
0 


This has been computed for 0 < 7 < 2 and the results are exhibited in 
Table I and Fig. 3. There is an initial unit discontinuity at the wave front 
7 = 0) and, for large 7, 6 = O(7-). 

Fig. 3 shows the comparative pressures at points on the parabolic 
cylinder and paraboloid of revolution for corresponding values of 7. At 
ny given time the abscissa represents the distance from the wave front 
measured along the x-axis, the unit of length being the latus rectum of the 


generating parabola. 
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TABLE I 


Excess pressure distribution along the contour 


for an incident wave F(t—2) 


- H(t—2). 





Paraboloid of 








T Parabolic cylinder | Revolution 
° I I 
orl 0°9092 0°8267 
o'2 078338 0°6959 
o°4 07166 o°5171 
06 0°6305 0°4049 
o'8 05650 0°3305 
are) 0°5138 0°2786 
I°2 0°4728 0°2409 
I*4 0°4392 o°2125 
16 O'4113 o*1904 
1°3 0° 3876 0°1727 
2°0 0° 3674 0°1583 
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wile: ~wraboloid of 
revolution 
4 ry 1 4. — * i 4 L 4. 
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. 
Fic. 3. Excess pressure for incident pulse F(t—2) H(t—z2). 

Of particular interest is the fact that the disturbance effect consists 
wholly of a cylindrical wave (or spherical wave in the case of the paraboloid) 
across which there is a discontinuity of pressure of magnitude 7/7 (or 
ng/n*). This should be compared with the well-known reflection pattern 
produced by a two-dimensional plane barrier, which consists of a plane 
reflected wave reproducing the discontinuities of the original pulse, 
together with a diffracted wave maintaining a continuity of pressure 
across its circular front (see, for example, Friedlander (6) and Chester (7)). 
The parabolic contour is such that the reflected and diffracted wave 
both have the same circular front with a discontinuity in pressure if the 
original wave has a discontinuity. This preserves the familiar doubling of 
the pressure at the point of contact of the wave-front with the barrier, but 
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in a fixed direction not parallel to the positive x-axis the discontinuity 
tends to zero inversely as the square root of the distance of the wave front 
from the focus in the case of the cylinder, and inversely as the distance 
for the body of revolution. 

In conclusion, I should like to thank Miss Nash for her assistance with 


the computation. 
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A BOUNDARY PROBLEM IN ORTHOTROPIC 
GENERALIZED PLANE STRESS 
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SUMMARY 
The results given represent a first step in the solution of certain boundary-value 
problems for aeolotropic plates. The effects of a curvilinear triangular hole on other- 
wise undisturbed stress distributions are considered. Approximate formulae are 
given for the stresses on and near the boundary of the hole, together with numerical 


solutions for certain hoop stresses for typical specimens of oak and spruce. 


1. Introduction 

A CONSIDERABLE extension of the theory of certain boundary-value 
problems of plane stress was made by Livens and Morris (1). They dealt 
with the effect of a single hole, both stressed and unstressed, on an other- 
wise uniform distribution of stress in an infinite aeolotropic plate, and gave 
detailed solutions of a number of typical problems when the hole is an 
ellipse. Assuming that the material of the plate has two directions of 
symmetry in its plane, and that these directions are taken to be parallel 
to the axes of x and y, the boundary of the hole considered is given by 
putting 7 = 0 in the net defined by the conformal transformation 


Zz 2(C) f ig| > a, ers) (C £+1n), (1.1) 
r=0 
where z= a+. (1.2) 


This transformation maps the area outside the hole, and extending to 
infinity, on to the space outside the unit circle z = a,e-*€. Livens and 
Morris point out that for n > 3, the methods given for the elliptic hole 
are not immediately applicable, since the potential functions involved in 
the solution, defined in terms of their variable z,, 


They suggest that this difficulty may be overcome by a modification in 


(t), become non-uniform. 


the method of approach, which will give rise to potentials which are 
uniform functions. 

2. This paper gives a first approach to the problem of the effect of a 
hole, whose shape is a curvilinear triangle, on certain uniform distributions 
of stress. This is the case when n = 3 in the transformation already quoted. 
The modified method suggested by Livens and Morris is employed. It will 
be seen from the analysis that approximations have to be made, but a good 
degree of accuracy is achieved in the final results. The stresses found are, 
owing to the methods used, valid only on the boundary of the hole, and in 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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its neighbourhood. Since the resulting expressions for the stresses are 
necessarily complicated, numerical results are given for several cases, 
using values of certain elastic constants for oak and spruce given by 
Green (3). 





3. The boundary curve 


We take as the boundary of the hole the curve given by 
ce-*€(14 per’), (3.1) 
where Cc c,e (3.2) 


ind Cy, wz are real. 

The boundary is thus a curvilinear triangle, with rounded ‘vertices’ 
curring where cos 3 |. The mid-points of the sides occur where 
Is BE 1, and the curve is flat at these points if » = }. We shall 
issume that « has this value for the purpose of approximation and the 
numerical evaluation of the results. The complex form of ¢ allows for any 
rientation of the hole in relation to the stress axes of the material. 

We may write the transformation 


z = ct(1+pl-), t = e-if, (3.3) 


so that the singular points fall inside the circle |t} = 1 when p < 3}, and 






this condition is satisfied if we take p = }. 
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4. The fundamental equations of the theory 

We use ££, 77, &) to denote the stress components appropriate to 
the (€,) coordinates. The complex stress combinations 22, 22 are given 
in terms of a potential function Q by means of the equations 


2z = ea+yy = +77 = —, (4.1) 
C2 
” — - —~ i eQ 
22 = rx—yy+ zy — (€E— n+ 2iéy)er? — ——, (4.2) 
a2 


where ¢ is the angle the outward normal to the curve 7 constant 
makes with the x-axis. 


The potential corresponding to a given undisturbed state of stress, Q,, 


may be written down, using equations (4.1), (4.2), since zz, 22 are known. 
It may be shown that on the curve 7 = constant 
0Q) 


2(q7— itn) e (4.3) 


We now consider two cases. 

1. For the unstressed hole, 60/eé is zero on the boundary, and so Q is 
a constant, which may be taken to be zero without loss of generality. In 
this case, in order to obtain Q, we have to add to Q, suitable terms which 
cancel it on the boundary and represent stresses which remain regular over 
the plate and tend to zero at infinity. 

2. For a stressed boundary, where (77) i€n) f(€), 

f 
Q = 2 | f(u)z’(u) du. (4.4) 


Here we need to add terms to Q, which will give the correct value of Q on 
the boundary. 

These additional terms are functions of a variable ¢,,,(z,,), which will now 
be defined; Q may be written 


—— : 
2 = 24(2,)+Q2(%) + 5 O44) + 5-44), (4.5) 


where Z. = 2+A,2 (4.6) 
and A,, A, are the y,, y. of Green’s analysis, and are constants less than 
unity for any given material. Consider the equation 


im =1t( > a,t7)+42,,t-( 5 ar), (4.1) 
r=0 r=0 


where ¢ = exp(—7é) when z,, corresponds to a point on the boundary of 
the hole. If we can solve the equation for ¢ in the form 


t= exp(—7€) tm(Zm)s (4.8) 
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we have a function ¢,, which reduces to exp(—7z£) on the boundary of the 
hole. It is this function which Livens and Morris use in their solution, 
interpreting both Q and Q, in terms of it. For n > 2, however, ¢,,, even 
if we can isolate it, is not a suitable function, because branch points of the 
;,-t transformation occur when |t| > 1, and any potential Q, a function 
of t,,, Will not, in general, be single-valued throughout the plane. 


5. Modification of the method for n > 2 


For n > 2, we cannot isolate the roots of equation (4.7) in simple form. 
7 


We know, however, that when the values of Q in (4.7) run from 0 to n+1, 
there will be » roots such that |t) > 1, when z,, corresponds to a point 


m 
' outside the hole, and they tend to infinity with z,,. There will be just one 


root ¢,, such that |¢ | when z,, corresponds to a point on the 


m0 mo m 


boundary curve. It is the reciprocal of t,,,. which we use as the basis of our 


modified method (cf. the details for the ellipse given by Livens and Morris). 


For » 3, the relation between z,, and ¢ is 


z,,(t) ct(1-+-pt 3) 4 r,, ét A] +- ut) (5.1) 


and we have two roots outside |{¢ 1, when z,, corresponds to a point 
t 


uniform functions, in combinations of the form 


outside the hole. If these roots are ¢ we use both of them to get 


m0» “m1? 


| l 
fs u 


) 
mo ml 


t 
the combinations used give functions which tend to zero at infinity. We 


where p is a positive integer. Then, since t,,, t,,; tend to infinity with z,,,, 
interpret each of these combinations in terms of t,,.) only, so that the value 
of the function on the boundary may also be expressed in terms of t,,9 


| which reduces to f e—S there. 


6. Expressions for the uniform functions 


We may write 


=™ 


| | l F z(t : 
(f) — dt, (6.1) 
} p } Dn ine o 
tno t l aUe . u im(e) Zm\tmo)s 
where C is a contour consisting of the unit circle || 1, the circle |t R, 
R being very large, and a radial line joining them, as shown in Fig. 2. 
If p is a positive integer, there is no contribution from the integral round 


| the circle at infinity, and we have to evaluate the integral for |¢ 1. 


5092.18 


P 
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. 
| 


3 l l | l 
anil r x) pp O82 m(t)—2m(tmo) T Pp | pai OBi=(!) Zt no)s dt T 


t 1 t 1 


” 


7 5(1/t)—Z(1/ting)] | 
P | poi 8)! I rn 


dt. (6.2 
2(t)—2(tno) | | 





Denoting the three expressions on the right of equation (6.2) by J, 1, J, 
respectively, we can show that 


l : 
i Q, f, . (6,3) 
y tno Pp 
We then have to evaluate 
F | Cult t-+-t 1} 

L=p log! LX, Mmolt + tmo) TT lat. (6.4) 

. ) . tla f 2 2 

“1° 1 . Ct nor | p(t T tino) nol s 
Now on |t 1, z(t), Z(1/t) are of the same order of magnitude. Also 
Xn 1. So we may expand the logarithms in powers of A,, for certain 


values of t,,9. It should be noted, however, that the presence of ¢,,,, in this 


expansion means that it cannot be valid when ¢,,,., becomes large. So the 
results obtained from this method of evaluation of the integral will give 
valid stresses only on the boundary of the hole, and throughout a finite 
area surrounding it. Since, however, it is the boundary values of the 
stresses which are of the greatest practical value, this limitation does not 
invalidate the method. 

We have 
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ve now pick out the terms in /, which, for a given value of p, involve 


i before integration. All other terms disappear on integration. Clearly 


eshall have to stop the expansion in powers of A,, at a certain stage, and 


1, 2, 3, 4, 5 correct to A> 


eshall give the results for p x 


7. We obtain the following results: 
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where 
66° .- 364 
(3) a 5 5 (3) es 4 4 
Am2 = —— AmB ’ an = “ct Xn 
(3 C5 (3) c° 4 6 
3) ‘ 3 p 5 (2,,4__1@ . "pf 
Amo —aAm ; Gm,-1 co Am( 3H 18u ) fs (7.6) 
3 , (3 (3) 
(3) é 9,5 3) : an hen te 
Am,—2 Xn “a OP 5 Am,—3 1, Am, ~~ Op. Xm co 
J 
] 1 1 
! a 7 (4) 4q¢ si 
(4) 14 TA = ££ “tM tno: (/.7) 
mo m1 qd 4 
where 
(4) 56 (4) got (4) 
vo Ps . 4 
Gm, thn GH ’; Amo Xm nal ’ Am,-4 I 
\ (7.8) 
@ | 
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For p > 5 we get, to our approximation, 


I z.. L (p) p 
yp | Dp An; ptmo: 
mo m1 


where a’? _, 1. 
8. Solution for unstressed boundary of the hole 

Simple tension and shear at infinity 

The uniform and undisturbed distribution of stress is assumed to 
arise from a simple tension 7' in the direction of the 2-axis and a uniform 
shear S parallel to the axes. 

These correspond to the stress components 


LX i yy 0; ry S. (8.1 


The conditions for the corresponding potential Q, are therefore 


Od eQ 
—*=T; aa! = TA Oi8: (8.2 
Cz Oz 
so that the appropriate form for Q, is 
Q Tz—(T+-2i8)z. (8.3) } 
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4s an illustration of the method of solution, we shall give the general 
xspression for 2 correct to A>,, and show how to evaluate the coefficients 
nit before making a further approximation for numerical purposes. 
(7.6 We have 
Q = Tz—-(T+28)z- 
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the boundary of the hole, t,,. = ¢ = t,¢, and so we get aset of equations 
for the coefficients A,,,, A,,o,---; Ams, by equating to zero the coefficient of 


every power of ¢ on the boundary. When this is done it is found that 


ed to | Ans: Amg are known as soon as A,,,, A,,2 are determined, and we may 

niforn | find A,,,, A,,o from the four equations given by the coefficients of ¢, ¢-1, 
| {2 { 2 

9. Before proceeding with the algebra of this solution we consider a 

(8.1) | further approximation. A glance at the coefficients a‘>), shows that powers 


} of u higher than the fifth occur in many of them. If we are neglecting 

terms of higher order than A°,, then clearly we need not, for numerical 
(8,2) | purposes, retain all the powers of w given. As a first approximation 
we shall obtain the solution correct to the fifth power of A,, and » com- 
bined. This apparently drastic approximation, for numerical purposes, is 
ustified, since it may be shown that additional powers of » have very little 


effect on the values of the hoop stress. 
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10. Let us consider the four equations derived by equating to zero the | 
coefficients of ¢, ¢~1, #, t-* in equation (8.4). To the order considered, and 
writing a = é/c, the equation from the coefficient of ¢ may be written 


MHP 4 Ape 


9 2 P 
a* a* 


Ao + (; + Az ay?) + ; As a'y?) A,, B,, -(10.1 
1 2 
where B, Tc. (10.2 


The equation from the coefficient of ¢-! is 


VB pw? 3 py? > 3 P 
(! a Jan 1 () a ee Pap amp? Ay, +AZa*py?A,, B,, (10.3 


ae ae 


where b, (T+ 2iS8)é+ 2u3a3(A? Mes 1-A2 Aas). (10.4 


1 


ton 


Taking the conjugates of equations (10.2), (10.3), we obtain two further 
equations 


1 Ay? L Asp? et ee ee 
(; + = Jan (; I af Jaa + Ay wPa* Ay) + Az pra? Ay, B;, (10.5) 
“4 “Wo 
r. w2 Neu? ee a 
ue Ay+ af Ay, +(1+Aj pa) Ay, +(14 Az wa) As, B,, (10.6 
where B, Te 


| aot 


9,3 
B, = (T—2i8)e+ as (Aj Aj2+4 V3 Ags) | 


Equations (10.1) to (10.4) may now be solved for A,,, Ao, and their 
conjugates. The solution, to the correct order, is 


(Ap—Ay)AAy = a2u?{ Ted? Ay + (T —2iS)c(AZ-++A, Ap—ABA,—A, AS—DZAB)} 


Ted, A,—(T-+ 2i8)e)A ! u{ = 8 } aun] — 
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2u3a5\3 , ~2u3a3A, A3 A oo (10.8) 
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where A f u(r 3) xi). (10.10) 
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Now consider the equation from the coefficient of t. Writing 


it is 


d,a23A*A,, +A, 023A*A,, + — A*A, + 


Similarly, the coefficient of ¢-* gives 


Mu” . ae —aere — 
_oF A*A,, a A*A,, +Aj p3a2A* A, +A3 p3a2A 
a°’ a 
A* Ay, I A* Ay, 


On substitution for A,,, As), 


(10.12), (10.13) the equations 


(10.11) 


*4,4 


TépA*. (10.13) 


and their conjugates, we obtain from 


(T’—2iS8)ca2y3A, A,(A,-+A.), (10.14) 


A, AA A; AAgo (T+2218 Epa, Ay 
3 
AA,.+- AAg, (T'+-2iS)é(Aj+A, A. A3)— Te 
a° 
pa*(T'—2i8)e(Az 


giving the solution for A,,, As» 


11. Thus, finally, we obtain the coefficients A,,, As, 





g 








, Ay, Ago. They are 








A*A,, a**(T' — 2iS)c(Aj+A, A, A? A, A, AF AR) a*u?T rj A,- - 
=— - ee = l — 
TX, A.— (T+ 2i8)é}A,4 p(o 3) f a>u?At 
2 a 
A*A,, a?.?(T' — 28 )e(Ay AS+-Ag AZ+-AZAZ—AZ—A, A.) —a2u? Tc, AZ- 
, Rte — 
TCA, A, t (T 21S)é|A,-4 pe ae 4 sas t a3u?rg 
. 2 a ‘ 
aA, (T'— 2i8)e[ pA, A,— 3a, (AZ+-A, A.+A§Z)]+ Tepa, 
ms _A3 3 
(7 T US )E ih 
as 
A*A (7 278 )e| wad, (AF A, Ao +-A¥) pA, Ay] T cpr, 
3.3 
-(T-+ 2i.8)e net 
a 
(11.1) 
To the order considered, we find that all the remaining coefficients 
A)g,.... A. are zero. 
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We thus obtain, to the correct order 


Q = Tz—(T+2i8)z4+ 


9 
- (1) 1 qi) 14 q@) 4 gq) | 
> Amlay atmo Im, 1!mo ! Ano | Ani tmo| ; 
m=1 
9 in 
> : A fq) 7 7-1 1) (1) 7 
1 m = = =/( i= 
a d [a m, tnd I a,- atmo T Got Gn tmo| T 
m=1 ~ 


i eg (2) ‘ 2 
[ow A nal ayy ating 4 a2 
m=1 
9 
‘ lit =(2 2 
! ra = ra 9 
ic oe [ae tmo T aol. (11.2) 


12. Determination of stresses on and near the boundary 


We have noted before the stress combinations 


LL YY $o+97 =- - | 
ij (12.1) 

xxe—yy + iay —(g€—n7- En)erid > —— 

oz 
We may notice also that 

Ano l Ot mo —_ l (12 2) 

Oz Am oz e(J - 2utnd) I Xn e(2 Dut mo— - tno) 
= mo ng z e —=—5 ] = (12.3) 

Xm Oz Cz C( l- 2ptng )4 Xn c(2 7 foo) 


-_ 


To the order considered we find that the stress combination £¢+-77 is 
given by 


uv 
uo 


el T.9, (1)__ (1) 2 9q(1) 3) “10 
7H = re T+2A,,(aQ—al_, ti 20;,2 bi -" + 


Ct. 
2A, (aW—a_, tip? —2a_, t;,3) “2 
Cz 
— . , ot 
-44,, a _.t,° 7 —14,4" _,t,,* a) (12.4) 
Oz wate eee 


(where only terms to the fifth order in A,, and » are taken). The explicit 
expression for this stress is extremely complicated, even on the boundary. 
We shall not give this explicit expression, but shall content ourselves with 
expressions for 


(1) the boundary hoop stress in the isotropic case, 
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BOUNDARY 
(2) a numerical evaluation of this hoop stress at certain points for 
typical specimens of oak and spruce, for a particular configuration. 


(1) For the isotropic case, A, > A, > 0, and we find the hoop stress on 


the boundary given by 


(7 (1—4p2)—2[T cos(2¢—20)— 2S sin 2(—8)][ 1 — 2 cos 3€]— 
- 2| 7 sin 2(€—6)+-28 cos 2(€ —6) |sin 3é 
, 1+ 4y2—4y cos 3& 





> 
(12.5) 
where we have now written 
ei0 


C Co 


For the case in which S = 0, and the plate is under a simple tension 7' in 


the direction of the x-axis, we obtain 


S T (1—4p?)—2T cos 2(€—6)(1— 2 cos 3€)— 27 sin 2(€—@)sin 3 


t&), , a 2 _ 
. 1—4y cos 3€+- 4.2 
(12.6) 
so that at the particular vertex € = 0, this stress is 
> T'(1+-2n)—2T cos 20 oa 
rz: H) (12.7) 


1—2p 


Equation (12.7) is Stevenson’s (2) result, remembering that his (7, &) 
correspond to our (—€, 7), and that his 8 should be replaced by @ for this 
result. 

(2) Green (3) considers specimens of oak and spruce, for which ,, A, are 


given below: 


Spruce . . A, = 0-608 A, = 0-111 
Oak . . Ay = 0-395 A, = 0-026 


Taking these values of A, and A, we give below the values of the hoop stress 
at certain points on the boundary curve. These are the three ‘vertices’ 
A, B, C and the mid-points D, E, F of the three ‘sides’. We take p = } 
and 6 = 0, so that the x-axis is parallel to the line joining the vertex € = 0 
to the mid-point € = 7, and we take also S = 0. The results are thus those 


for simple tension 7' in the direction of the x-axis. 


Table of values of €€/T at points of boundary (see Fig. 1) 
— en 
F,E | BC | D 


4°33 
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13. All-round tension at infinity 


The undisturbed distribution of stress now arises from stress components 


re = T = yy; xy = 0. (13.1) 
These give éé mn = T; En 0, (13.2) 


so that we have an all-round tension 7' along any boundary curve. 
The conditions for the corresponding potential Q, are 


€Qo 2Qy 


~“0 _ oT, ~° _ 0, (13.3) 
Cz Ce 
so that the appropriate form for the potential Q, is 
Q, 2T 2. (13.4) 
Using the same method as before, and equating to zero the coefficients of 
t, ¢~1, #2, t-? on the boundary, we have four equations: 


A f 
IT ¢ 4 a) | 4 a) 4 11 aa "Wipe 21 a | 
- 4441 41 1 449 9) i. —2 2-1 
Mi di | 
7 i 7 | 
(1) (1) <* Ri (2) 21 (1) “*12 =(2) **22 =(2) | 
0 = Ay ay? 1 +Ag ay” 1+ ie ae i a i 
1 2 1 a | 
0 An aa) | Ay, am) 4 Aja x2) Ayo =) | 
2s 2,-271 Le 2,—2 
A, A, A, A, 
- A A, ' " 
2T cu Ay, @)) 9+ Ag ap 9+ a+ aQ+A a? .+Asga? 
r, r, } 


(13.5) 
We proceed to solve in exactly the same way as in the previous problem 
and finally obtain, to the order required, 
(A,- A,)A oo 2T cprs, (A, r,)Aj. 2T cpa, 
(A.—A,)AA,g, = 2T CA, A,—2T ca*pAr3 d, 
(A,—A,)AA,, = 27 ca®u*r? A,— 2T CA, A, 


the remaining constants being zero to the correct order. We may now, in 


(13.6) 


theory, evaluate the stresses on the boundary, and at points near it. 
For the isotropic case, we find the hoop stress on the boundary is given 
by 7 2\ 
‘ = 27 (1—4p?) 
(E& \n 0 


; 1—4y cos 3+ 42’ 


UT, 


(13.7) 


which is Stevenson’s result. 


For the same samples of spruce and oak the results are given below 


Table of values of éé T at points of boundary (see Fig. 1) 
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14. Solution for hole under constant boundary stress 
We assume a constant pressure P and a constant shear Q to be applied 
to the boundary of the hole, the plate being under zero stress at infinity. 


0. —~ — 


So when 7 m7 P. En Q. (14.1) 


We write mn —vé1 P—iQ = —py. (14.2) 


Then on the boundary of the hole we have from (4.4) 


. 2 ( Poz' (u) du, (14.3) 


so that Q) 


‘ 2p )c(t+-pt-*). (14.4) 


Our potential 2 must thus reduce to this value on the boundary and tend 
to zero at infinity. We now assume for 2 the expression given in equation 
8.4), but omitting the terms in z and Z. On writing down the boundary 
conditions, we find that the equations obtained are those of the previous 
problem of all-round stress, with 7' replaced by py. The solutions for 
Ay, A jp, 
combinations may be evaluated at any point on or near the boundary of 


the hole. 


A,,, Ao. therefore follow immediately, and, in theory, the stress 


15. Hoop stress on the boundary 
We shall give two results:+ 
(1) the value of the hoop stress in the isotropic case, 
(2) numerical results for samples of spruce and oak for py = py = P, 
}, and @ 0. 
(1) The isotropic case. Here A, A, 0 and we find the hoop stress 


given by ; _ 
oF. SufP cos 3€—Q sin 3—2pP} 


(££),-o = 
| —4py cos 3€+ 4:2 (15.1) 


which is Stevenson’s result. 
(2) For the configuration 6 0, p +, De = Be P, we find the 
following numerical results for spruce and oak, given below. 


Table of values of €&/P at points of boundary (see Fig. 1) 


Considerations of space have led us to omit a third result, an explicit expression 


tor the hoop stress, to the order considered, for Po Po P, and Po Po d. 
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THE STABILITY OF A COMPRESSED ELASTIC RING 
AND OF A FLEXIBLE HEAVY STRUCTURE SPREAD 
BY A SYSTEM OF ELASTIC RINGS 
By G. W. H. STEVENS (Air Ministry, London, S.W.1) 
Received 30 November 1950] 


SUMMARY 
The stability of a compressed elastic ring has been studied by a method which 
in be extended to solve the problem of the stability of a flexible heavy structure 
spread by a system of hoops as in a crinoline skirt. The original work by Levy, which 
was developed by Timoshenko and Love, cannot be generalized to problems in 
which the compressing forces are affected by the deformation of the ring. 

It is shown that the load at which a ring will buckle depends not only upon the 
magnitude of the load but also upon its first derivative relative to the radial distance. 
\ positive derivative causes the ring to buckle at a higher load. When this result 
s applied to a cone of heavy and loosely draped fabric spread by a rigid hoop of 


radius 7, and a larger and flexible hoop of radius 7, below it, both hoops being in 


horizontal planes, then various modes of buckling other than oval are possible 
according to the relative magnitudes of r, and r,. It is found that oval buckling 
hanges to three-wave buckling when 7,/r, 13/24, three-wave changes to four- 
wave when 7,/7; 97/120, and as 7, and rz approach nearer to equality the buckled 


form progressively changes to more waves. When applied to a structure spread by 
many horizontal hoops of which the top one is rigid and oval, it is found that all 


other hoops, if each is designed to the criterion qr? = 3EI, will have the same 
rbsolute deviation from circularity as the rigid hoop. If any one hoop is designed 
so that qr? > 3EI, then the oval shape of the rigid hoop is magnified on all flexible 
hoops. 


1.1. Introduction 

THE stability of an elastic ring subjected to a uniform radial pressure was 
first investigated by Levy (1) about seventy years ago. This problem is 
described in the standard books (2, 3) on elastic stability with proofs similar 
to that of the original investigator. The method of proof, however, is based 
upon a geometrical theorem which is only valid when the pressure is con- 
stant, and cannot be extended easily to the case of a ring in which the 
radial pressure changes with the radial deflexion. Such a case occurs in 
a flexible heavy structure like the crinoline skirt whose mechanics have, 
apparently, never been analysed. Although there is a practical application 
for this analysis since certain plastic materials of limited elastic strength, 
but otherwise very satisfactory, have recently become available for crinoline 
hoops, the analysis has also helped to provide a better understanding of 
the stability and deformation of elastic rings. Moreover, the crinoline skirt 
is a structure which can be represented mathematically with some pre- 
cision, which is impossible in many aeronautical structures; it is a structure 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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in which large deformations are important since a flexible structure, and 
not a rigid one, is to be designed. Therefore, the use of relationships which 
assume small deformations must be watched very carefully throughout 
the work. 

This investigation was begun when it was realized that the account of 
the buckling of a ring given by Timoshenko (2) was not only approximate, 
but incorrect in principle. It seemed necessary, therefore, to start again 
from first principles. However, Fairthorne (4), some years ago, had 
realized that the standard method of stressing a ring was unsound and that 
the problem had not been solved for large deflexions. Moreover, he recog- 
nized that the published treatments were unsound in principle but thought 
that a more rigorous solution would be intractable. More recently Biezeno 
and Koch (5, 6) in Holland have given the buckling of rings much atten- 
tion, but their approach is not quite the same as the one described in this 
paper. By describing compressive systems as those characterized by zero 
shear force and zero bending moment in the undeformed state they have 
restricted the generality of their work. The practical experience that stable 
elastic rings can be formed out of material whose initial curvature differs, 
within limits, from that of the ring, suggests that this restriction should be 
avoided. Moreover, these investigators do not transform one of the basic 
equations by integration with respect to distance around the ring, a step 
which appears to simplify the problem, and add further complication by 
introducing harmonic functions at an early stage. 

In the first of the following sections the general principles of the analysis 
are discussed, pointing out certain fundamental differences between this 
and previous approaches. The next two sections develop the stability 
equation for a single ring with reference to a generalized loading function. 
The last main section shows how the hooped skirt provides an example 
of a structure in which the loading depends upon the displacement, to 
which the results of previous sections are applied. 


1.2. List of symbols 
1.2.1. Geometric 


s = distance along the circumference of the hoop, 


r = radial ordinate at an angular ordinate % of the deformed ring, 
7, = radius of the undeformed circular ring, 

r,, = radius of the undeformed circular pth ring, 

Cc radius of the undeformed coil of material from which a ring is 


formed, 
p = radius of curvature at a point 7, on the ring, 
radial displacement outward from the mean radius, 
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d = vertical spacing between the rings, 
§ = vertical inclination of the fabric in a hooped structure. 
..2.2, Mechanical 
E = Young’s Modulus, 
EI = flexural rigidity, 


M bending moment, 

M, = bending moment in a circular ring, 
S = shear force normal to a ring element, 
C compression force in a ring element, 


7 = distributed load per unit length normal to an element of the 
ring, 


t = distributed load per unit length tangential to an element of the 


ring, 
ly 7» = distributed normal load per unit length on an undeformed ring, 
T,, = tension in fabric above the pth hoop, 
W, = weight of the structure below the pth hoop. 


12.3. Ope rational 
k number of mode, 
A = d/du& the differential coefficient, 
F (A) = operating function for the pth ring. 
2. Definition of the analytical problem 
The ring or hoop in this problem is considered to be uniform in section 
round its circumference; it is also assumed to be thin so that the deflexions 


lue to bending moment swamp those due to shear and compression. 
Elastic bending obeying Hooke’s Law is assumed so that the flexural 
rigidity can be denoted by EJ. 

This problem, like others in elastic stability, is concerned with the effect 
f applied forces upon the shape of a structure, a ring in this case. The 
riterion of instability or of deformation should be equally apparent 
whether a kinematic or a geometric quantity is taken as the principal 
wiable. The quantities which would satisfy this role are the pair which 
rovide the basic relationship between the kinematic effect and the 
geometry. This basic relationship in problems concerning thin rings is a 
roportionality between changes in bending moment and of curvature in 
nelement of the ring. If the radius of curvature of the element is c before 
tis stressed and p afterwards, then the magnitude of the bending moments 
pplied at the ends of the element is EJ(1/p—1/c). The equation will be 


written introducing a negative sign before M, as Timoshenko (2) does; thus 


Bil: _ _M. (1) 


p C 
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Care must be taken, however, with the sign convention because the correct 
sign for c in relation to M affects the second-order terms of the equation to 
be developed. 

Although the curvature and the angular ordinate % could, in theory, 
determine the shape of the ring, curvature is rather inconvenient for 
graphing. As most previous investigators have appreciated, r, %, andw 
(radial displacement) coordinates are more convenient. Curvature can 
be changed to the radial displacement by use of the relationship in 


ee , } ] ] 
oe = Ce ee (2) 
dip? rm % Pp 


which assumes that the deformation of the ring is small compared with »,. 


differential geometry, 


Here, w is taken positively when directed outwards. The curvature can be 
eliminated between equations (1) and (2) and we obtain 
2 
aw ee Mro (3) 
dif * EI’ 
which equation is still essentially one for transforming variables. 8. Timo- 
shenko bases his treatment upon this equation, though, because of the 
nature of equation (2), he has not used it correctly. This particular 
equation must be irrelevant in the initial stages of the proof. His approach 
leads to a second-order differential equation in w for the stability con- 
dition for a ring under constant pressure. No clear lead is given of the 
steps to generalize to the case where the pressure varies with displacement. 
The logical method, however, is to examine the equilibrium of an element 
of a ring subjected to the components of force and to the turning moments. 
This method has been adopted by others, notably Biezeno and Koch (5), 
but the manipulation of the basic equations in this paper is different from 


that of these previous workers. 


3. The stability equation 


3.1. The equilibrium of an element of a ring 

An element ds is shown in Fig. 1 dissected from a section of a thin ring 
by cuts at A and B. The centre of curvature is shown to the right of the 
positive direction of measurement of s along the ring. The radius of curva- 
ture and the compression force C are both taken as positive quantities s0 
that a distributed normal force directed towards the centre is negative. 
Then, to conform with equation (1), an anti-clockwise moment applied at 
the positive end of the element will be in a positive sense. 

At the cut A the section of the ring towards the origin is subjected to 
a bending moment M, a shear force S, and a compression force C with 
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equal and opposite values on the other side of the cut. At the cut B these 
uantities are increased by dM, dS, and dC respectively. Ifa distributed 
wrmal force g and a distributed tangential force ¢ are introduced, then 
ll the major force and force derivatives acting upon the element ds are 
stated so that the necessary equations to be 
satisfied when an element is in static equili- 
rium can be written down. These equations 
vill be referred to axes fixed in the element 


nd. to first order, they are 


dS ( t) 
q ( 
Ag p 
0 centre of 
r lateral equilibrium: Lees. a 
curvature 
Ss j d¢ 5) 
(2 
p ds 


r longitudinal equilibrium 


dM 


S (6) 
ds 





for rotational equilibrium. 
In the manipulation of the above equa- Fic. 1. 
ms, the obvious first step is to use 
equation (6) to eliminate S from each of the other two equations. 
Moreover, equation (1) can be used to eliminate 1/p from equations (4) and 


When these eliminations are made an interesting equation is obtained 


t, (7) 
Cc EI} ds ds 


which, subject to the restriction that c is constant, can be integrated term 


from equation (5), f A dc 


vy term. It is by this step that the present treatment differs from that of 


iezeno and Koch (5). On integration with respect to s the following 
equation is obtained: 

M M- 
Cc 2H 1 


[he integral | tds can be significant if the applied forces are no longer 


C = ( tds+Q,. (8) 


rmal to the ring after deformation. When these forces act centrally 


then ¢ = g(dw/ds) and q,w would be a significant term. It is, therefore, 
mportant to pay attention in any particular problem to the way in which 
these forces act during the deformation of the ring. However, in the 
pplication discussed in this paper they remain substantially normal to 


he ring after deformation and the tangential force can be neglected as a 


second-order effect in many cases. Why this is so will be discussed later. 
j2,18 


Q 
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Incidentally, it should be noted that the average value of [ tds around 
the ring is zero, otherwise the ring would rotate. j 

Equation (8) also shows that the compression force C no longer remains 
constant in a deformed ring. The arbitrary constant C,, the value of ( 
when 1 is zero, does not necessarily represent the compression force in 
the undeformed ring. The interpretation of this constant will be found 
later. 


3.2. The equation for large deformations 
Progress can now be made from equation (4). The coefficient d?.M ds 
is substituted for dS/ds: the force C and the curvature 1/p are eliminated 


using equations (8) and (1) respectively. As a result of these operations 


(' M 
c 7 d, 


\ 


the following equation is obtained: 
d?M nM, we 
ds? - c " TH 
a2 0 3? 3 
PM | (: ar. VM VM ( 


pe 
ds? a BI 2cHI  2(EI)? c 


or q.- (9) 
Equation (9) is clearly capable of providing stability criteria. If we con 
sider the case of a ring whose radius is the same as that of the unstressed 
coil of material, loaded by a uniform radial pressure which is uninfluenced 
by displacement, C,/c is found to be equal to gq. For small deformations, 
that is WM <C,c, the equation 
om -(: a eye 0 (10) 
ds? ce OKI 
shows that the circular ring becomes unstable when C, is greater than 
3EL1I/c* or when q is greater than 3#J//c?. The proof would follow that 
given in standard textbooks (2, 8). 

So far, small deformations have not been assumed and, subject to the 
applied forces remaining normal to the ring after deformation and to the 
curvature being uniform before stressing, equation (9) is valid for all 
deformations which obey Hooke’s law. However, the analysis cannot 
proceed far without assuming some restriction to the amplitude of the 
deformation. Equation (2) assumes small amplitudes because, inter alia, 
l dw\? 
2, dis) 


is substituted for WV there is no justification for retaining the power terms 


a term has been neglected relative to w. Consequently, when w 


in WM without the above second-order term in w. The second-order terms 
in w and in M are both comparable and additive so that their effect must 
be watched in a ring assembly in which certain rings can have large 
deformations. 
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Although the power terms in M will not be retained for the main task, 


which the first-order term is zero is conceivable. Such a problem is that 
concerning the draping qualities of an elastic curtain gathered along the 
supporting edge. In that problem both 1/¢ and C, are zero in the draped 


portion. 


3.3. The stability equation for small displace ments 

If the stability of a ring of radius 7 is to be examined, small displace 
ments are adequate. The element ds can then be replaced by 7 dy. More 
ver, equation (3) can now be used to transform M into w. The first step, 

wever, is to show how the arbitrary constant of integration, Co, is related 


to the compression on the undeformed ring. If the distributed radial load 


the ring in the undeformed condition is q,, then the compression will 


€,%. There is also a bending moment M, because 


1] M, 


—— (11) 
C El 


Equation (11) is the boundary condition to be applied to 


d2M ( 372M? = or2 M3 C 
2 o\y 0 1 0 21-0 _ gg). 12 
r EI) 2E 1 2(E I)? il c | _ 


when d?M (djs? is zero all round the ring. Substituting VM, for M, it is 


found that M, 
C. - ‘ ) (13) 
/ 


» ‘ 
0 & 


showing that C, only corresponds to the compression force in the un 
deformed ring when M, is zero 

For small deviations of curvature there will be only small deviations of 
bending moment so that 1 can be replaced by M,+AM. When this 
substitution is made in equation (12) and the value of C, from equation 
13) is also used, then neglecting powers of A.M, it is found that 


( 


a ' | ro Io 2(¢ : 
—(AM)+(1 ET pe r2(dy—9), (14) 


(iw 


which shows that the stability of the ring is unchanged by using a material 
Whose unstressed curvature is different from that of the ring. It must be 
remembered, however, that this result assumes that the ring has an infinite 
stiffness in twist, for, in practice, a combined instability in bending and 
twist can occur if the radius c is much smaller than the radius %. Further 
more, a limited and irregular variation of ¢ from 7%, although deforming 


the hoops, would not prevent the formation of stable hoops. 
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Having established the fact that the stability is independent of the initia] 
radius of the material, the problem may be simplified for further investiga- 
tion by making ¢ equal to %. The quantity C, can now be described as the 
compression force in the undeformed ring and the small difference AM can 
be replaced by M in equation (14). Moreover, a stage has now been reached 
in the analysis when the radial displacement w can be introduced as the 
principal variable. By substituting for w from equation (3) it is found that 

d4w (2 dos dw v Yo r a f, owt ra Ls 
dif ea aT} ait | r sre) do 4 ET’ (15) 


at the same time making a note that g can be a function of w. 


4. The generalized stability of a ring under radial loads 

Equation (15) is a differential equation which can be used to study the 
generalized stability of a ring. Four boundary con- 
ditions are required to accompany it. The point of 
maximum displacement will also be the point of 
maximum curvature, and this point is selected as the 
origin for %. Thus when % = 0, 


dw/dé = 0 and dw/ds? = 0. 


These differential coefficients must also vanish again 





- when % = 27, but if only oval deformation of the 
Kia. 2 ring is being studied, they will also vanish when 
ws har (see Fig. 2). 


Now, q(w) can be expressed as a series in w by Taylor’s theorem 
because w is small compared with 75, 


' er 
q\w) Jot do + 5510 T--- (16) 


Because of earlier assumptions concerning small displacements, there is no 
need to retain more than the first two terms in this expansion so that the 
right-hand side of equation (15) can be written —réq,w/EI. The roots of 
equation (15) are clearly modified by this term, so that the stability of a 
ring is modified by the first derivative of the distributed radial load. This 
particular result is not clearly stated by Biezeno and Koch (5). Although 
the implications will not be considered in this paper, one must not forget 
that g, and qj could be functions of . 


The equation providing stability criteria has now become 


d*w I on tere) aw F | (1 ! ry To wT O\ 0. (17) 
did El dias | qo El | 


and since only periodic deformations are of interest it is easier to deal with 
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the equation by an operational calculus procedure than by classical 
methods. A periodic solution w = w,cos ky will therefore be considered, 
and the conditions under which such solutions are independent of the 
ymplitude w, will be found. 


In operator form, equation (17) can be written 


/ q,.r2 » a \q. xe 
}y4 | (2 lo 0\,21 44 (1 04 1o"0\ 0. (18) 
El ' dq | ET) 
where A represents d/di. For the type of solution selected A? k?, so 
that 8 , 3 
kA (2 | eels se (1 - "oy 10\Io"o 0 (19) 
El qo | El 


is the condition which will determine whether the particular loading on the 
ring causes the disappearance of resistance to a deformation wy, cos ky. 
The alternative condition for centrally directed forces has the factor 
l+r%qo/Go reduced to %%4/qp. This can be verified by replacing C, by 
(,—q)w in equation (12), retaining only first-order terms in w. 


Condition (19) can be rewritten 


, FL, %% 


q, = (e@—1)—4 20.1 
1 aT i] ( ) 
for normally acting loads and 
q (k2—1)? EI \ "% qo (20.2) 
' % k? 


for centrally acting loads. 

Equations (20.1) and (20.2) show that a positive derivative increases the 
stability of the system and a negative one reduces it. The further study 
of this result will be made with reference to the hooped skirt as a practical 
structure, the mechanics of which have been described qualitatively else- 


where (9). 


5. The stability of a hooped skirt 

5.1. The stress analysis 

The stability of a hooped skirt probably constitutes the best practical 
problem in which the distributed loads are affected by the displacement. 
The example of the structure for stress analysis is a heavy flexible skirt 
whose full shape is obtained by a system of equally spaced hoops each in 
a horizontal plane. The diameters of the hoops increase towards the hem 
ofthe skirt. Fig. 3 is an exaggerated section illustrating the stress distribu- 
tion at the (p 


l)th, pth, and (p+1)th hoops in a set of n hoops. At the 
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pth hoop the fabric of the skirt will be in tension under the weight of 


all the lower material. A symmetrical design is 
assumed so that this weight can be regarded as 
distributed uniformly around the hoop. The 





Ip 
th | . 4 fabric will be sloping outwards and downwards 
2 e | with the vertical component of fabric tension 
7, equal to the distributed weight. Thus the vertical 
p-t" 6, component of tension over unit length of hoop of 
radius 7 can be written W,/27r,, where W,, is the 
(7 weight of the structure hanging below the pth 
hoop. If @, is the inclination of the fabric to the 
vertical in the pth zone which is bounded by the 

Fic. 3. pth and (p+1)th hoops then, when 
co, > Gs» 


the fabric will bend over the pth hoop and exert a pressure on this hoop. 
If 7', is the tension above the pth hoop, then 
= sig W. 
T’, cos @, T,,_,cos6,_, E . 
/ / } I Dar 
p 
The horizontal components of the two tensions are unequal and provide 


the distributed load q(w,,) on the hoop 
q(w,,) T,, sin 6, T, , Sin 6,,_;. (21) 


) 


If the tension is replaced by the weight, then 


W 
‘ ‘ ‘ ) >») 
q(v,) (tan @,—tan 6, ‘= , (22.1) 
2rr,, 
W,, ' 9 
(,—8, 1)=— (for small angles) (22.2) 
<7Yr 
p 


If d is the vertical spacing between the hoops, and the radius vectors of the 
deformed hoops are 7,,+-w,,, 44 +Wpiaers in the common vertical plane 
at ordinate %, then 


rt+w,—r w 
tan@, p Pi a (23) 
d 
W, 
, ” we : , ', ee p 24) 
und Uy) = (%y+- Wy —% 41 — Wp 4a —"p-1 — Mp-r tM th js 7 (24 
2rr,, 


When all the w’s are zero the value of ¢ = q, for the undeformed ring is 
obtained. Therefore equation (24) can be rewritten 
W 2w,— w - 
q(w,,)—q (2 uv w ) P q P pti p= (25) 
p p p p-1 pti? 9. ] > eo. " 
2rr,, ¢ or —¥% ? 
p p p+1 p-l 


showing that the distributed load on the pth hoop depends upon the 





er 








defe 


equ 


bee 


in ¢ 


whi 
of { 
sup 
doe 
mal 
will 
first 
wit! 
asi 


the 

5 

} 
str 
of t 
the 
hoo 
the 


be ; 


ant 


use 


Thi 
dey 
pre 
anc 
tha 
ina 


bec 


rds 
ion 


cal 


» of 


the 
oth 
the 
the 


Op. 


ide 


) 


the 


une 








re ree 








STABILITY OF A COMPRESSED ELASTIC RING 231 


deformation of the (p—1l1)th, pth, and (p+1)th hoops. By assuming 
equally spaced hoops, an unnecessary complication in formulation has 
been avoided in the basic analysis; unevenly spaced hoops would result 
in different coefficients for each of the w’s. 

The tension in the fabric of the skirt has been assumed in the above 
stress analysis to act in the vertical plane at ordinate w. If it were restrained 
to this plane, then an allowance should be made for the tangential forces, 
which allowance is necessary in the case of a single hoop spreading a cone 
of fabric supported at its apex. The case is different, however, if the 
support is a rigid hoop. The tension is no longer centrally restrained; it 
does not act strictly in the plane % and the applied forces tend to act nor 
mally to the hoops. The tangential forces experienced by the lower hoop 
will be second order and, therefore, will be neglected since those of the 
first order act upon the support hoop. It has been observed in a skirt 
without a rigid supporting hoop that the uppermost flexible hoop buckles 
is if the compression force in the flat arc is excessive. The assumption that 


the forces act normally to the lower hoops thus appears justifiable. 


5.2. The case of two hoops 

Before developing a matrix for the solution of the deformations of a 
structure spread by a set of » hoops, it is worth while examining the case 
of two hoops. In equation (25), p now relates to the first hoop and, since 
the material will tend to hang vertically below this hoop, the (jp—1)th 
hoop can be looked upon as a fictitious hoop having the same shape as 
the first hoop; therefore w,,_, w,. The (p+1)th, or second hoop, will 
be assumed rigid so that it has no deformation. It follows then that 


W, uw 44 uw (26) 
2rrd 1 —" 


ind that q 4, (%—%). This is a simple expression for q’ which can be 


used in equation (19) to give 
3 
4111 (k?— 1)*k?—1 —,(7,—2) 1}. (27 
El 
This equation shows that the mode of deformation of the lower ring 
depends upon the relative size of the two rings and is purely a geometric 
property of the rings. The mode of deformation is characterized by k 
ind, for a particular combination of 7, and 7,, the appropriate value of & is 
that which gives q, 73, EI the lowest positive value, negative values being 
inadmissible. It is seen that, as 7, and 7 approach equality, 7,/(7,—7%) 
becomes large and / must be larger in order to keep the fraction positive. 


Che simplest mode of buckling corresponds to k equal to 2, which occurs 
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when the distributed load is independent of the displacement. It will be 
appreciated, however, that a value of 7,/(7,—7,) will occur for which a 


value of either 2 or 3 for k gives the same value of q,73/E7. Thus, when 


9 64 


: ; . ° (28) 
3—1,(7,—12) 8—7,(7,—12) 


a condition is obtained which gives the numerical value of the ratio of , 
to 7, at which value either two-wave or three-wave deformation of the 
ring is possible. It is a transition ratio below which only two-wave de- 
formation occurs and just above which three-wave deformation occurs. 
The value of the ratio is 75/7, }3. Three-wave deformation has been veri- 
fied in a qualitative experiment. 

As the size of the rigid hoop approaches nearer to that of the flexible 
hoop the number of waves in the buckled flexible hoop increases one by 
one. Transition points occur at values of the ratio 7/7, given by the 
general condition 


(k2—1)2 ((k-+1)2—1)2 a 
k?@—1—1,(r,—1%)-  (kK+-1)?—1—2,(4.—2) _ 
from which it follows that 
)]-2_1 9) 
: 2h" +2k—1 (30) 
r, (k—1)k(k-+ 1)(k+2) 


Moreover, when the value of 7/(7,—7) obtained from equation (29) is 
substituted into equation (27), the value of ¢, r?/#I is given by 2k?+2k—1 
at the transition from the kth to the (k-+-1)th mode. 

Equation (30) has been used to compute the following sequence of 
transitions. 





Transition rer; rise! 
Two to three waves | 13/24 0°542 II 
Three to four waves 97/120 0808 23 
Four to five waves 107/120 0892 39 
Five to six waves 779/840 0'928 59 
Six to seven waves 1§97/105e 0°953 83 


These transitions should be verifiable by experiment, at least the first two 
or three of them, but it may be found that corresponding values of 1/%, 
and q,7j/EI quoted above are respectively higher and lower in practice 
because of the neglect of the tangential forces during deformation. In 
Fig. 4 these transitions show up as discontinuities in the relation between 
r/r, and q,7?/EI. It should be noted that the above results do not depend 
upon the stoutness of the flexible hoop. 
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0 2 - 6 8 1:0 


5.3. The case of n hoops 


Equation (25) giving the load q(w,,) on the pth hoop is 


2w w 


dy» q\w ,) e p p-1 p (25) 
/ I / Dr r r 
“2 p-l p+1 
However, at the lowest hoop this equation reduces to 
w w 
. 1 2 ¢ 

q q\W’}) 44 . —* (31) 

Y > 


These values of g,,—q(w,) can be substituted into equation (15) which is 


uso rewritten in operator form. The stability equation for the pth hoop 


becomes 


i" \ Ply ) ’ 29 
| F(A) +2 |w, Wyss tWy-r (32) 
where 
(2r r r \E] q r3 q ioe 
F (Xr) ee p-l pire }y4 119.1 4p Pp 2 | J]. 1P Pp ‘ 
: wn \ VEL, EI, | 


ind F(A) is the operator function. 

The stability of a heavy flexible structure, such as a skirt extended by a 
system of horizontal elastic hoops, can be represented by a set of (n—1) 
equations, all similar to equation (32), except for the one which relates to the 


lowest hoop. The matrix for the coefficients of these equations is as follows: 


[1 —[F (A)+2] l 0 0'0 0 0 
v | F,-1(A)+ 2] 0'0 0 0 
v V0 l F_(A)+2]1:0 0 0 
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If the shape of the supporting rigid hoop, the nth in this case, is defined 
relative to the mean circle, the deformations of the other hoops are 
calculable. In the case of a simple oval shape whose deviation from the 
circle is represented by w,, w, cos 2%, the deformation w,, will be of the 
same form and the values of the operating functions F(A) can be evaluated 


p 


by putting i? 4 in each one: F(A) when k 2 will be denoted by 
F(2j). If q,73/ EL, = 3 for every hoop except the nth, which is rigid, ‘ 
case of special interest is obtained in which F,(2j) F(2j)) = 0. The 
matrix then becomes 


l 2 | 0 0:0 0 0 0 
0 l 2 l 0:0 0 0 0 
0 0 l 2 1:6 0 0 0 
0 0 0 0 Oo} 3 2 l 0 
0 0 0 0 0:0 | 2 l 
0 0 0 0 0'0 0 l l 








[t is seen that all major determinants arising from this matrix are unity 
so that each hoop would have the deformation of the support. Since the 
deformation is absolute and not relative, it may appear insignificant on 
the lowest of a series of progressively larger hoops. 

It is not proposed to enter upon the rather complicated calculations of 
particular assemblies. Instead, a rather qualitative discussion of the trends 
in the shape of a skirt in which the hoops are at or near the critical stiffness 
is made with reference to the preceding matrix. 

If F(2j) is negative for any hoop, then the deformation of all hoops 
must be greater than that of the supporting stiff hoop. However, the 
amplitudes remain finite so that stable yet non-circular hoops are obtained 
when certain of these operator functions are negative. In fact, very in- 
teresting shapes are possible in which the effect of higher-order operator 
functions with negative values should be considered. This is particularly 
necessary when the shape of the rigid hoop contains a significant compo 
nent of one of these high orders. Because of this, a lead is given to a method 
of analysis of the folds of an ordinary draped skirt; the transition from the 
crinoline to the ordinary skirt could be illustrated by a suitable selection 
of numerical examples. 

A theory based upon small deformations can show that the lowest hoops 
of the skirt can have large oval deformations. It is, therefore, important to 
look back at the basic equations for the second-order terms which, although 
mentioned in the earlier discussion, have been neglected. Qualitatively, 


the trend of error of the small deformation theory is in under-estimating 
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‘he maximum inward displacement. The terms that have been neglected 


| have a negative sign when w is positive. Therefore, if a substitution 
w = uw, cos 2s is made into the equations with the second-order terms 
etained, a cos*2ys term with a coefficient involving the ratio w,/r, arises. 


) 


1 


is a harmonic whose effect does not appear to be as important as that 
fthe neglected tangential forces. A precise account of second-order terms 
wld have little practical value so long as the possible variation of the 
irvature l/c around the hoop is neglected. Although such a variation 
s been neglected to simplify the algebra, it can be very important 
shen the stiffness of the hoop is nearly critical. Near the critical stiff- 
ess a cos 2 modulation of the curvature l/c of amplitude equal to 1/c 


in cause a multi-hooped structure to cave in prematurely on one side. 


6. Concluding discussion 
\Ithough this analysis was started to assist in the design of a flexible 
ped skirt, it has also helped to clarify the theory of the elastic stability of 
ngs. The novel analytical results obtained could be verified quantitatively 
fairly simple experiments. The analysis has not been too encumbered by 
pproximations and the progress up to equation (9) has not assumed small 
leformations. That equation could be used to examine by a series solution 
the large deformations of a circular ring and would provide an alternative 
ethod of obtaining certain results quoted by Fairthorne (4). 
It 


there are localized loads. Such loads produce discontinuities in the shear 


s probable that the method of this paper becomes complicated if 


force. Some well-known results for struts and ares can be deduced from 
equation (9) by applying the appropriate boundary conditions. It should 
e noted, however, that the integral approach is better suited to solve strut 
roblems than the differential approach which, as has been shown, solves 
roblems of rings more directly 

Shear and compression forces have been fully accounted for, but 
leformations due to them have been neglected by considering thin rings. 
\ generalized distributed load has been considered on somewhat simpler 
ines than those adopted by Biezeno and Koch, although only the case of 

non-rotating ring with a normal loading which is independent of the 
rientation has been worked out in detail. For a uniform distributed load 
sensitive to displacement, it has been shown that the stability also depends 
ipon the first derivative with respect to the radial distance, the ring 
becoming more stable if the loads tend to increase with an outward 
deformation. Owing to this modification to the stability, the modes of 
nstability of a flexible ring supported below a smaller rigid ring differ 


iccording to the relative sizes of the rings. 
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Finally, it has been valuable to have an application in which th 
qualitative shape and deformations were of primary importance and the 
stresses of secondary importance. This point led the author to adopt 
the method used in this paper and to disregard the accepted methods 
for stress calculations. 
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UNIQUENESS THEOREMS OF TWO-DIMENSIONAL 
ELASTICITY THEORY 


By R. TIFFEN (Birkbeck College, London) 
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SUMMARY 


his paper is concerned with elastic material occupying a multiply-connected 
nand in a state of plane strain or generalized plane stress. For a finite material 
s shown that a solution giving specified stresses over all boundaries is effectively 
nique, i.e. the stresses are uniquely determined throughout the material and the 


wements are unique except for rigid body displacements, and that a solution 


ng specified displacements over all boundaries is unique. The theorems are stated 


hat they retain their significance if the boundary loading is discontinuous, and 
proofs make no appeal to special forms of the complex potential functions. The 
therto neglected problem of material for which the exterior contour is partly at 
ty and partly in the finite part of the plane is considered, and a solution giving 

d stresses over all boundaries in the finite part of the plane and specified 
resses at infinity of orders higher than o(r~1), where r is the distance from the 
gin, is proved to be effectively unique. The simple requirement of vanishing 
tresses at infinity is shown to be, in general, insufficient to define a unique elastic 
roblem, whilst on the other hand stresses at infinity which are o(r~!) may not be 


bitrarily specified. A solution giving specified displacements over all finite bounda- 
and prescribed finite and infinite displaccments at infinity is shown to be unique, 


t an arbitrary assignment of infinitesimal displacements at infinity is not per- 
ssible. The special cases of material extending to infinity in all directions and of 
terial occupying a half-plane are considered in detail. 


1. Introduction 
[HE uniqueness theorem given by Kirchhoff (1) in 1859 applies only to finite 
naterial and refers to three-dimensional problems. The investigations of 
Muschelisvili (2), (3) are concerned with two-dimensional problems and are 
onfined to the following cases 
i) Material bounded by n+-1 closed contours C,, C),...,C,.,, where C,,., 
encloses ¢ a5 4 , 
ii) The limiting case of the above as C,, ,, tends to infinity in all directions. 
For infinite material Muschelisvili’s treatment requires the material to ex- 
tend to infinity in all directions and not to be surrounded by a contour partly 
tinfinity and partly in the finite part of the plane, as in the case of the elastic 
half-plane or the elastic strip. There is evidence (4), (5) suggesting that 
Muschelisvili has pursued his investigations further but, if so, the relevant 
papers are not readily available. Below we give proofs of uniqueness which 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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apply to more general limiting forms of C,, ,,. We also discuss the important 
questions of what must be postulated if part of the contour is at infinity and 
of how much may be specified without rendering the problem insoluble 
Neglect of these two questions has led to vagueness with regard to the postu. 
lation of elastic problems when infinite material is involved. Also, it is not 
clear, in general, to which finite problems such cases approximate. The 
mixture of real and complex variable theory used by Muschelisvili is replaced 
by proofs entirely in terms of complex variables, with a consequent gain ii 
simplicity. Information with regard to the uniformity or non-uniformity 
of the complex potentials is obtained from the equations defining the 
physical problem without appealing to special forms of the potential fune- 


tions. In the only case in which we appeal to the logarithmic function a 








simple proof is given that it is the only function possible. The statements 
of the theorems are framed so that they retain their significance in those 
cases in which the potential functions have singularities on the boundary 
contours. These cases have hitherto been excluded from uniqueness 
theorems although they are of the greatest importance. 


2. Fundamental equations 

The notation used is that used by Stevenson (6), (7) who has shown that 
if complex coordinates z = x«+-iy, Z = x—iy are used, the solution for the 
complex displacement D and the stress combinations ©, ® in plane strain 
are given in terms of two functions Q(z), w(z) of a single complex variable 
(the complex potentials) by 


8uD Su(utir) KQ(z)—2Q!(Z)—o' (2), (] 
20 = 2(axr+yy) = Q'(z)+0'(2), (2 
2 2(ca—yy— xy) = 2Q"(z)+w"(z), (3 


assuming no body force or couple. Here , is the rigidity modulus and 
K 3—4n, where 7 is Poisson’s ratio, and accents denote derivatives. 

The same formulae hold for the mean displacements and the mean stresses 
“x, yy, xy of generalized plane stress, provided that 7 is replaced by the 
modified Poisson ratio o, where (l—o)(1-4 n) 1, so that x 53—4o. 
Throughout this paper the analysis is made in terms of generalized plane 
stress, although it is clear that this may equally well be interpreted in terms 
of plane strain. 

Continuity of the stresses and displacements requires the functions (2), 
w(z) and their first and second derivatives to be continuous at all interior 
points of the material, i.e. in the open region obtained by excluding the 
points on the boundary contours. Assuming this continuity, it can be shown 
from equations (2), (3) that the resultant foree X--iY and the moment V 
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hout the origin of the stresses acting across an are AB in the material 


proceeding along the are from A to B the forces act on the material on the 


eft of AB) are given by (Stevenson, 6) 


; =\)B i 


[47 A (z, 2 A(z. Z).— A(z, 2), 1, (4) 
[AG,2)]% = “IAG. 2p ; 
vhere aie. Z €)(z) 20)’ (z) 1 w'(z) (5) 
B 
nd t.\ re | = dA re| 2Z02'(z) 1 2q'(z) w(z) |". (6) 
1 
From (2) and (3) we obtain the useful stress combination 
2(0—®) = 4(Fy+ixy) = Q'(z)4+0'(2)4+20"(z)+w"(z). (7) 
The function A(z,Z) introduced above is a non-analytic complex potential 


f considerable importance in two-dimensional elasticity theory. 

We now give four lemmas useful in the main theorems. 

Lemma I. The function A(z,2) is constant along any stress-free arc in the 
terial. This result has been given, in a different notation, by Stevenson (7). 


The three lemmas given below follow from results obtained by Steven- 


son (6) 
Lemma II. The function 
8uD = KQ(z)—20'(2)—@'(Z) (8) 
ro along any arc in the material at all points of which the displace ments 


Lt} MMA II] F r non dislow ational proble ms 


Cy, Q'(z) = Cy,[«Q(z)—@'(2)| = 9. (9) 


Cy, denotes the cyclic va round any closed contour 6 in the open 


' 1 ° 
70n dé fined Dy the mate rial. 


Proof. Uniformity of the stresses requires 


Cy,[Q’(z)+0'(2)| = 0 (10) 

d Cy,f2Q"(z)+-w"(z)] = 0. (11) 
Uniformity of the displacements requires 

Cy;,{ Q(z) —20' (2) —a' (2) | = 9. (12) 

Differentiating (12) with respect to z and using commutativity of the 


perators Cy ind é/éz (see (8), p. 572) 
Cy,| eQ’(z)—Q'(2)|] = 0. (13) 
Equations (10) and (13) give 
Cyg 22'(z) = 0 (14) 


ind (12) then gives Cy;,| «Q(z)—o'(Z)| = 0. (15) 
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LemMMA IV. /f all boundaries of the material are stress-free and there are no 
external forces in the interior of the material,+ then 
Cy Q(z) = Cygw'(z) = 0 (16 
for all closed contours © in the open region defined by the material. 
Proof. For equilibrium we must have, from (4), (5), (14), that 
Cy,{ Q(2) + w'(z)| 0: (17 
equations (15), (17) then give (16). 
ANALYTIC THEOREM. The two-dimensional Gauss gradient theorem may 
be stated , . 
| én ds grad ¢d dS, (18) 
Ss 


( 
where d = 4(2,y). 6b/ Cx, 6b/Cy are continuous inside and on the contour ( 
enclosing a region S, s denotes distance along C and n is a unit outward 
normal vector to C. In terms of complex variables n ds is replaced by —i d: 

» c = yy 
and grad ¢ by 2 —¢(z,2). Then 
OZ 
° ; Pg Bes: 
d(z.2) dz 2 | - (2,2) dS (19 
Oz 


« . 


C S 


which implies d(z, Zz) dz 27 © dz. z) dS. (20) 
° : Ce 

Equations (19), (20) are sometimes referred to as the complex Stokes’s 

theorem (9). Clearly equations (18), (19), (20) may be applied when ¢ is 

a complex function of # and y. 

General theory requires that in the interior of elastic material the stresses 
and displacements must be continuous. This has been ensured by continuity 
of Q, Q’, Q", w', w" in the open region defined by the material. Previous 
writers have assumed, in addition, that Q(z), Q’(z), w’(z) are continuous in 
the closed region obtained by addition of the boundary points of the material 
to the open region. When the latter conditions are satisfied Lemmas I, I] 
[1Il, 1V apply to the boundary contours of the material as well as to interior 


arcs. 


3. Uniqueness theorem for finite elastic material with specified 
stresses over all boundaries 
Let elastic material occupy the region between contours C,, Cy,..., C, and a 
finite contour c, 4 surrounding C1; a C and let the stresses be specified 
over all boundaries. If Q(z), w(z) are complex potentials satisfying the given 


boundary conditions, then the only additional potentials Q4(z), w(z) which do 


+ This is assumed throughout this work unless otherwise stated. 
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iolate the boundary conditions and are such that Q,(z), Q6(z), wo(z) are 


ntinuous in the closed re gion de fined by the mate rial are 

(2,(z) wz+-a, w(z) = Bz+y, (21) 
here cis a real constant and x, B, y are complex constants. 
Remarks on the above theorem. The additional potentials give a rigid body 


splacement and hence give rise to no stresses. Note that we do not neces- 


arily suppose that Q(z), Q’(z), w’(z) are continuous in the closed region, but 


ve asserted that, whether these functions are continuous or not, the 
tentials Q(z), w(z) are unique to the extent that no additional potentials 


satisfying the extra conditions of continuity may be added, other than those 


siven in (21). The functions Q(z), Q’(z), w’(z) may have singularities, inside 
ndon the boundary, corresponding to isolated forces in the interior of the 
iterial, to isolated boundary forces, or to finite discontinuities in the 
undary loading. Such singularities, however, are supposed to be neces- 
sarily common to any two particular solutions Q,(z), w,(z) and Q,(z), w.(z) 
the same problem, and the difference functions Q,(z), wo(z) are free from 
suchsingularities. Further, singularities specifying the degree of dislocation 
: multiply-connected material would be common to any two particular 
lutions with specified dislocations, and the difference solution would be 
free from dislocations. Again, we may include body forces, for two particular 
utions with the same body forces give a difference solution with zero body 
oun 
With the above remarks as to what is specified in a given elastic problem, 
the uniqueness theorem resulting from (21) is very general. 
Proof of the theorem. The entire contour C,, C,,..., C,,.,, together with the 
uts necessary to make the region simply connected (traversed once in each 
lirection) will be referred to as ( 


Consider the integral 


| | [024(z)Ag(z, Z) dz — Q4(Z)Ag(Z, z) dz], (22) 

where the functions involved refer to the additional potent ials defined above 

ind are assumed to satisfy the extra conditions of continuity which make 
> 


he integrand of (22) continuous in the closed region. These potentiats leave 


the boundaries stress-free, so that, by Lemma I, 


AA2. 2) K, on C, (r SS en n+1) (23) 
vhere K,, K,,..., K,,,, are complex constants. Thus 
[= Ki 4, Cye, Q(z) — Ky as C¥e, ., Qo(Z) 
S[K Cyp, O,(z) — K, Cye, OQ9(2)]A4 


r=] 
contributions from the integrals along the cuts. (24) 
R 
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In view of the uniformity of Q,(z), the integrals round C,, C,,.... C141 are 
zero, whilst uniformity of the integrand in the closed region gives zero con. 
tributions from the cuts (see Lemmas I, III, IV). Hence 

I= 0. (25 
Transforming (22) to a surface integral by (19), (20) 


. 


2i | 


. 


’ 


%o(2) = Ag(2, 2) (2) “ 


(2) — Ag(Z,2)| dS, 
Cz 


— 


¢ 





where S is the region occupied by the material. Thus 
I 27 | [Q4(z) + Q) z) |? das. (26 
S 

The transformation is justified by continuity in the closed region of the 
functions Q5(z), A g(z,Z) relating to the additional potentials. Since the 
integrand in (26) is real and positive and / 0, we have 

0% (2) +0;(2) = 0 (27) 
and Q)(z) ic, where c is a real constant. Hence 

Q),(z) icz+«a, (28) 
where « is a complex constant, and on C, 


w,(z) = K,—4a; (29) 
since w,(z) is analytic and uniform in the open region and is continuous in 
the closed region, we may apply Cauchy’s integral theorem, so that, for any 
point z in the open region, 


F ] wy(C) 5- 
Wo(2) = : | Oe 
2m J C—z 
2 
l ’ K Rae — (K—é,. : : 
- | | nsa ge * | tr ac nk 
2mi| . C—z a1 J) C-z '| 
Gicas r=1 
Hence w,(z) = K,,.,—a% = B (say) (30) 


and, by continuity, 8 is the value taken round all contours. Thus 
wo(z) = Be+y. (31) 


4. Extension of the theorem to the case of infinite material 

Infinite material must be regarded as the limiting form of finite material. 
The theorem as proved clearly holds if all or part of C,,,, tends to infinity, 
provided that the stresses at infinity are specified exactly. Such information 
is not, in general, available and we shall see that an arbitrary assignment 
of the values of the infinitesimal stresses at infinity is not permissible. For 
infinite material we distinguish two general cases according to the orders 
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C.,., are | of magnitude of the total lengths of those parts of C,,,, which are made to 
zero con- | tend to infinity. 

Uniqueness of solution when C,,,, has sections at infinity of total length O(z) 
(25 Let us restrict the additional potentials Q,(z), w9(z) so that they give zero 
stresses over the finite part of the boundary and in addition we have 


Q(z) o(z-1). wh(z) o(z-), (32) 


i.e. lim 2Q)(z) = 90, lim zw(z) = 0 (33) 
on those parts of C,,,, which are to tend to infinity. 
Referring to these additional potentials, A,(z,Z) is bounded everywhere 





n of the so that the integrand in (22) is o(z~!) on those parts of C,,,, at infinity. If 
nce the | these parts of the contour can be chosen so that their total length is O(z) at 
| most, then the limiting value of the integral J is zero. The remainder of the 
- | proof is the same as above except that the possible additional potentials 
seeds Q(z) =a,  w(z) = Bet+y (34) 
29) | in view of the orders specified at infinity. Also there will be an evanescent 
term in the equations leading to (30). We are thus led to the following: 

: Statement of the theorem when C,, ,, has sections at infinity of total length O(z) 
” Let Q(z), w(z) be complex potentials which (i) satisfy boundary conditions 
lous in of stress over all boundaries in the finite part of the plane, (ii) have specified 
or any terms in Q'(z), w"(z) of degrees equal to and higher than z~—' on those parts of 
C,., which are to tend to infinity. The given stresses and the given terms at 
infinity must be such that the limiting form of the material experiences zero 
resultant force. The only possible additional potentials Q)(z), wo(z), not violating 
these conditions, for which Qo(z), Q5(z), wo(z) are continuous in the closed region 
and of the orders given in (32), are the elementary potentials given in (34). 

Remarks on the above theorem 

(30 From (2), (3) itis evident that we are demanding knowledge of the stresses 
at infinity other than those which are o(z~'). Application of (4) to the 
contour which is to tend to infinity shows that the terms of the orders defined 
in (32), which give stresses at infinity of orders o(z~!) can make only an 


infinitesimal contribution to the resultant force. However, we have to con- 


erial. sider only terms in Q’(z), w"(z) which are o(z~-?) before we can be sure of 
nity obtaining infinitesimal contributions to the moment NV from the limiting 
ation part of the contour at infinity. There is, however, no need to specify all the 
ment terms at infinity involved in the moment NV, as those of order lower than 

For that given in (32) automatically take their unique values when the higher- 








rders order terms take their correct values. Moreover, it is not permissible to 
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specify arbitrarily those terms at infinity which give infinitesimal contribu- 
tions to both the resultant force and the moment, although these terms can- 
not affect the equilibrium of the material. This explains why the limiting 
form of a boundary-value problem for finite material when some of the 
boundaries tend to infinity is often unsatisfactory. 


The importance of the above uniqueness theorem may be illustrated by 
{ : y 


a: call 
the following special case. Let z z(C) (S Be , where € €+-in, and 
21 
consider the region 7 > 0 bounded by the parabola 7 = 0. The potentials 
(2, (C) b w,(C) 1C3/6+ 5ia?Z/2 
and (2,(C) tC, w(C) ¢3/2—a?C/2 


both leave the boundary » = 0 stress-free and, at infinity, 

Q,(), Q4(¢) = O() = O(2!) 

w,(C). w»(C) O(C) O(2°). 
Thus both pairs of potentials give stresses at infinity which are O(¢-'), ice. 
O(z~*), and may be superimposed on a solution to a problem of specified 
stresses over a parabolic boundary balanced by vanishing stresses at infinity 
(see e.g. (10)). This superposition would, however, make a profound differ- 
ence to the stress distribution in the finite part of the plane, for the hoop 


stresses on the boundary 7 = 0 are 
i 9 = 9 
e, a, oe = 
- +a? sia +a? 


We now consider in detail two particular limiting forms of C,,,,. Our aim 
is to justify the demands which we are making with regard to the stresses 
specified at infinity. 


5. (1) The boundary value problem when (,,, tends to infinity in 
all directions 
In this case elastic material occupies the entire region outside any circle Q 
with centre at the origin and of sufficiently large radius. In this region 2'(z) 
is analytic and uniform, by Lemma IIT. and may be expanded in a Laurent 
series about the origin. Thus 


Q'(z) = A, tid, +P P24 6 (2-1) (35) 


and hence Q(z) (A,+7A,)z+(p,+ip,)logz+ E+0(1), (36) 


where £ is a complex constant and A,, Ay, p;, py are real constants. Higher 
powers of z have been omitted as they would give infinite stresses at infinity. 
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By Lemma IIT, «Q(z)— 6’ (2) is uniform outside Q,. Hence, if Q is any contour 


enclosing’ GV, 


Cyy @ (2) 2rik(p,+tps) (37) 
Ol Cyp w*(z) 27rik(p,—tpo) (38) 
so that Cyg| w (z)+-«(py ip,)log z| 0. (39) 
Hence w’(z)+-«(p,—ip,)logz is analytic and uniform outside Y and may be 


expanded in a Laurent series about the origin. Hence 
w'(z)+«(p,—ip,)logz = (B,+7B,)z+ F+o0(1), (40) 
where F is a complex constant and B,, B, are real constants. Thus 
w (Zz) k(p,—ip,)log z+-(B,+1B,)z- F-+o0(1), (41) 
where we have again omitted powers which would give infinite stresses. 
From (1), (2) the finite stresses at infinity are given by 
Lory By, brx 2A,—B,, 4iy Bb. (42) 
Statement of any problem would, therefore, imply immediate knowledge 
gives a rigid body displacement 


9S 


of the unique values of A,, B,, B,, whilst A 
and may be put equal to zero. Application of (4) to the limiting form of Q 
gives for the resultant force X+-iY due to the remaining terms, 


» 


p, = vX, DP. = vY where v —-, (43) 
m(«x+1) 

Thus when terms of higher degrees have been specified, the terms in Q’(z) 
w"(z) which are O(z~!) are uniquely determined by the necessity for zero 
resultant force on the material, i.e. p,, p, are known. The remaining terms 
are of the orders given in (32) and are, therefore, unique save for terms of 
the forms given in (34). Thus in this case the constants A,, Ay, B,, B, may 
be chosen arbitrarily provided that they, together with the applied stresses 
over the finite boundaries, give rise to a finite resultant, but no other terms 
may be specified arbitrarily save for terms giving a rigid body displacement. 
In this case, which was considered by Muschelisvili (2), there is no possi- 
bility of vanishing terms in 02’(z), w”(z) of higher degrees than z-!. This 
results from the Laurent expansion for Q’(z), given by (35), the validity of 
which is due to the fact that elastic material is assumed to extend to infinity 
in all directions. If the latter is not the case, evanescent stresses of higher 


degree than z~! can, as we have seen, occur. 


6. (2) The boundary value problem when C, ,, tends tothe boundary 
of a half-plane 
Let the contours C;, C,,..., C,, be in the finite part of the plane and the 
contour C!,,, consist of the real axis between the points x - R and the 
semicircle, centre the origin, radius R, above the real axis. We are concerned 
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with the limiting form of C,,,, as R tends to infinity. Let the stresses be 
given over C,, C,....,C,, and let external stresses be applied, in addition, over 
y = 0so that 


YYy-0 = Y9 =flx), zy = g(x) (44) 
where lim af (a) L, < ©, lim x h(a) L. < @. (45) 


Also let Q’(z), w"(z) be O(z-") at infinity in y > 0 so that, for large |z|, we 
may write 
Q’(2) = PAT Ps 1. 6 (2-1), w"(z) = 42 1 9 (2-1), (46) 


~ ~ 


We then assert that the stresses are unique. Terms in (’(z), w”(z) of higher 
degree than assumed in (46) may be introduced without seriously altering 
the proof. It is worth noting that, if the stresses over the real axis and the 
contours C;, C3,..., C,, have a non-zero force resultant, the degrees assumed 
in (46) for Q’(z), w”(z) at infinity are the lowest possible. If the above stresses 
give zero resultant force but have a non-zero moment about the origin, then 
the lowest possible degrees for 2’(z), w”(z) at infinity are 


Q)'(z) O(z *), w"(z) = U2 a) (47) 


Proof of theorem. We require to show that p,, 5, 7,, Jo are uniquely deter- 
yo 1 1 P2> W1> I2 {uel 
mined. The general theorem then shows that all the terms which are smaller 
at infinity are also unique. From (46), (7) we have, for large |z|, 
he oes +1, . P,—iD. ; 
4(yy+ixy) Pit Pe, Pr a" d 2__ (p,+ips) 


~ = ~ 


as (1+%92) +-o(z-). 


Hence, for y = 0 and for large |z|, 
4x] f (x) +-ib(x)] = (pyt+1)+%(q2—P2)+0 (1). 

From (45) we have 

Pit Gs Ly, q2—P2 = Ly. (48) 
Let the resultant of the stressesover thereal axis and the contours(, ,Q),...,C, 
be —(X+iY). For equilibrium of the material the force resultant of the 
stresses over the semicircle as R tends to infinity must tend to X +7Y. Hence, 
from (4), 








X—i¥ = 7[(p,—1)—i(pa+ 42) +-0 (0). (49) 
From (48), (49), 
2X JL, [, 2X 2Y LL, 2Y Lf, 
Pi er ee 11 eee, Pe — 9 ie a 


Hence the constants p,, Po. 7,, Jo are uniquely determined in terms of speci- 
1» P2> 11: C2 juen) 
fied stresses. In many problems we assume, for the purpose of application 
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{the Fourier integral formula, that f(x), 4(x) satisfy the conditions 


( f(x)| dx 0, ( d(x) dx < a 


nd L, L, = @. 
this case P,, Po, 91, Jo are given by 
Py Vy 2X /z, Pe qo 2Y /x. 
\n almost identical theorem can be proved for elastic material which can 


mapped conformally on to a half-plane. 


7, Uniqueness of solution when C,,,, has sections at infinity of 

finite total length 

In this case we are concerned with material for which we specify exactly 
the stresses over all but a finite section of C,, ,, at infinity. A typical example 
sthe infinite strip with specified stresses over the straight boundaries. Let 
isrestrict the additional potentials Q,(z), w»(z) so that they give zero stresses 
ver the finite parts of the boundaries and, in addition, 

(2) (z) ia), wo(2) O(z-1) 

for large |z|. In this case the integrand of (22) is O(z—logz) at infinity so 
that J = 0. Thus there is no need for a special consideration of terms in 


('(z), w”(z) which are O(z-). 


8. Uniqueness theorem for elastic material with specified dis- 
placements over all boundaries 
Let C,, Cg,..., Cy 4 be finite closed contours with C 
ind let complex potentials Q(z), w(z) give specified displacements over all con- 


2 Y Y Y 
n+ enclosing C,, Cy,..., Cy; 
tours. The only additional potentials Q,(z), wo(z) not violating the boundary 

miitions and with Q,(z), Q6(z), wo(z) continuous in the closed region defined 
hy the material are 

Q),(z) ¥ Wo(Z) Kaz - B, (50) 
vhich gue zero di splace ments eve rywhe re. 


Proof. The contour C is as defined above. Let 


J = 8p | [0Q5(z)Dg(Z, z) dz —Q5(Z) Do(z, 2) dz], (51) 
o 
where D,(z, Z), Q5(z) refer to the additional potentials; by Lemma IT, D,(z, 2) 
vanishes on C, C,,..., C,,,,. Also Q5(z) and D,(z,Z) are continuous and 
uniform in the closed region defined by the material so that the integrals 
long the cuts are zero. Thus 


J = @. (52) 
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Transforming to a surface integral by using (19), (20), we have 
sui | loti) ~ = 0 (z)£ 5 
J 162 (29(z) = Dy(2, z) + Q6(2) — Dy(z, 2)| dS 

C 


» ww 


2% | [{0Q6(z)}? + {Q4(Z) PP — 2«Q5(z)Q4(Z)] dS, (53 
where S is the region occupied by the material. The transformation js 
justified by continuity and uniformity of Q5(z) and D,(z,2Z) in the closed 
region S. Let ; , } 

' 0),,(z) T1+1T. (54 


where 7, and 7, are real functions of x and y. Equation (53) then becomes 


47 [ [ 7Z(«—-1)+-73(«+1)] dS. (55 


The elastic constant « is 3—4y in plane strain and 3—4o in generalized 
plane stress. Since 7 A/2(A+p) and (1+ 7)(1—o) 1 it follows that 
K 1. Thus, if we are to include the case x 1 (which was not considered 
by Muschelisvili), we conclude from (55) that 


To 0, 7} c(k), (56) 
(9, if K ] - 
where e(x) =: Su (01 
\c, a real constant, if « 3 
Hence (),(z) c(k)z+a, (58) 
where « is a complex constant. On the contours C;, C)...., C,, ., the function 


D,(z,Z) = 0, so that 

| c(«)z+a«]—e(«)z—@o(Z) = 9. 
i.e. wo(Z) = K[e(«)Z+&]—e(k)z, 
so that, for 1, we have, 

w(z) = Ka (59) 
on Cj, Cy,..., C,+,. By Cauchy’s integral theorem (59) holds throughout the 
material and 8uDp «[e(K)24 x] —[e(«)z+ Ka] 0. 
Thus the additional potentials give zero displacements. 


9. Extension of the theorem to the case of infinite material with 

sections of C,,,, at infinity of total length O(z) 

If C,,,, tends partly or wholly to infinity the above proof of uniqueness 
still holds but requires complete knowledge of the displacements at infinity. 
We shall show that these may not be specified arbitrarily. If we restrict 
()(z), w9(z) asin (32), then D,is bounded at infinity (see (1)) and the integrand 
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of J in (51) is o(z~*) on those parts of C,,,, which are to tend to infinity. 
Hence. if we may choose those parts of the contour so that their total length 
is O(z), their contribution to J is zero. The remainder of the proof is as 
before. The theorem may now be stated as follows: 

Let Q(z), w(z) give specified displacements over C,, Cy,..., C,, and over those 
parts of C,., in the finite part of the plane, and give specified terms of orders 


higher than those in (32) on those parts of C,,,, which are to tend to infinity. 


n 
Then the only additional potentials Q(z), we(z) which do not violate the above 
conditions, and are such that Q,(z), Qo(z), wo(z) are continuous in the closed 


9 


reqion and are of the orders specified in (32), are 
Q),.(z) X, w,(z) = Kaz+B, (60) 
where cand B are comple r constants. 


Thus in the statement of the problem we require at infinity only terms of 
orders higher than those given in (32) and, consequently, when solving the 
problems, need satisfy only the terms of these orders. 

Terms of lower degree must not be arbitrarily specified. This explains why 
i: result obtained by letting the boundary of the material tend to infinity 
ifter solving a finite boundary-value problem of given displacements is often 
satisfactory. Forexample, if we solve the problem of the annulus bounded 
by circles of radiia and b (a < b) with specified displacements over the inner 
boundary and zero displacements over the outer one, and then let b tend to 
infinity, all the coefficients in the potentials tend to zero. 


We now consider two particular limiting forms of C,, ,,. 


10. (1) The boundary-value problem when C, ,, tends to infinity in 

all directions 

Equations (36), (41) of the corresponding problem with boundary con- 
ditions of stress apply to the present case. Uniqueness of solution follows 
if B, F, A,, Ay, B,, By, p,, po are specified. This implies that all finite and 
infinite displacements at infinity are given, the displacements being con- 
sistent with zero resultant force on the material. Apart from p,, p. the above 
constants may be specified arbitrarily provided that they give a finite result- 
ant and p,, p, then take the values required for equilibrium. The remaining 


terms take unique values automatically. 


11. (2) The boundary-value problem when (,,, tends to the 
boundary of a half-plane 
Let C,, Cy,..., C, lie in the finite part of the plane and C, ,, consist of the 
real axis between the points x R and the semicircle, centred at the 


origin, radius R, above the real axis. Let the displacements be specified over 
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C,, Cy,..., C, and, in addition, let the displacements over y = 0 be given by 


U9 =wW=f(x), v= g(x), (61) 


where, for large |x 


f (x) A log\x|+ B+-Cargx+o(1), (62) 
d(x) Diog\x|+ E+ F argx+o(1), (63) 
A, Byouw, F are real constants, and arg 2 is 0 or 7 according as x is positive or 


negative. Let us suppose that, in addition, we require (’(z), w"(z) to be 
O(z-*) at infinity in y > 0. For large values of |z| we write 

P 1, +195 e qd, +14. ; 

Q’(z) = P24 o(2 1) w"(z) = BOO fe, o(z-), (64) 


~ 


so that 


Q(z) = (p,+ip,)logz+ M-+o(1), w'(z) = (q,+7tq,)logz+-N+0(1). (65) 
We have no assurance that the forms specified give a possible solution, for 
the plate may not then be in equilibrium. The solution is possible if (48) is 
satisfied. However, we will show that if a solution of the above form is 
known, then it is the unique solution to the problem. We show first that the 
constants ;, Po, 71, J2 are uniquely determined. From (65) the displace- 
ments at infinity are given by 
8uD = K(p,+tp,)log z—(p,—ips) = — (q,—tq2)log 7+-«M —N+0(1). 
Thus at infinity, on the real axis, 
8uD° = «(py +ip,)log |x| —(p,—ip,)— (q,—t2)log |x| + 
+-arg x|«(ip,—po)+ (iq, +42) |+«M—N-+0 (1). (66) 
From (65), (66), 
KPy—h =A, KAth=F,  q—Kpp=C,  qet+Kp, = D, 
so that 
2Kp, A+F, 2q, = F—A, 2K Do D—C, 242 C+D. 
Thus p,, Po, 9, J2 are given in terms of the displacements specified over 
y = 0. The forms chosen for the latter displacements include important 
special cases. If the specified displacements tend to zero at infinity, the 
constants );, Ps, J1, Jz are all zero and the displacements are everywhere 
bounded. The uniqueness theorem may now be completed as in the finite 
case. However, in the important case of the half-plane with no interior 
contours, a simpler integral may be employed and a shorter proof given. 
Consider the integral 
K = [ Q4(2)[8-D,(Z,2)+(e—2)Q4(2)] dz, (67) 
G 


where the contour C consists of the real axis between the points a LR 
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nd the semicircle above the real axis. The integrand refers to the additional 
notentials defined above. Clearly the integrand vanishes on the real axis 
nd is o (z~1) on the semicircle. Thus 


lim K 0. (68) 


fransforming to a surface integral by (19) 





K = 2 | Q4(z)] 8 = Diz.) Qi(z)| dS 
2Ki | [Q4(z)Q6(Z)] dS. 
Thus K = 2x«i | |Q4(z)|? dS. (69) 


Ss 
[he above transformation is justified by uniformity and continuity of Q5(z) 


ind D,(z,Z) in the region y 0. Thus we have 


Q(z) = 0 and Q,(z) v, a complex constant, in y > 0 (70) 
and 8uD, Ka @o(Z). 
But D, 0 on ¥ 0 so that. on y 0, we have 
@,(Z) KM, (71) 


nd by Cauchy’s integral theorem, (71) holds throughout y > 0. Thus 
D, = 0 throughout y > 0. 


12. Sections of C,,,, at infinity of finite total length 


Let the additional potentials Q,(z), w»(z) give zero displacements over all 


n 


boundaries except those at infinity where they satisfy the conditions 

(2, (z) O(z-*), w(z) = O(z 1), (72) 
Onthose parts of C,,,, which are to tend to infinity, Dis O(log z), so that the 
integrand of J in (51) is O(z-!logz) and these sections of C,,,, give zero 
imiting contribution to J. The remainder of the proof is similar to the 
finite case. We note that a special discussion of the terms in 2’(z), w”(z) 
which are O(z-1) is not necessary when C,,,, has only a finite section at 


infinity. 
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ON THE COMPRESSION OF A CUBE BETWEEN 
RIGID ROUGH PLATES 
By J. M. PRENTIS (Imperial College, London) 
teceived 13 March 1951] 


SUMMARY 
[he problem is to determine approximately the mode of deformation of a cube 
elastic material compressed between rigid plates sufficiently rough to prevent any 
tive movement between the plates and the ends of the specimen. A method of 
ving this type of problem has been given by Prager and Synge (1) and is applied 


to determine the relationship between the true Young’s modulus and the 
parent modulus obtained by a test of the type described above. As two other 
ers have recently been published (2, 3) applying the same method to similar 


blems, the work is given in outline only. 


1. Introduction 

lnE method of analysis evolved by Prager and Synge for the approximate 
solution of elastic boundary solution problems requires the selection of 
1 number of comparatively simple stress states which, in the aggregate, 
ipproximate to the natural stress conditions. The method of selecting the 
tificial states separates them into two classes, called the associated and 
the complementary. In one class the restrictions imposed by the equili- 
brium equations and the boundary conditions on stress are relaxed, whilst 
in the other class the compatibility equations and the imposed boundary 
onditions on displacements are ignored. Consideration of these two 
classes enables upper and lower bounds to be placed upon the total strain 
energy of the body. Ifa sufficient number of approximate states are taken 
n each class these bounds can be narrowed, allowing the total strain 


energy to be calculated to any desired degree of accuracy. 


2. Following the method used by Edelman (3), the problem Py (see 
Fig. 1 for details) may be resolved into two components. We shall con- 
sider P, in which the cube is subjected to a uni-axial compression, the 
lateral expansion being unrestricted, and P such that P+P, = P,. 

3. Symmetry allows us to confine our attention to one-eighth of the 
body, 0 < a, y <b, and 0 < z <h,b and h being the breadth and height 
of the specimen respectively 


The boundary conditions ol P are: 


u(O, y, 2) 0 U(X, Y, 2 h) oax ) 
u(x, 0,2) = 0 u, (x,y,z = h) cay } for displacements, (1) 
u(x.y, 0) 0 U(X, Y,z h) 0 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 2 (1952)] 
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o being Poisson’s ratio, and 
oS 


Dew, Deg = O for ze = 6 or 0 
Pry Py: = 9fory=borod|] . 
‘intial wey for stresses. (2 
o., = 0 for z = 6 
Py, = O for y = b 
FP 
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Fic. 1. 


4. This is a displacement boundary condition problem and _ for its 
solution Prager and Synge require a series of states: 

(a) The completely associated state S*, the deformations of which must 

satisfy equations (1) and the strains must satisfy the compatibility 


equations, 9 ; 
UW; s+ U; 2 5 j, (3) 


where uw; ; stands for éu;/éu;. 

(6) Asequence of states S$}, S),..., Sj,,, called the homogeneous associated 
states which must also satisfy equations (3) whilst the w’ are required 
to vanish wherever they are prescribed on the boundary by equa- 
tions (1). 

(c) A sequence of states S{, S3,..., S%, called the complementary states 
selected to satisfy equations (2) and the equilibrium equations, 


Pisa QO. (4) 


Except that they must satisfy the required conditions given above, the 
states are all quite arbitrary. 


5. Defining the scalar product of two states as 


(S.S’) = [ e,, pi; de, 





in 


fi 
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Prager and Synge show that 


n m 
s% 69”)2 « Y2 « * §* > (S* YT’ /2 5 
> (S*.1,) S (S*.S*)— > (S*.1,)%, (5) 
q=1 p=1 
Pp 
: where I’ (and I”) are orthonormal sequences S c,S}. and the c, are 


hn CT 
r=1 


chosen to satisfy the conditions, 


(I, 1.) = 8, 


§.. being the Kronecker delta, and where S? = (S.S) 2xthe total 


strain energy of the state P. 





__ 6. Selection of the states 
It is necessary to specify Poisson’s ratio and the relative dimensions of 
the specimen. We have taken o Land b =h l. 
& The associated states selected are given in Table I, where for convenience 
in computation # and a are temporarily put equal to unity. The strains 
m TABLE | 
The associated states 
u Uy wu 
S 0 0 (1 4) 
= S r2*(1—2* yz*(1 —z*) 0 
or Its S 0 0 2(1 —2z*)(a+y) 
S/ 1 l 4 y222(1 4) 0 
S 0 0 2(1 —z*)(x? + y?) 
must Ox r24/3 yz4/3 0 
bilit; . : i 
9, § corresponding to each state are obtained from equations (3) and the 
io 
stresses from 3 R 3 
Pij t\©xr Tt Cuy + €22)055 I 4Cij- 
iated The complementary states are most easily selected by use of Maxwell’s 
ate stress functions us, (7 x, y, 2), the stresses being obtained from these 
functions as follows: 
qua Ob, Cubs, 
Pra > T ¢ = etc., 
Cz" oy” 
ates 9 
tat C =. 
anc i. » etc. 
oxoy 
(4 Thus defined, the stresses automatically satisfy the equilibrium equations. 
the lhe only restrictions placed on the choice of the #’s is that equations (2) 


must be satisfied. The strains are obtained from the stresses by using 


Ci 303 — 3\ Pra Pyyt Paz)ij- 
In selecting the complementary states, deformations are not considered, 
so that the compatibility equations are not, in general satisfied, and 
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no corresponding displacements exist. The functions used are given in® 
Table IT. 
TaBLe II 


The complementary states 


ub P 


0 
0 
0 
(1 —a*)?(1 —y*)*(1 —2*) 


7. Computation of the apparent 


Using the selected functions it is found that 


4 
> (S*.1/)? = 0-060697a2E = L, say 
1 


q 
and (S*.S*)— > (S*.1)? = 0-063000a?E = U, say. 
p=1 
The strain energies of P, and P may be superimposed by simple addition 
so that S? S248? It is easily seen that S? 1 Fa and that S? 
Substitution in equation (5) gives 
L+a?7k < 8: 
Putting }F, the average stress, equal to f and, since h = 1, a (= a/h) =¢, 
the mean e,,, we obtain 
0-94278f/e > BE > 0-94073f/e, 
or, averaging, E = 0-94176f/e with a maximum possible error of +01 
per cent. 
Thus, in a test of this type the true Young’s modulus is 94-2 per cent. 
of the apparent modulus as given by the overall relative movement of the 


plates. 
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